Engineering topics

Petrosian, Vahe

Publications and source records attributed to Petrosian, Vahe.

At least 19 records

DOE OSTI2024

Inferring the Redshift of More than 150 GRBs with a Machine-learning Ensemble Model

Gamma-ray bursts (GRBs), due to their high luminosities, are detected up to a redshift of 10, and thus have the potential to be vital cosmological probes of early processes in the Universe. Fulfilling this potential requires a large sample of GRBs with known redshifts, but due to observational limitations, only 11% have known redshifts (z). There have been numerous attempts to estimate redshifts via correlation studies, most of which have led to inaccurate predictions. To overcome this, we estimated GRB redshift via an ensemble-supervised machine-learning (ML) model that uses X-ray afterglows of long-duration GRBs observed by the Neil Gehrels Swift Observatory. The estimated redshifts are strongly correlated (a Pearson coefficient of 0.93) and have an rms error, namely, the square root of the average squared error $\langle$Δz 2 $\rangle$, of 0.46 with the observed redshifts showing the reliability of this method. The addition of GRB afterglow parameters improves the predictions considerably by 63% compared to previous results in peer-reviewed literature. Finally, we use our ML model to infer the redshifts of 154 GRBs, which increase the known redshifts of long GRBs with plateaus by 94%, a significant milestone for enhancing GRB population studies that require large samples with redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2016

Data-Driven Radiative Hydrodynamic Modeling of the 2014 March 29 X1.0 Solar Flare

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in Mg II hk, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DSTIBIS) in H(alpha) 6563A and Ca II 8542A, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multithreaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated Mg II line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of approx. 800 km below the transition region can turn the simulated Mg II line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

Da Costa, Fatima Rubio↗

NASA NTRS2009

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗

NASA NTRS2005

Magnetic Variations Associated With Solar Flares

A report summarizes an investigation of helioseismic waves and magnetic variations associated with solar flares, involving analysis of data acquired by the Michelson Doppler Imager (MDI) aboard the Solar and Heliocentric Observatory (SOHO) spacecraft, the Yohkoh spacecraft, and the Ramaty High Energy Solar Spectroscopic Imager (RHESSI) spacecraft. Reconstruction of x-ray flare images from RHESSI data and comparison of them with MDI magnetic maps were performed in an attempt to infer the changes in the geometry of the magnetic field. It was established that in most flares observed with MDI, downward propagating shocks were much weaker than was one observed in the July 9, 1996 flare, which caused a strong helioseismic response. It was concluded that most of the observed impulsive variations result from direct impact of high-energy particles. Computer codes were developed for further study of these phenomena.

Petrosian, Vahe↗

NASA NTRS2005

Studies of Particle Acceleration, Transport and Radiation in Impulsive Phase of Solar Flares

Solar activity and its most prominent aspect, the solar flares, have considerable influence on terrestrial and space weather. Solar flares also provide a suitable laboratory for the investigation of many plasma and high energy processes important in the magnetosphere of the Earth and many other space and astrophysical situations. Hence, progress in understanding of flares will have considerable scientific and societal impact. The primary goal of this grant is the understanding of two of the most important problems of solar flare physics, namely the determination of the energy release mechanism and how this energy accelerates particles. This is done through comparison of the observations with theoretical models, starting from observations and gradually proceeding to theoretically more complex situations as the lower foundations of our understanding are secured. It is generally agreed that the source of the flare energy is the annihilation of magnetic fields by the reconnection process. Exactly how this energy is released or how it is dissipated remains controversial. Moreover, the exact mechanism of the acceleration of the particles is still a matter of debate. Data from many spacecrafts and ground based instruments obtained over the past decades have given us some clues. Theoretical analyses of these data have led to the standard thick target model (STT) where most of the released energy goes into an (assumed) power law spectrum of accelerated particles, and where all the observed radiations are the consequence of the interaction of these particles with the flare plasma. However, some theoretical arguments, and more importantly some new observations, have led us to believe that the above picture is not complete. It appears that plasma turbulence plays a more prominent role than suspected previously, and that it is the most likely agent for accelerating particles. The model we have developed is based on production of a high level of plasma waves and turbulence in the reconnection region above a flare loop. This turbulence accelerates particles stochastically which radiate some of their energy in this region but carry most of their energy to the footpoints of the loop, where they lose all their energy and radiate bulk of the observed radiation as in the traditional thick target model. In the past we have worked on various aspects of this model. We have evaluated the interaction rates of the plasma waves with electrons and ions, developed theoretical frameworks for the acceleration, transport and radiative processes, and produced numerical codes for the investigation of these processes. The goal of this grant has been further development and testing of this new paradigm, with emphases on the relative acceleration of electrons and ions and on a comprehensive investigation of the turbulence generation, cascade, and damping processes. We review several pieces of important evidence that we have uncovered indicating the crucial roles of turbulence, in and we describe accomplishments during the past two years of this grant.

Petrosian, Vahe↗

NASA NTRS2004

Analysis of RXTE data on Clusters of Galaxies

This grant provided support for the reduction, analysis and interpretation of of hard X-ray (HXR, for short) observations of the cluster of galaxies RXJO658--5557 scheduled for the week of August 23, 2002 under the RXTE Cycle 7 program (PI Vahe Petrosian, Obs. ID 70165). The goal of the observation was to search for and characterize the shape of the HXR component beyond the well established thermal soft X-ray (SXR) component. Such hard components have been detected in several nearby clusters. distant cluster would provide information on the characteristics of this radiation at a different epoch in the evolution of the imiverse and shed light on its origin. We (Petrosian, 2001) have argued that thermal bremsstrahlung, as proposed earlier, cannot be the mechanism for the production of the HXRs and that the most likely mechanism is Compton upscattering of the cosmic microwave radiation by relativistic electrons which are known to be present in the clusters and be responsible for the observed radio emission. Based on this picture we estimated that this cluster, in spite of its relatively large distance, will have HXR signal comparable to the other nearby ones. The planned observation of a relatively The proposed RXTE observations were carried out and the data have been analyzed. We detect a hard X-ray tail in the spectrum of this cluster with a flux very nearly equal to our predicted value. This has strengthen the case for the Compton scattering model. We intend the data obtained via this observation to be a part of a larger data set. We have identified other clusters of galaxies (in archival RXTE and other instrument data sets) with sufficiently high quality data where we can search for and measure (or at least put meaningful limits) on the strength of the hard component. With these studies we expect to clarify the mechanism for acceleration of particles in the intercluster medium and provide guidance for future observations of this intriguing phenomenon by instrument on GLAST. The details of the nonthermal particle population has important implications for the theories of cluster formation, mergers and evolution. The results of this work were first presented at the High Energy Division meeting of the American astronomical Society at Mt. Tremblene, Canada (Petrosian et al. 2003). and in an invited review talk at the General Assembly of the International Astronomical Union at Sydney, Australia (Petrosian, 2003). A paper describe the observations, the data analysis and its implication is being prepared for publication in the Astrophysical Journal.

Petrosian, Vahe↗

NASA NTRS2001

Publicly Available Numerical Codes for Modeling the X-ray and Microwave Emissions from Solar and Stellar Activity

We have posted numerical codes on the Web for modeling the bremsstrahlung x-ray emission and the a gyrosynchrotron radio emission from solar and stellar activity. In addition to radiation codes, steady-state and time-dependent Fokker-Planck codes are provided for computing the distribution and evolution of accelerated electrons. A 1-D hydrodynamics code computes the response of the stellar atmosphere (chromospheric evaporation). A code for modeling gamma-ray line spectra is also available. On-line documentation is provided for each code. These codes have been developed for modeling results from the High Energy Solar Spectroscopic Imager (HESSI) along related microwave observations of solar flares. Comprehensive codes for modeling images and spectra of solar flares are under development. The posted codes can be obtained on NASA/Goddard's HESSI Web Site at http://hesperia.gsfc.nasa.gov/hessi/modelware.htm. This work is supported in part by the NASA Sun-Earth Connection Program.

Holman, Gordon D.↗

NASA NTRS1995

On the correlation of angular position with time of occurrence of gamma-ray bursts

Evidence indicating that a large fraction of gamma-ray bursts are repeaters would provide strong support for noncosmological origin of these sources. Wang & Lingenfelter have claimed existance of a correlation between angular position and time of occurrence of bursts. We perform statistical tests and find a marginal evidence for nearby bursts occurring within 4 to 5 days of each other in the BATSE 1B catalog. This evidence is present also in the 2B catalogs, which in addition, shows some marginal evidence for bursts repetition at longer time delays up to the total length of the observations.

Petrosian, Vahe↗

NASA NTRS1994

Gamma-ray emission and electron acceleration in solar flares

Recent observations have extended the spectra of the impulsive phase of flares to the GeV range. Such high-energy photons can be produced either by electron bremsstrahlung or by decay of pions produced by accelerated protons. In this paper we investigate the effects of processes which become important at high energies. We examine the effects of synchrotron losses during the transport of electrons as they travel from the acceleration region in the corona to the gamma-ray emission sites deep in the chromosphere and photosphere, and the effects of scattering and absorption of gamma rays on their way from the photosphere to space instruments. These results are compared with the spectra from so-called electron-dominated flares, observed by GRS on the Solar Maximum Mission, which show negligible or no detectable contribution from accelerated protons. The spectra of these flares show a distinct steepening at energies below 100 keV and a rapid falloff at energies above 50 MeV. Following our earlier results based on lower energy gamma-ray flare emission we have modeled these spectra. We show that neither the radiative transfer effects, which are expected to become important at higher energies, nor the transport effects (Coulomb collisions, synchrotron losses, or magnetic field convergence) can explain such sharp spectral deviations from a simple power law. These spectral deviations from a power law are therefore attributed to the acceleration process. In a stochastic acceleration model the low-energy steepening can be attributed to Coulomb collision and the rapid high-energy steepening can result from synchrotron losses during the acceleration process.

Petrosian, Vahe↗

NASA NTRS1994

The dynamics of charged particles in turbulent astrophysical plasmas

We consider the resonant interaction of energetic charged particles and transverse plasma wave propagating parallel and/or antiparallel to the uniform magnetic field B(sub 0) in an underlying background plasma of density n. The coupling of the plasma waves and the energetic particles will be controlled by the ratio n/(the absolute value of B(sub 0)(exp 2). A variation of this ratio leads to a strong variation of the dynamics of the energetic particles. By taking into account the whole transverse plasma branch for the resonant interaction we discuss the influence of the background plasma density, the background magnetic field, the cross helicity, and the magnetic helicities on the dynamics of charged particles in astrophysical plasmas. It is shown that low-energy electrons can be accelerated efficiently by the higher electromagnetic waves and short-wavelength whistlers for low values of the ratio n/(the absolute value of B(sub 0)(exp 2), which means for low values of the ratio of plasma frequency to gyrofrequency.

Dung, Rudiger↗

NASA NTRS1994

Do gamma-ray bursts originate from an extended Galactic Halo of high-velocity neutron stars?

The gamma-ray burst brightness distribution is inhomogeneous and the distribution on the sky is nearly isotropic. These features argue against an association of gamma-ray bursts with those Galactic objects that are known to exhibit a strong concentration toward the Galactic center or plane. The observed statistical properties indicate a cosmological origin. Circumstantial evidence suggests that neutron stars are involved in the burst phenomenon. Here we consider Population II neutron stars in an extended Galactic Halo (EGH) as an alternative to cosmological scenarios. The Burst and Transient Source Experiment (BATSE) data indicate a small deviation from isotropy near the 2 sigma level of statistical significance. If confirmed for an increasing number of bursts, these anisotropies could rule out cosmological scenarios. On the other hand, EGH models require small anisotropies like those observed by BATSE. We consider simple distribution models to determine the generic properties such halos must have to be consistent with the observations and discuss the implications of the corresponding distance scale on burst models.

Hartmann, Dieter H.↗

NASA NTRS1993

Statistics of arcs in clusters of galaxies

Samples of gravitational lens events in clusters show many large arcs compared to arclets, relative to what can be obtained by idealized singular lens models. We describe the probability of image magnification for point sources and for simple but more realistic gravitational lensing models that include a finite core size and an ellipticity. In addition, we explore the changes in the probability distribution of image magnifications, distortions, and angular extents for sources of different sizes as the parameters of the lenses are varied. A finite core in spherically symmetric lens models introduces a discontinuity in the probability distribution at which the relative number of highly magnified images is increased. In elliptical lenses, this discontinuity and its effect are replaced by a continuous increase in the probability of obtaining high-magnification images relative to singular spherically symmetric models. We also find that the finite size of the source causes a further increase in the expected number of images just below the maximum possible magnification.

Bergmann, Anton G.↗

NASA NTRS1993

The distribution of flare parameters and implications for coronal heating

Assuming that solar hard X-ray flares are caused by bremsstrahlung from a distribution of accelerated electrons impinging on a thick target, we determine the relationship between the total flare energy released in the corona and observable flare parameters. It is shown that for the more numerous smaller flares to be energetically more important than larger ones, there must be an anticorrelation between observed flux and one or more of the photon spectral index, the total flare duration, the ratio of the total energy that goes into heating the corona to the energy that goes into acceleration of electrons, and the reciprocal of the low-energy cutoff of the accelerated electron spectrum. Using a new statistical test which properly takes account of data truncation due to observational selection effects, we find that neither spectral index nor duration show an anticorrelation with flux in ISEE 3/ICE hard X-ray data. Using a nonparametric method, we obtain single distributions of peak flux, spectral index, duration, and fluence, corrected for selection biases. The functional forms of these distributions do not support coronal heating by numerous small flares without very specific behavior by yet unobserved parameters.

Lee, Theodore T.↗

NASA NTRS1993

Interpretation of gamma-ray burst source count statistics

Ever since the discovery of gamma-ray bursts, the so-called log N-log S relation has been used for determination of their distances and distribution. This task has not been straightforward because of varying thresholds for the detection of bursts. Most of the current analyses of these data are couched in terms of ambiguous distributions, such as the distribution of Cp/Clim, the ratio of peak to threshold photon count rates, or the distribution of V/Vmax = (Cp/Clim) exp -3/2. It is shown that these distributions are not always a true reflection of the log N-log S relation. Some kind of deconvolution is required for obtaining the true log N-log S. Therefore, care is required in the interpretation of results of such analyses. A new method of analysis of these data is described, whereby the bivariate distribution of Cp and Clim is obtained directly from the data.

Petrosian, Vahe↗

NASA NTRS1992

A simple test of independence for truncated data with applications to redshift surveys

This paper presents an easily applied permutation test for H0, closely related to Lyden-Bell's (1971) estimate of the marginal distribution of truncated data. The test is applied to two redshift-magnitude surveys, one of galaxies and one of quasars. Assuming statistical independence, testing H0 amounts to testing validity of the cosmological model. Segal's (1986) chronomatic cosmological model is rejected under H0. On the other hand, for the quasar sample H0 is rejected strongly in a conventional cosmological model (and in a chronomatic model as well) indicating either incorrectness of the models or, as is more commonly assumed, indicating strong luminosity evolution.

Efron, Bradley↗

NASA NTRS1992

Stochastic acceleration of electrons. I - Effects of collisions in solar flares

Stochastic acceleration of thermal electrons to nonrelativistic energies is studied under solar flare conditions. We show that, in turbulent regions, electron-whistler wave interactions can result in the acceleration of electrons in times comparable to or shorter than the Coulomb collision time. The kinetic equation describing the evolution of the electron energy distribution including stochastic acceleration by whistlers and energy loss via Coulomb interactions is solved for an initial thermal electron energy spectrum. In general, the shape of the resulting electron distributions are characterized by the energy E(c) where systematic energy gain by turbulence equals energy loss due to Coulomb collisions. For energies less than E(c), the spectra are steep (quasi-thermal) whereas above E(c), the spectra are power laws. We find that hard X-ray spectra computed using the electron distributions obtained from our numerical simulations are able to explain the complex spectral shapes and variations observed in impulsive hard X-ray bursts. In particular, we show that the gradual steepening observed by Lin et al. (1981) could be due to a systematic increase in the density of the plasma (due to evaporation) and the increasing importance of collisions instead of the appearance of a superhot thermal component.

Hamilton, Russell J.↗

NASA NTRS1992

On the continuum radiation of gamma-ray bursts

I describe the qualitative features of a model for gamma-ray bursts assuming that they are produced by galactic neutron stars. The model is based on the assumption that electrons are continuously accelerated high up in the magnetosphere to relativistic energies and large pitch angles. These electrons emit a self-absorbed synchrotron spectrum at energies greater than mc-squared and, as a result, lose all of the energy associated with their motion perpendicular to the magnetic field. They thus slide along the field lines hitting the high density crust where by collision with ions acquire non-zero pitch angle (or are excited to higher Landau levels) and emit either cyclotron or bremsstrahlung photons.

Petrosian, Vahe↗

NASA NTRS1991

Evolution of the luminosity function of quasar accretion disks

Using an accretion-disk model, accretion disk luminosities are calculated for a grid of black hole masses and accretion rates. It is shown that, as the black-hole mass increases with time, the monochromatic luminosity at a given frequency first increases and then decreases rapidly as this frequency is crossed by the Wien cutoff. The upper limit on the monochromatic luminosity, which is characteristic for a given epoch, constrains the evolution of quasar luminosities and determines the evolultion of the quasar luminosity function.

Caditz, David M.↗