Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stochastic processes”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,783 records · Page 99

Filtering Solar-Like Oscillations for Exoplanet Detection in Radial Velocity Observations

Cool main-sequence, subgiant, and red giant stars all show solar-like oscillations, pulsations that are excited and intrinsically damped by near-surface convection. Many overtones are typically excited to observable amplitudes, giving a rich spectrum of detectable modes. These modes provide a wealth of information on fundamental stellar properties. However, the radial velocity (RV) shifts induced by these oscillations can also be problematic when searching for low-mass, long-period planets; this is because their amplitudes are large enough to completely mask such minute planetary signals. Here we show how fine-tuning exposure times to the stellar parameters can help efficiently average out the solar-like, oscillation-induced shifts. To reduce the oscillation signal to the RV precision commensurate with an Earth analog, we find that for cool, low-mass stars (near spectral type K), the necessary exposure times may be as short as ∼4 minutes, while for hotter, higher-mass stars (near spectral type F, or slightly evolved), the required exposure times can be longer than 100 minutes. We provide guideline exposure durations required to suppress the total observed amplitude due to oscillations to a level of 0.1 m s{sup −1}, and a level corresponding to the Earth-analog reflex amplitude for the star. Owing to the intrinsic stochastic variability of the oscillations, we recommend in practice choosing short exposure durations at the telescope and then averaging over those exposures later, as guided by our predictions. To summarize, as we enter an era of 0.1 m s{sup −1} instrumental precision, it is critical to tailor our observing strategies to the stellar properties.

79 ASTRONOMY AND ASTROPHYSICS↗

The Dynamics of Interstellar Asteroids and Comets within the Galaxy: An Assessment of Local Candidate Source Regions for 1I/’Oumuamua and 2I/Borisov

The low velocity of interstellar asteroid 1I/’Oumuamua with respect to our galaxy’s local standard of rest implies that it is young. Adopting the young age hypothesis, we assess possible origin systems for this interstellar asteroid and for 2I/Borisov, though the latter’s higher speed means it is unlikely to be young. First, their past trajectories are modeled under gravitational scattering by galactic components (“disk heating”) to assess how far back one can trace them. The stochastic nature of disk heating means that a back integration can only expect to be accurate to within 15 pc and 2 km s{sup −1} at −10 Myr, dropping steeply to 400 pc and 10 km s{sup −1} at −100 Myr, sharply limiting our ability to determine a precise origin. Nevertheless, we show that ’Oumuamua’s origin system likely is currently within 1 kpc of Earth, in the local Orion Arm. Second, we back integrate ’Oumuamua’s trajectory in order to assess source regions, emphasizing young systems and moving groups. Though disk heating allows for only a statistical link to source regions, ’Oumuamua passed through a considerable subset of the Carina and Columba moving groups when those groups were forming. This makes them perhaps the most plausible source region, if ’Oumuamua was ejected during planet formation or via intra-cluster interactions. We find three stars in the Ursa Major group, one brown dwarf, and seven other stars to have plausible encounters with 2I/Borisov, within 2 pc and 30 km s{sup −1}. These encounters’ high relative speeds mean none are likely to be the home of 2I/Borisov.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Accuracy of the ALMA Flux Calibration in the Time Domain and across Spectral Windows

A diverse array of science goals requires accurate flux calibration of observations with the Atacama Large Millimeter/submillimeter array (ALMA); however, this goal remains challenging due to the stochastic time-variability of the “grid” quasars ALMA uses for calibration. In this work, we use 343.5 GHz (Band 7) ALMA Atacama Compact Array observations of four bright and stable young stellar objects over seven epochs to independently assess the accuracy of the ALMA flux calibration and to refine the relative calibration across epochs. The use of these four extra calibrators allows us to achieve an unprecedented relative ALMA calibration accuracy of ~3%. On the other hand, when the observatory calibrator catalog is not up to date, the Band 7 data calibrated by the ALMA pipeline may have a flux calibration poorer than the nominal 10%, which can be exacerbated by weather-related phase decorrelation when self-calibration of the science target is either not possible or not attempted. We also uncover a relative flux calibration uncertainty between spectral windows of 0.8%, implying that measuring spectral indices within a single ALMA band is likely highly uncertain. We thus recommend various methods for science goals requiring high flux accuracy and robust calibration, in particular, the observation of additional calibrators combined with a relative calibration strategy, and observation of solar system objects for high absolute accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗

SCExAO/MEC and CHARIS Discovery of a Low-mass, 6 au Separation Companion to HIP 109427 Using Stochastic Speckle Discrimination and High-contrast Spectroscopy

We report the direct imaging discovery of a low-mass companion to the nearby accelerating A star, HIP 109427, with the Subaru Coronagraphic Extreme Adaptive Optics (SCExAO) instrument coupled with the Microwave Kinetic Inductance Detector Exoplanet Camera (MEC) and CHARIS integral field spectrograph. CHARIS data reduced with reference star point spread function (PSF) subtraction yield 1.1–2.4 μm spectra. MEC reveals the companion in Y and J band at a comparable signal-to-noise ratio using stochastic speckle discrimination, with no PSF subtraction techniques. Combined with complementary follow-up L {sub p} photometry from Keck/NIRC2, the SCExAO data favors a spectral type, effective temperature, and luminosity of M4–M5.5, 3000–3200 K, and log{sub 10}(L/L{sub ⊙})=−2.28{sub −0.04}{sup +0.04}, respectively. Relative astrometry of HIP 109427 B from SCExAO/CHARIS and Keck/NIRC2, and complementary Gaia–Hipparcos absolute astrometry of the primary favor a semimajor axis of 6.55{sup +3.0} {sub −0.48} au, an eccentricity of 0.54{sub −0.15}{sup +0.28}, an inclination of 66.7{sub −14}{sup +8.5} degrees, and a dynamical mass of 0.280{sub −0.059}{sup +0.18} M {sub ⊙}. This work shows the potential for extreme AO systems to utilize speckle statistics in addition to widely used postprocessing methods to directly image faint companions to nearby stars near the telescope diffraction limit.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

LBT Reveals Large Dust Particles and a High Mass-loss Rate for K2-22 b

The disintegrating planet candidate K2-22 b shows periodic and stochastic transits best explained by an escaping debris cloud. However, the mechanism that creates the debris cloud is unknown. The grain size of the debris as well as its sublimation rate can be helpful in understanding the environment that disintegrates the planet. Here, we present simultaneous photometry with the g band at 0.48 μm and K {sub S} band at 2.1 μm using the Large Binocular Telescope. During an event with very low dust activity, we put a new upper limit on the size of the planet of 0.71 R {sub ⊕} or 4500 km. We also detected a medium depth transit that can be used to constrain the dust particle sizes. We find that the median particle size must be larger than about 0.5–1.0 μm, depending on the composition of the debris. This leads to a high mass-loss rate of about 3 × 10{sup 8} kg s{sup −1}, which is consistent with hydrodynamic escape models. If they are produced by some alternate mechanism such as explosive volcanism, it would require extraordinary geological activity. Combining our upper limits on the planet size with the high mass-loss rate, we find a lifetime of the planet of less than 370 Myr. This drops to just 21 Myr when adopting the 0.02 M {sub ⊕} mass predicted from hydrodynamical models.

47 OTHER INSTRUMENTATION↗

THE BURSTY STAR FORMATION HISTORIES OF LOW-MASS GALAXIES AT 0.4 < z < 1 REVEALED BY STAR FORMATION RATES MEASURED FROM H β AND FUV

We investigate the burstiness of star formation histories (SFHs) of galaxies at 0.4 < z < 1 by using the ratio of star formation rates (SFRs) measured from H β and FUV (1500 Å) (H β -to-FUV ratio). Our sample contains 164 galaxies down to stellar mass ( M {sub *}) of 10{sup 8.5} M {sub ⊙} in the CANDELS GOODS-N region, where Team Keck Redshift Survey Keck/DEIMOS spectroscopy and Hubble Space Telescope /WFC3 F275W images from CANDELS and Hubble Deep UV Legacy Survey are available. When the ratio of H β - and FUV-derived SFRs is measured, dust extinction correction is negligible (except for very dusty galaxies) with the Calzetti attenuation curve. The H β -to-FUV ratio of our sample increases with M {sub *} and SFR. The median ratio is ∼0.7 at M {sub *} ∼ 10{sup 8.5} M {sub ⊙} (or SFR ∼ 0.5 M {sub ⊙} yr{sup −1}) and increases to ∼1 at M {sub *} ∼ 10{sup 10} M {sub ⊙} (or SFR ∼ 10 M {sub ⊙} yr{sup −1}). At M {sub *} < 10{sup 9.5} M {sub ⊙}, our median H β -to-FUV ratio is lower than that of local galaxies at the same M {sub *}, implying a redshift evolution. Bursty SFH on a timescale of a few tens of megayears on galactic scales provides a plausible explanation for our results, and the importance of the burstiness increases as M {sub *} decreases. Due to sample selection effects, our H β -to-FUV ratio may be an upper limit of the true value of a complete sample, which strengthens our conclusions. Other models, e.g., non-universal initial mass function or stochastic star formation on star cluster scales, are unable to plausibly explain our results.

79 ASTRONOMY AND ASTROPHYSICS↗

THE ASTEROID BELT AS A RELIC FROM A CHAOTIC EARLY SOLAR SYSTEM

The orbital structure of the asteroid belt holds a record of the solar system’s dynamical history. The current belt only contains ∼10{sup −3} Earth masses yet the asteroids’ orbits are dynamically excited, with a large spread in eccentricity and inclination. In the context of models of terrestrial planet formation, the belt may have been excited by Jupiter’s orbital migration. The terrestrial planets can also be reproduced without invoking a migrating Jupiter; however, as it requires a severe mass deficit beyond Earth’s orbit, this model systematically under-excites the asteroid belt. Here we show that the orbits of the asteroids may have been excited to their current state if Jupiter’s and Saturn’s early orbits were chaotic. Stochastic variations in the gas giants’ orbits cause resonances to continually jump across the main belt and excite the asteroids’ orbits on a timescale of tens of millions of years. While hydrodynamical simulations show that the gas giants were likely in mean motion resonance at the end of the gaseous disk phase, small perturbations could have driven them into a chaotic but stable state. The gas giants’ current orbits were achieved later, during an instability in the outer solar system. Although it is well known that the present-day solar system exhibits chaotic behavior, our results suggest that the early solar system may also have been chaotic.

79 ASTRONOMY AND ASTROPHYSICS↗

BOOSTING LY α AND He ii λ 1640 LINE FLUXES FROM POPULATION III GALAXIES: STOCHASTIC IMF SAMPLING AND DEPARTURES FROM CASE-B

We revisit calculations of nebular hydrogen Ly α and He ii λ 1640 line strengths for Population III (Pop III) galaxies, undergoing continuous, and bursts of, star formation. We focus on initial mass functions (IMFs) motivated by recent theoretical studies, which generally span a lower range of stellar masses than earlier works. We also account for case-B departures and the stochastic sampling of the IMF. In agreement with previous work, we find that departures from case-B can enhance the Ly α flux by a factor of a few, but we argue that this enhancement is driven mainly by collisional excitation and ionization, and not due to photoionization from the n = 2 state of atomic hydrogen. The increased sensitivity of the Ly α flux to the high-energy end of the galaxy spectrum makes it more subject to stochastic sampling of the IMF. The latter introduces a dispersion in the predicted nebular line fluxes around the deterministic value by as much as a factor of ∼4. In contrast, the stochastic sampling of the IMF has less impact on the emerging Lyman Werner photon flux. When case-B departures and stochasticity effects are combined, nebular line emission from Pop III galaxies can be up to one order of magnitude brighter than predicted by “standard” calculations that do not include these effects. This enhances the prospects for detection with future facilities such as the James Webb Space Telescope and large, ground-based telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗

M31N 2008-12a—the remarkable recurrent nova in M31: panchromatic observations of the 2015 eruption

The Andromeda Galaxy recurrent nova M31N 2008-12a had been observed in eruption 10 times, including yearly eruptions from 2008 to 2014. With a measured recurrence period of P{sub rec}=351±13 days (we believe the true value to be half of this) and a white dwarf very close to the Chandrasekhar limit, M31N 2008-12a has become the leading pre-explosion supernova type Ia progenitor candidate. Following multi-wavelength follow-up observations of the 2013 and 2014 eruptions, we initiated a campaign to ensure early detection of the predicted 2015 eruption, which triggered ambitious ground- and space-based follow-up programs. In this paper we present the 2015 detection, visible to near-infrared photometry and visible spectroscopy, and ultraviolet and X-ray observations from the Swift observatory. The LCOGT 2 m (Hawaii) discovered the 2015 eruption, estimated to have commenced at August 28.28 ± 0.12 UT. The 2013–2015 eruptions are remarkably similar at all wavelengths. New early spectroscopic observations reveal short-lived emission from material with velocities ∼13,000 km s{sup −1}, possibly collimated outflows. Photometric and spectroscopic observations of the eruption provide strong evidence supporting a red giant donor. An apparently stochastic variability during the early supersoft X-ray phase was comparable in amplitude and duration to past eruptions, but the 2013 and 2015 eruptions show evidence of a brief flux dip during this phase. The multi-eruption Swift/XRT spectra show tentative evidence of high-ionization emission lines above a high-temperature continuum. Following Henze et al. (2015a), the updated recurrence period based on all known eruptions is P{sub rec}=174±10 days, and we expect the next eruption of M31N 2008-12a to occur around 2016 mid-September.

79 ASTRONOMY AND ASTROPHYSICS↗

REVISITING EVIDENCE OF CHAOS IN X-RAY LIGHT CURVES: THE CASE OF GRS 1915+105

Nonlinear time series analysis has been widely used to search for signatures of low-dimensional chaos in light curves emanating from astrophysical bodies. A particularly popular example is the microquasar GRS 1915+105, whose irregular but systematic X-ray variability has been well studied using data acquired by the Rossi X-ray Timing Explorer . With a view to building simpler models of X-ray variability, attempts have been made to classify the light curves of GRS 1915+105 as chaotic or stochastic. Contrary to some of the earlier suggestions, after careful analysis, we find no evidence for chaos or determinism in any of the GRS 1915+105 classes. The dearth of long and stationary data sets representing all the different variability classes of GRS 1915+105 makes it a poor candidate for analysis using nonlinear time series techniques. We conclude that either very exhaustive data analysis with sufficiently long and stationary light curves should be performed, keeping all the pitfalls of nonlinear time series analysis in mind, or alternative schemes of classifying the light curves should be adopted. The generic limitations of the techniques that we point out in the context of GRS 1915+105 affect all similar investigations of light curves from other astrophysical sources.

79 ASTRONOMY AND ASTROPHYSICS↗

The MOSDEF survey: AGN multi-wavelength identification, selection biases, and host galaxy properties

We present results from the MOSFIRE Deep Evolution Field (MOSDEF) survey on the identification, selection biases, and host galaxy properties of 55 X-ray, IR, and optically selected active galactic nuclei (AGNs) at 1.4<z<3.8. We obtain rest-frame optical spectra of galaxies and AGNs and use the BPT diagram to identify optical AGNs. We examine the uniqueness and overlap of the AGNs identified at different wavelengths. There is a strong bias against identifying AGNs at any wavelength in low-mass galaxies, and an additional bias against identifying IR AGNs in the most massive galaxies. AGN hosts span a wide range of star formation rates (SFRs), similar to inactive galaxies once stellar mass selection effects are accounted for. However, we find (at ∼2–3σ significance) that IR AGNs are in less dusty galaxies with relatively higher SFR and optical AGNs in dusty galaxies with relatively lower SFR. X-ray AGN selection does not display a bias with host galaxy SFR. These results are consistent with those from larger studies at lower redshifts. Within star-forming galaxies, once selection biases are accounted for, we find AGNs in galaxies with similar physical properties as inactive galaxies, with no evidence for AGN activity in particular types of galaxies. This is consistent with AGNs being fueled stochastically in any star-forming host galaxy. We do not detect a significant correlation between SFR and AGN luminosity for individual AGN hosts, which may indicate the timescale difference between the growth of galaxies and their supermassive black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

Stochastic Fermi Energization of Coronal Plasma during Explosive Magnetic Energy Release

The aim of this study is to analyze the interaction of charged particles (ions and electrons) with randomly formed particle scatterers (e.g., large-scale local “magnetic fluctuations” or “coherent magnetic irregularities”) using the setup proposed initially by Fermi. These scatterers are formed by the explosive magnetic energy release and propagate with the Alfvén speed along the irregular magnetic fields. They are large-scale local fluctuations ( δB / B ≈ 1) randomly distributed inside the unstable magnetic topology and will here be called Alfvénic Scatterers (AS). We constructed a 3D grid on which a small fraction of randomly chosen grid points are acting as AS. In particular, we study how a large number of test particles evolves inside a collection of AS, analyzing the evolution of their energy distribution and their escape-time distribution. We use a well-established method to estimate the transport coefficients directly from the trajectories of the particles. Using the estimated transport coefficients and solving the Fokker–Planck equation numerically, we can recover the energy distribution of the particles. We have shown that the stochastic Fermi energization of mildly relativistic and relativistic plasma can heat and accelerate the tail of the ambient particle distribution as predicted by Parker and Tidman and Ramaty. The temperature of the hot plasma and the tail of the energetic particles depend on the mean free path ( λ {sub sc}) of the particles between the scatterers inside the energization volume.

79 ASTRONOMY AND ASTROPHYSICS↗

The Stochastic X-Ray Variability of the Accreting Millisecond Pulsar MAXI J0911–655

In this work, I report on the stochastic X-ray variability of the 340 Hz accreting millisecond pulsar MAXI J0911–655. Analyzing pointed observations of the XMM-Newton and NuSTAR observatories, I find that the source shows broad band-limited stochastic variability in the 0.01--10 Hz range with a total fractional variability of ∼24% rms in the 0.4--3 keV energy band that increases to ∼40% rms in the 3–10 keV band. Additionally, a pair of harmonically related quasi-periodic oscillations (QPOs) are discovered. The fundamental frequency of this harmonic pair is observed between frequencies of 62 and 146 mHz. Like the band-limited noise, the amplitudes of the QPOs show a steep increase as a function of energy; this suggests that they share a similar origin, likely the inner accretion flow. Based on their energy dependence and frequency relation with respect to the noise terms, the QPOs are identified as low-frequency oscillations and discussed in terms of the Lense–Thirring precession model.

79 ASTRONOMY AND ASTROPHYSICS↗

Reconstruction of Galaxy Star Formation Histories through SED Fitting:The Dense Basis Approach

We introduce the dense basis method for Spectral Energy Distribution (SED) fitting. It accurately recovers traditional SED parameters, including M {sub *}, SFR, and dust attenuation, and reveals previously inaccessible information about the number and duration of star formation episodes and the timing of stellar mass assembly, as well as uncertainties in these quantities. This is done using basis star formation histories (SFHs) chosen by comparing the goodness-of-fit of mock galaxy SEDs to the goodness-of-reconstruction of their SFHs. We train and validate the method using a sample of realistic SFHs at z = 1 drawn from stochastic realizations, semi-analytic models, and a cosmological hydrodynamical galaxy formation simulation. The method is then applied to a sample of 1100 CANDELS GOODS-S galaxies at 1 10{sup 9.5} M{sub ⊙}. About 40% of the CANDELS galaxies have SFHs whose maximum occurs at or near the epoch of observation. The dense basis method is scalable and offers a general approach to a broad class of data-science problems.

79 ASTRONOMY AND ASTROPHYSICS↗

Shock Acceleration Model for the Toothbrush Radio Relic

Although many of the observed properties of giant radio relics detected in the outskirts of galaxy clusters can be explained by relativistic electrons accelerated at merger-driven shocks, significant puzzles remain. In the case of the so-called Toothbrush relic, the shock Mach number estimated from X-ray observations (M{sub X}≈1.2--1.5) is substantially weaker than that inferred from the radio spectral index (M{sub rad}≈2.8). Toward understanding such a discrepancy, we here consider the following diffusive shock acceleration (DSA) models: (1) weak-shock models with M{sub s}≲2 and a preexisting population of cosmic-ray electrons (CRe) with a flat energy spectrum, and (2) strong-shock models with M{sub s}≈3 and either shock-generated suprathermal electrons or preexisting fossil CRe. We calculate the synchrotron emission from the accelerated CRe, following the time evolution of the electron DSA, and the subsequent radiative cooling and postshock turbulent acceleration (TA). We find that both models could reproduce reasonably well the observed integrated radio spectrum of the Toothbrush relic, but the observed broad transverse profile requires the stochastic acceleration by downstream turbulence, which we label “turbulent acceleration” or TA to distinguish it from DSA. Moreover, to account for the almost uniform radio spectral index profile along the length of the relic, the weak-shock models require a preshock region over 400 kpc with a uniform population of preexisting CRe with a high cutoff energy (≳40 GeV). Due to the short cooling time, it is challenging to explain the origin of such energetic electrons. Therefore, we suggest the strong-shock models with low-energy seed CRe (≲150 MeV) are preferred for the radio observations of this relic.

79 ASTRONOMY AND ASTROPHYSICS↗

Facilitating Follow-up of LIGO–Virgo Events Using Rapid Sky Localization

We discuss an algorithm for accurate and very low-latency (<1 s) localization of gravitational-wave (GW) sources using only the relative times of arrival, relative phases, and relative signal-to-noise ratios for pairs of detectors. The algorithm is independent of distances and masses to leading order, and can be generalized to all discrete (as opposed to stochastic and continuous) sources detected by ground-based detector networks. Our approach is similar to that of BAYESTAR with a few modifications, which result in increased computational efficiency. For the LIGO two-detector configuration (Hanford+Livingston) operating in O1 we find a median 50% (90%) localization of 143 deg{sup 2} (558 deg{sup 2}) for binary neutron stars. We use our algorithm to explore the improvement in localization resulting from loud events, finding that the loudest out of the first 4 (or 10) events reduces the median sky-localization area by a factor of 1.9 (3.0) for the case of two GW detectors, and 2.2 (4.0) for three detectors. We also consider the case of multi-messenger joint detections in both the gravitational and the electromagnetic radiation, and show that joint localization can offer significant improvements (e.g., in the case of LIGO and Fermi /GBM joint detections). We show that a prior on the binary inclination, potentially arising from GRB observations, has a negligible effect on GW localization. Our algorithm is simple, fast, and accurate, and may be of particular utility in the development of multi-messenger astronomy.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiple Stellar Populations of Globular Clusters from Homogeneous Ca-CN Photometry. II. M5 (NGC 5904) and a New Filter System

Using our ingeniously designed new filter systems, we investigate multiple stellar populations of the red giant branch (RGB) and the asymptotic giant branch (AGB) in the globular cluster (GC) M5. Our results are the following. (1) Our cn{sub JWL} index accurately traces nitrogen abundances in M5, while other color indices fail to do so. (2) We find bimodal CN distributions in both RGB and AGB sequences, with number ratios between CN-weak (CN-w) and CN-strong (CN-s) of n(CN-w):n(CN-s) = 29:71 (±2) and 21:79 (±7), respectively. (3) We also find a bimodal photometric [N/Fe] distribution for M5 RGB stars. (4) Our cn{sub JWL}–[O/Fe] and cn{sub JWL}–[Na/Fe] relations show clear discontinuities between the two RGB populations. (5) Although small, the RGB bump of CN-s is slightly brighter, ΔV{sub bump} = 0.07 ± 0.04 mag. If real, the difference in the helium abundance becomes ΔY = 0.028 ± 0.016, in the sense that CN-s is more helium enhanced. (6) Very similar radial but different spatial distributions with comparable center positions are found for the two RGB populations. The CN-s RGB and AGB stars are more elongated along the NW–SE direction. (7) The CN-s population shows a substantial net projected rotation, while that of the CN-w population is nil. (8) Our results confirm the deficiency of CN-w AGB stars previously noted by others. We show that it is most likely due to stochastic truncation in the outer part of the cluster. Finally, we discuss the formation scenario of M5.

79 ASTRONOMY AND ASTROPHYSICS↗

Precipitation and Release of Solar Energetic Particles from the Solar Coronal Magnetic Field

Most solar energetic particles (SEPs) are produced in the corona. They propagate through complex coronal magnetic fields subject to scattering and diffusion across the averaged field lines by turbulence. We examine the behaviors of particle transport using a stochastic 3D focused transport simulation in a potential field source surface model of coronal magnetic field. The model is applied to an SEP event on 2010 February 7. We study three scenarios of particle injection at (i) the compact solar flare site, (ii) the coronal mass ejection (CME) shock, and (iii) the EUV wave near the surface. The majority of particles injected on open field lines are able to escape the corona. We found that none of our models can explain the observations of wide longitudinal SEP spread without perpendicular diffusion. If the perpendicular diffusion is about 10% of what is derived from the random walk of field lines at the rate of supergranular diffusion, particles injected at the compact solar flare site can spread to a wide range of longitude and latitude, very similar to the behavior of particles injected at a large CME shock. Stronger pitch-angle scattering results in a little more lateral spread by holding the particles in the corona for longer periods of time. Some injected particles eventually end up precipitating onto the solar surface. Even with a very small perpendicular diffusion, the pattern of the particle precipitation can be quite complicated depending on the detailed small-scale coronal magnetic field structures, which could be seen with future sensitive gamma-ray telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗