BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

L. Lin

L. Lin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEphysics.comp-phastro-ph.COastro-ph.GAastro-ph.SRphysics.plasm-phastro-ph.IMhep-latphysics.space-phastro-phComputational Engineering, Finance, and Sciencecond-mat.mtrl-scicond-mat.str-elphysics.flu-dyncond-mat.mes-hallmath.NANumerical Analysisq-fin.GNq-fin.RM

Catalog footprint

What is connected

32works
19topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

32 published item(s)

preprint2022arXiv

GOODS-ALMA 2.0: Starbursts in the main sequence reveal compact star formation regulating galaxy evolution prequenching

Compact star formation appears to be generally common in dusty star-forming galaxies (SFGs). However, its role in the framework set by the scaling relations in galaxy evolution remains to be understood. In this work we follow up on the galaxy sample from the GOODS-ALMA 2.0 survey, an ALMA blind survey at 1.1mm covering a continuous area of 72.42arcmin$^2$ using two array configurations. We derived physical properties, such as star formation rates, gas fractions, depletion timescales, and dust temperatures for the galaxy sample built from the survey. There exists a subset of galaxies that exhibit starburst-like short depletion timescales, but they are located within the scatter of the so-called main sequence of SFGs. These are dubbed starbursts in the main sequence and display the most compact star formation and they are characterized by the shortest depletion timescales, lowest gas fractions, and highest dust temperatures of the galaxy sample, compared to typical SFGs at the same stellar mass and redshift. They are also very massive, accounting for $\sim 60\%$ of the most massive galaxies in the sample ($\log (M_{\rm{*}}/M_{\odot}) > 11.0$). We find trends between the areas of the ongoing star formation regions and the derived physical properties for the sample, unveiling the role of compact star formation as a physical driver of these properties. Starbursts in the main sequence appear to be the extreme cases of these trends. We discuss possible scenarios of galaxy evolution to explain the results drawn from our galaxy sample. Our findings suggest that the star formation rate is sustained in SFGs by gas and star formation compression, keeping them within the main sequence even when their gas fractions are low and they are presumably on the way to quiescence.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Broadband X-ray Burst Spectroscopy of the FRB-Emitting Galactic Magnetar

Magnetars are young, magnetically-powered neutron stars possessing the strongest magnetic fields in the Universe. Fast Radio Bursts (FRBs) are extremely intense millisecond-long radio pulses of primarily extragalactic origin, and a leading attribution for their genesis focuses on magnetars. A hallmark signature of magnetars is their emission of bright, hard X-ray bursts of sub-second duration. On April 27th 2020, the Galactic magnetar SGR J1935+2154 emitted hundreds of X-ray bursts in a few hours. One of these temporally coincided with an FRB, the first detection of an FRB from the Milky Way. Here we present spectral and temporal analyses of 24 X-ray bursts emitted 13 hours prior to the FRB and seen simultaneously with the NASA NICER and Fermi/GBM missions in their combined energy range, 0.2 keV-30 MeV. These broadband spectra permit direct comparison with the spectrum of the FRB-associated X-ray burst (FRB-X). We demonstrate that all 24 NICER/GBM bursts are very similar temporally, albeit strikingly different spectrally, from FRB-X. The singularity of the FRB-X burst is perhaps indicative of an uncommon locale for its origin. We suggest that this event originated in quasi-polar open or closed magnetic field lines that extend to high altitudes.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

GOODS-ALMA 2.0: Source catalog, number counts, and prevailing compact sizes in 1.1 mm galaxies

Submillimeter/millimeter observations of dusty star-forming galaxies with the Atacama Large Millimeter/submillimeter Array (ALMA) have shown that dust continuum emission generally occurs in compact regions smaller than the stellar distribution. However, it remains to be understood how systematic these findings are. Studies often lack homogeneity in the sample selection, target discontinuous areas with inhomogeneous sensitivities, and suffer from modest $uv$ coverage coming from single array configurations. GOODS-ALMA is a 1.1mm galaxy survey over a continuous area of 72.42arcmin$^2$ at a homogeneous sensitivity. In this version 2.0, we present a new low resolution dataset and its combination with the previous high resolution dataset from the survey, improving the $uv$ coverage and sensitivity reaching an average of $σ= 68.4μ$Jy beam$^{-1}$. A total of 88 galaxies are detected in a blind search (compared to 35 in the high resolution dataset alone), 50% at $S/N_{peak} \geq 5$ and 50% at $3.5 \leq S/N_{peak} \leq 5$ aided by priors. Among them, 13 out of the 88 are optically dark or faint sources ($H$- or $K$-band dropouts). The sample dust continuum sizes at 1.1mm are generally compact, with a median effective radius of $R_{e} = 0"10 \pm 0"05$ (a physical size of $R_{e} = 0.73 \pm 0.29$kpc at the redshift of each source). Dust continuum sizes evolve with redshift and stellar mass resembling the trends of the stellar sizes measured at optical wavelengths, albeit a lower normalization compared to those of late-type galaxies. We conclude that for sources with flux densities $S_{1.1mm} > 1$mJy, compact dust continuum emission at 1.1mm prevails, and sizes as extended as typical star-forming stellar disks are rare. The $S_{1.1mm} < 1$mJy sources appear slightly more extended at 1.1mm, although they are still generally compact below the sizes of typical star-forming stellar disks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

JINGLE -- IV. Dust, HI gas and metal scaling laws in the local Universe

Scaling laws of dust, HI gas and metal mass with stellar mass, specific star formation rate and metallicity are crucial to our understanding of the buildup of galaxies through their enrichment with metals and dust. In this work, we analyse how the dust and metal content varies with specific gas mass ($M_{\text{HI}}$/$M_{\star}$) across a diverse sample of 423 nearby galaxies. The observed trends are interpreted with a set of Dust and Element evolUtion modelS (DEUS) - incluidng stellar dust production, grain growth, and dust destruction - within a Bayesian framework to enable a rigorous search of the multi-dimensional parameter space. We find that these scaling laws for galaxies with $-1.0\lesssim \log M_{\text{HI}}$/$M_{\star}\lesssim0$ can be reproduced using closed-box models with high fractions (37-89$\%$) of supernova dust surviving a reverse shock, relatively low grain growth efficiencies ($ε$=30-40), and long dus lifetimes (1-2\,Gyr). The models have present-day dust masses with similar contributions from stellar sources (50-80\,$\%$) and grain growth (20-50\,$\%$). Over the entire lifetime of these galaxies, the contribution from stardust ($>$90\,$\%$) outweighs the fraction of dust grown in the interstellar medium ($<$10$\%$). Our results provide an alternative for the chemical evolution models that require extremely low supernova dust production efficiencies and short grain growth timescales to reproduce local scaling laws, and could help solving the conundrum on whether or not grains can grow efficiently in the interstellar medium.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A gPAV-Based Unconditionally Energy-Stable Scheme for Incompressible Flows with Outflow/Open Boundaries

We present an unconditionally energy-stable scheme for approximating the incompressible Navier-Stokes equations on domains with outflow/open boundaries. The scheme combines the generalized Positive Auxiliary Variable (gPAV) approach and a rotational velocity-correction type strategy, and the adoption of the auxiliary variable simplifies the numerical treatment for the open boundary conditions. The discrete energy stability of the proposed scheme has been proven, irrespective of the time step sizes. Within each time step the scheme entails the computation of two velocity fields and two pressure fields, by solving an individual de-coupled Helmholtz (including Poisson) type equation with a constant pre-computable coefficient matrix for each of these field variables. The auxiliary variable, being a scalar number, is given by a well-defined explicit formula within a time step, which ensures the positivity of its computed values. Extensive numerical experiments with several flows involving outflow/open boundaries in regimes where the backflow instability becomes severe have been presented to test the performance of the proposed method and to demonstrate its stability at large time step sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

An Energy-Stable Scheme for Incompressible Navier-Stokes Equations with Periodically Updated Coefficient Matrix

We present an energy-stable scheme for simulating the incompressible Navier-Stokes equations based on the generalized Positive Auxiliary Variable (gPAV) framework. In the gPAV-reformulated system the original nonlinear term is replaced by a linear term plus a correction term, where the correction term is put under control by an auxiliary variable. The proposed scheme incorporates a pressure-correction type strategy into the gPAV procedure, and it satisfies a discrete energy stability property. The scheme entails the computation of two copies of the velocity and pressure within a time step, by solving an individual de-coupled linear equation for each of these field variables. Upon discretization the pressure linear system involves a constant coefficient matrix that can be pre-computed, while the velocity linear system involves a coefficient matrix that is updated periodically, once every $k_0$ time steps in the current work, where $k_0$ is a user-specified integer. The auxiliary variable, being a scalar-valued number, is computed by a well-defined explicit formula, which guarantees the positivity of its computed values. It is observed that the current method can produce accurate simulation results at large (or fairly large) time step sizes for the incompressible Navier-Stokes equations. The impact of the periodic coefficient-matrix update on the overall cost of the method is observed to be small in typical numerical simulations. Several flow problems have been simulated to demonstrate the accuracy and performance of the method developed herein.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Fast optical absorption spectra calculations for periodic solid state systems

We present a method to construct an efficient approximation to the bare exchange and screened direct interaction kernels of the Bethe-Salpeter Hamiltonian for periodic solid state systems via the interpolative separable density fitting technique. We show that the cost of constructing the approximate Bethe-Salpeter Hamiltonian scales nearly optimally as $\mathcal{O}(N_k)$ with respect to the number of samples in the Brillouin zone $N_k$. In addition, we show that the cost for applying the Bethe-Salpeter Hamiltonian to a vector scales as $\mathcal{O}(N_k \log N_k)$. Therefore the optical absorption spectrum, as well as selected excitation energies can be efficiently computed via iterative methods such as the Lanczos method. This is a significant reduction from the $\mathcal{O}(N_k^2)$ and $\mathcal{O}(N_k^3)$ scaling associated with a brute force approach for constructing the Hamiltonian and diagonalizing the Hamiltonian respectively. We demonstrate the efficiency and accuracy of this approach with both one-dimensional model problems and three-dimensional real materials (graphene and diamond). For the diamond system with $N_k=2197$, it takes $6$ hours to assemble the Bethe-Salpeter Hamiltonian and $4$ hours to fully diagonalize the Hamiltonian using $169$ cores when the brute force approach is used. The new method takes less than $3$ minutes to set up the Hamiltonian and $24$ minutes to compute the absorption spectrum on a single core.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electron tunneling through a single magnetic barrier in HgTe topological insulator

Electron tunneling through a single magnetic barrier in a HgTe topological insulator has been theoretically investigated. We find that the perpendicular magnetic field would not lead to spin-flip of the edge states due to the conservation of the angular moment. By tuning the magnetic field and Fermi energy, the edge channels can be transited from switch-on states to switch-off states and the current can be transmitted from unpolarized states to totally spin polarized states. These features offer us and efficient way to control the topological edge state transport, and pave a way to construct the nanoelectronic devices utilizing the topological edge states.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

eXTP -- enhanced X-ray Timing and Polarimetry Mission

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hexagonal phase stabilization and magnetic orders of multiferroic Lu$_{1-x}$Sc$_x$FeO$_3$

Hexagonal LuFeO$_3$ has drawn a lot of research attention due to its contentious room-temperature multiferroicity. Due to the unstability of hexagonal phase in the bulk form, most experimental studies focused on LuFeO$_3$ thin films which can be stabilized by strain using proper substrates. Here we report on the hexagonal phase stabilization, magnetism, and magnetoelectric coupling of bulk LuFeO$_3$ by partial Sc-substitution of Lu. First, our first-principles calculations show that the hexagonal structure can be stabilized by partial Sc substitution, while the multiferroic properties including the noncollinear magnetic order and geometric ferroelectricity remain robustly unaffected. Therefore, Lu$_{1-x}$Sc$_x$FeO$_3$ can act as a platform to check the multiferroicity of LuFeO$_3$ and related materials in the bulk form. Second, the magnetic characterizations on bulk Lu$_{1-x}$Sc$_x$FeO$_3$ demonstrate a magnetic anomaly (probable antiferromagnetic ordering) above room temperature, $\sim425-445$ K, followed by magnetic transitions in low temperatures ($\sim167-172$ K). In addition, a magnetoelectric response is observed in the low temperature region. Our study provides useful information on the multiferroic physics of hexagonal $R$FeO$_3$ and related systems.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The wind nebula around magnetar Swift J1834.9-0846

We report on the analysis of two deep XMM-Newton observations of the magnetar Swift J1834.9-0846 and its surrounding extended emission taken in March 2014 and October 2014, 2.5 and 3.1 years after the source went into outburst. The magnetar is only weakly detected in the first observation with an absorption corrected flux $F_{\rm 0.5-10 keV}\approx4\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$, and a $3σ$ upper limit during the second observation of about $3\times10^{-14}$ erg s$^{-1}$ cm$^{-2}$. This flux level is more than 3 orders of magnitude lower than the flux measured at the outburst onset on September 2011. The extended emission, centered at the magnetar position and elongated towards the south-west, is clearly seen in both observations; it is best fit by a highly absorbed power-law (PL), with a hydrogen column density of $N_{\rm H}=8.0\times10^{22}$ cm$^{-2}$ and PL photon index $Γ=2.2\pm0.2$. Its flux is constant between the two observations at $F_{\rm 0.5-10 keV}=1.3\times10^{-12}$ erg s$^{-1}$ cm$^{-2}$. We find no statistically significant changes in the spectral shape or the flux of this extended emission over a period of 9 years from 2005 to 2014. These new results strongly support the extended emission nature as a wind nebula and firmly establish Swift J1834.9-0846 as the first magnetar to show a surrounding wind nebula. Further, our results imply that such nebulae are no longer exclusive to rotation-powered pulsars and narrow the gap between these two sub-populations of isolated neutron stars. The size and spectrum of the nebula are compatible with those of pulsar-wind nebulae but its radiative efficiency $η_{\rm X}=L_{\rm X}/\dot{E}\approx0.1$ is markedly high, possibly pointing to an additional wind component in Swift J1834.9-0846.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Short Gamma-Ray Bursts with Extended Emission Observed with Swift/BAT and Fermi/GBM

Some short GRBs are followed by longer extended emission, lasting anywhere from ~10 to ~100 s. These short GRBs with extended emission (EE) can possess observational characteristics of both short and long GRBs (as represented by GRB 060614), and the traditional classification based on the observed duration places some of them in the long GRB class. While GRBs with EE pose a challenge to the compact binary merger scenario, they may therefore provide an important link between short and long duration events. To identify the population of GRBs with EE regardless of their initial classifications, we performed a systematic search of short GRBs with EE using all available data (up to February 2013) of both Swift/BAT and Fermi/GBM. The search identified 16 BAT and 14 GBM detected GRBs with EE, several of which are common events observed with both detectors. We investigated their spectral and temporal properties for both the spikes and the EE, and examined correlations among these parameters. Here we present the results of the systematic search as well as the properties of the identified events. Finally, their properties are also compared with short GRBs with EE observed with BATSE, identified through our previous search effort. We found several strong correlations among parameters, especially when all of the samples were combined. Based on our results, a possible progenitor scenario of two-component jet is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Five Year Fermi/GBM Magnetar Burst Catalog

Since launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has detected many hundreds of bursts from magnetar sources. While the vast majority of these bursts have been attributed to several known magnetars, there is also a small sample of magnetar-like bursts of unknown origin. Here we present the Fermi/GBM magnetar catalog, giving the results of the temporal and spectral analyses of 440 magnetar bursts with high temporal and spectral resolution. This catalog covers the first five years of GBM magnetar observations, from July 2008 to June 2013. We provide durations, spectral parameters for various models, fluences and peak fluxes for all the bursts, as well as a detailed temporal analysis for SGR J1550-5418 bursts. Finally, we suggest that some of the bursts of unknown origin are associated with the newly discovered magnetar 3XMM J185246.6+0033.7.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolution of the Major Merger Galaxy Pair Fraction at z < 1

We present a study of the largest available sample of near-infrared selected (i.e., stellar mass selected) dynamically close pairs of galaxies at low redshifts ($z<0.3$). We combine this sample with new estimates of the major-merger pair fraction for stellar mass selected galaxies at $z<0.8$, from the Red Sequence Cluster Survey (RCS1). We construct our low-redshift $K-$band selected sample using photometry from the UKIRT Infrared Deep Sky Survey (UKIDSS) and the Two Micron All Sky Survey (2MASS) in the $K-$band ($\sim 2.2~μ$m). Combined with all available spectroscopy, our $K-$band selected sample contains $\sim 250,000$ galaxies and is $> 90\%$ spectroscopically complete. The depth and large volume of this sample allow us to investigate the low-redshift pair fraction and merger rate of galaxies over a wide range in $K-$band luminosity. We find the major-merger pair fraction to be flat at $\sim 2\%$ as a function of $K-$band luminosity for galaxies in the range $10^8 - 10^{12} L_{\odot}$, in contrast to recent results from studies in the local group that find a substantially higher low-mass pair fraction. This low-redshift major-merger pair fraction is $\sim 40-50\%$ higher than previous estimates drawn from $K-$band samples, which were based on 2MASS photometry alone. Combining with the RCS1 sample we find a much flatter evolution ($m = 0.7 \pm 0.1$), in the relation $f_{\rm{pair}} \propto (1+z)^m$, than indicated in many previous studies. These results indicate that a typical $L\sim L^*$ galaxy has undergone $\sim 0.2-0.8$ major mergers since $z=1$ (depending on the assumptions of merger timescale and percentage of pairs that actually merge).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Time resolved spectroscopy of SGR J1550-5418 bursts detected with Fermi/GBM

We report on time-resolved spectroscopy of the 63 brightest bursts of SGR J1550-5418, detected with Fermi/Gamma-ray Burst Monitor during its 2008-2009 intense bursting episode. We performed spectral analysis down to 4 ms time-scales, to characterize the spectral evolution of the bursts. Using a Comptonized model, we find that the peak energy, E_peak, anti-correlates with flux, while the low-energy photon index remains constant at -0.8 up to a flux limit F~10^-5 erg s-1 cm-2. Above this flux value the E_peak-flux correlation changes sign, and the index positively correlates with flux reaching 1 at the highest fluxes. Using a two black-body model, we find that the areas and fluxes of the two emitting regions correlate positively. Further, we study here for the first time, the evolution of the temperatures and areas as a function of flux. We find that the area-kT relation follows lines of constant luminosity at the lowest fluxes, R^2 \propto kT^-4, with a break at higher fluxes ($F>10^-5.5 erg s-1 cm-2). The area of the high-kT component increases with flux while its temperature decreases, which we interpret as due to an adiabatic cooling process. The area of the low-kT component, on the other hand, appears to saturate at the highest fluxes, towards R_max~30 km. Assuming that crust quakes are responsible for SGR bursts and considering R_max as the maximum radius of the emitting photon-pair plasma fireball, we relate this saturation radius to a minimum excitation radius of the magnetosphere, and put a lower limit on the internal magnetic field of SGR J1550-5418, B_int>~4.5x10^15 G.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Current Sheet Statistics in Three-Dimensional Simulations of Coronal Heating

In a recent numerical study [Ng et al., Astrophys. J. {\bf 747}, 109, 2012], with a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on a series of runs in which random photospheric motions are imposed for hundreds to thousands of \al time in order to obtain converged statistical values. The heating rate found in these simulations saturate to a level that is independent of the Lundquist number. This scaling result was also supported by an analysis with the assumption of the Sweet-Parker scaling of the current sheets, as well as how the width, length and number of current sheets scale with Lundquist number. In order to test these assumptions, we have implemented an automated routine to analyze thousands of current sheets in these simulations and return statistical scalings for these quantities. It is found that the Sweet-Parker scaling is justified. However, some discrepancies are also found and require further study.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Experimental observation of ferrielectricity in multiferroic DyMn2O5: revisiting the electric polarization

The electric polarization and its magnetic origins in multiferroic RMn2O5, where R is rare-earth ion, are still issues under debate. In this work, the temperature-dependent electric polarization of DyMn2O5, the most attractive member of this RMn2O5 family, is investigated using the pyroelectric current method upon varying endpoint temperature of the electric cooling, plus the positive-up-negative-down (PUND) technique. It is revealed that DyMn2O5 at low temperature does exhibit the unusual ferrielectricity rather than ferroelectricity, characterized by two interactive and anti-parallel ferroelectric sublattices which show different temperature-dependences. The two ferroelectric sublattices are believed to be generated from the symmetric exchange-striction mechanisms associated with the Mn-Mn spin interactions and Dy-Mn spin interactions, respectively. The path-dependent electric polarization reflects the first-order magnetic transitions in the low temperature regime. The magnetoelectric effect is mainly attributed to the Dy spin order which is sensitive to magnetic field. The present experiments may be helpful for clarifying the puzzling issues on the multiferroicity in DyMn2O5 and probably other RMn2O5 multiferroics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

SN2009ip: Constraints on the progenitor mass-loss rate

Some supernovae (SNe) show evidence for mass-loss events taking place prior to their explosions. Measuring their pre-outburst mass-loss rates provide essential information regarding the mechanisms that are responsible for these events. Here we present XMM-Newton and Swift X-ray observations taken after the latest, and presumably the final, outburst of SN 2009ip. We use these observations as well as new near infra-red and visible light spectra, and published radio and visible light observations to put six independent order-of-magnitude constrains on the mass-loss rate of the SN progenitor prior to the explosion. Our methods utilize: the X-ray luminosity, the bound-free absorption, the H alpha luminosity, the SN rise-time, free-free absorption, and the bolometric luminosity of the outburst detected prior to the explosion. Assuming spherical mass-loss with a wind density profile, we estimate that the effective mass-loss rate from the progenitor was between 10^-3 to 10^-2 solar masses per year, over a few years prior to the explosion, with a velocity of ~1000 km/s. This mass-loss rate corresponds to a total circum stellar matter mass of ~0.04 solar masses, within 6x10^15 cm of the SN. We note that the mass-loss rate estimate based on the H alpha luminosity is higher by an order of magnitude. This can be explained if the narrow line H alpha component is generated at radii larger than the shock radius, or if the CSM has an aspherical geometry. We discuss simple geometries which are consistent with our results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Energy Distribution of Nanoflares in Three-Dimensional Simulations of Coronal Heating

In a recent computational campaign [Ng et al., Astrophys. J. 747, 109, 2012] to investigate a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on a series of runs in which random photospheric motions are imposed for hundreds to thousands of Alfvén time in order to obtain converged statistical values. Using this collection of numerical data, we have performed additional statistical analysis related to the formation of current sheets and heating events, or nanoflares [Parker, Astrophys. J. 330, 474, 1988]. While there have been many observations of the energy distribution of solar flares, there have not been many results based on large-scale three-dimensional direct simulations due to obvious numerical difficulties. We will present energy distributions and other statistics based on our simulations, calculated using a method employed in [Dmitruk & Gómez, Astrophys. J., 484, L83, 1997]. We will also make comparisons of our results with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-Lundquist Number Scaling in Three-Dimensional Simulations of Parker's Model of Coronal Heating

Parker's model is one of the most discussed mechanisms for coronal heating and has generated much debate. We have recently obtained new scaling results in a two-dimensional (2D) version of this problem suggesting that the heating rate becomes independent of resistivity in a statistical steady state [Ng and Bhattacharjee, Astrophys. J., 675, 899 (2008)]. Our numerical work has now been extended to 3D by means of large-scale numerical simulations. Random photospheric footpoint motion is applied for a time much longer than the correlation time of the motion to obtain converged average coronal heating rates. Simulations are done for different values of the Lundquist number to determine scaling. In the high-Lundquist number limit, the coronal heating rate obtained so far is consistent with a trend that is independent of the Lundquist number, as predicted by previous analysis as well as 2D simulations. In the same limit the average magnetic energy built up by the random footpoint motion tends to have a much weaker dependence on the Lundquist number than that in the 2D simulations, due to the formation of strong current layers and subsequent disruption when the equilibrium becomes unstable. We will present scaling analysis showing that when the dissipation time is comparable or larger than the correlation time of the random footpoint motion, the heating rate tends to become independent of Lundquist number, and that the magnetic energy production is also reduced significantly.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

SGR J1550-5418 bursts detected with the Fermi Gamma-ray Burst Monitor during its most prolific activity

We have performed detailed temporal and time-integrated spectral analysis of 286 bursts from SGR J1550-5418 detected with the Fermi Gamma-ray Burst Monitor (GBM) in January 2009, resulting in the largest uniform sample of temporal and spectral properties of SGR J1550-5418 bursts. We have used the combination of broadband and high time-resolution data provided with GBM to perform statistical studies for the source properties. We determine the durations, emission times, duty cycles and rise times for all bursts, and find that they are typical of SGR bursts. We explore various models in our spectral analysis, and conclude that the spectra of SGR J1550-5418 bursts in the 8-200 keV band are equally well described by optically thin thermal bremsstrahlung (OTTB), a power law with an exponential cutoff (Comptonized model), and two black-body functions (BB+BB). In the spectral fits with the Comptonized model we find a mean power-law index of -0.92, close to the OTTB index of -1. We show that there is an anti-correlation between the Comptonized Epeak and the burst fluence and average flux. For the BB+BB fits we find that the fluences and emission areas of the two blackbody functions are correlated. The low-temperature BB has an emission area comparable to the neutron star surface area, independent of the temperature, while the high-temperature blackbody has a much smaller area and shows an anti-correlation between emission area and temperature. We compare the properties of these bursts with bursts observed from other SGR sources during extreme activations, and discuss the implications of our results in the context of magnetar burst models.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Accounting of the Dust-Obscured Star Formation and Accretion Histories Over the Last ~11~Billion Years

(Abridged) We report on an accounting of the star formation and accretion driven energetics of 24um detected sources in GOODS North. For sources having infrared (IR; 8-1000um) luminosities >3x10^12 L_sun when derived by fitting local SEDs to 24um photometry alone, we find these IR luminosity estimates to be a factor of ~4 times larger than those estimated when the SED fitting includes additional 16 and 70um data (and in some cases mid-infrared spectroscopy and 850um data). This discrepancy arises from the fact that high luminosity sources at z>>0 appear to have far- to mid-infrared ratios, as well as aromatic feature equivalent widths, typical of lower luminosity galaxies in the local Universe. Using our improved estimates for IR luminosity and AGN contributions, we investigate the evolution of the IR luminosity density versus redshift arising from star formation and AGN processes alone. We find that, within the uncertainties, the total star formation driven IR luminosity density is constant between 1.15 < z < 2.35, although our results suggest a slightly larger value at z>2. AGN appear to account for <18% of the total IR luminosity density integrated between 0< z < 2.35, contributing <25% at each epoch. LIRG appear to dominate the star formation rate (SFR) density along with normal star-forming galaxies (L_IR < 10^11 L_sun) between 0.6 < z < 1.15. Once beyond z >2, the contribution from ultraluminous infrared galaxies ULIRGs becomes comparable with that of LIRGs. Using our improved IR luminosity estimates, we find existing calibrations for UV extinction corrections based on measurements of the UV spectral slope typically overcorrect UV luminosities by a factor of ~2, on average, for our sample of 24um-selected sources; accordingly we have derived a new UV extinction correction more appropriate for our sample.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fermi/GBM observations of the ultra-long GRB 091024: A burst with an optical flash

In this paper we examine gamma-ray and optical data of GRB 091024, a gamma-ray burst (GRB) with an extremely long duration of T90~1020 s, as observed with the Fermi Gamma-Ray Burst Monitor (GBM). We present spectral analysis of all three distinct emission episodes using data from Fermi/GBM. Because of the long nature of this event, many ground-based optical telescopes slewed to its location within a few minutes and thus were able to observe the GRB during its active period. We compare the optical and gamma-ray light curves. Furthermore, we estimate a lower limit on the bulk Lorentz factor from the variability and spectrum of the GBM light curve and compare it with that obtained from the peak time of the forward shock of the optical afterglow. From the spectral analysis we note that, despite its unusually long duration, this burst is similar to other long GRBs, i.e. there is spectral evolution (both the peak energy and the spectral index vary with time) and spectral lags are measured. We find that the optical light curve is highly anti-correlated to the prompt gamma-ray emission, with the optical emission reaching the maximum during an epoch of quiescence in the prompt emission. We interpret this behavior as the reverse shock (optical flash), expected in the internal-external shock model of GRB emission but observed only in a handful of GRBs so far. The lower limit on the initial Lorentz factor deduced from the variability time scale ($Γ_{min}=195_{-110}^+{90}$)is consistent within the error to the one obtained using the peak time of the forward shock ($Γ_0=120$) and is also consistent with Lorentz factors of other long GRBs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Large-Scale High-Lundquist Number Reduced MHD Simulations of the Solar Corona Using GPU Accelerated Machines

We have recently carried out a computational campaign to investigate a model of coronal heating in three-dimensions using reduced magnetohydrodynamics (RMHD). Our code is built on a conventional scheme using the pseudo-spectral method, and is parallelized using MPI. The current investigation requires very long time integrations using high Lundquist numbers, where the formation of very fine current layers challenge the resolutions achievable even on massively parallel machines. We present here results of a port to Nvidia CUDA (Compute Unified Device Architecture) for hardware acceleration using graphics processing units (GPUs). In addition to a brief discussion of our general strategy, we will report code performance on several machines which span a variety of hardware configurations and capabilities. These include a desktop workstation with commodity hardware, a dedicated research workstation equipped with four Nvidia C2050 GPUs, as well as several large-scale GPU accelerated distributed memory machines: Lincoln/NCSA, Dirac/NERSC, and Keeneland/NICS.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rest-frame properties of 32 gamma-ray bursts observed by the Fermi Gamma-Ray Burst Monitor

Aims: In this paper we study the main spectral and temporal properties of gamma-ray bursts (GRBs) observed by Fermi/GBM. We investigate these key properties of GRBs in the rest-frame of the progenitor and test for possible intra-parameter correlations to better understand the intrinsic nature of these events. Methods: Our sample comprises 32 GRBs with measured redshift that were observed by GBM until August 2010. 28 of them belong to the long-duration population and 4 events were classified as short/hard bursts. For all of these events we derive, where possible, the intrinsic peak energy in the $νF_ν$ spectrum (\eprest), the duration in the rest-frame, defined as the time in which 90% of the burst fluence was observed (\tninetyrest) and the isotropic equivalent bolometric energy (\eiso). Results: The distribution of \eprest has mean and median values of 1.1 MeV and 750 keV, respectively. A log-normal fit to the sample of long bursts peaks at ~800 keV. No high-\ep population is found but the distribution is biased against low \ep values. We find the lowest possible \ep that GBM can recover to be ~ 15 keV. The \tninetyrest distribution of long GRBs peaks at ~10 s. The distribution of \eiso has mean and median values of $8.9\times 10^{52}$ erg and $8.2 \times 10^{52}$ erg, respectively. We confirm the tight correlation between \eprest and \eiso (Amati relation) and the one between \eprest and the 1-s peak luminosity ($L_p$) (Yonetoku relation). Additionally, we observe a parameter reconstruction effect, i.e. the low-energy power law index $α$ gets softer when \ep is located at the lower end of the detector energy range. Moreover, we do not find any significant cosmic evolution of neither \eprest nor \tninetyrest.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Challenge of Data Reduction for Multiple Instruments on Stratospheric Observatory For Infrared Astronomy (SOFIA)

SOFIA presents a number of interesting challenges for the development of a data reduction environment which, at its initial phase, will have to incorporate pipelines from seven different instruments. Therefore, the SOFIA data reduction software must run code which has been developed in a variety of dissimilar environments which will only increase in future generations of instrumentation. We investigated three distinctly different situations for performing pipelined data reduction in SOFIA: automated data reduction after data archival at the end of a mission, re-pipelining of science data with updated calibrations or optimum parameters, and the interactive user-driven local execution and analysis of data reduction by an investigator. These different modes would traditionally result in very different software implementations of algorithms used by each instrument team, in effect tripling the amount of data reduction software that would need to be maintained by SOFIA. We present here a unique approach for enfolding all the instrument-specific data reduction software in the observatory framework and verifies the needs for all three reduction scenarios as well as the standard visualization tools. The SOFIA data reduction structure would host the different algorithms and techniques that the instrument teams develop in their own environment. Ideally, duplication of software is minimized across the system because instrument teams can draw on software solutions and techniques previously delivered to SOFIA by other instruments. With this approach, we minimize the effort for analyzing and developing new software reduction pipelines for future generation instruments. We also explore the potential benefits of this approach in the portability of the software to an ever-broadening science audience, as well as its ability to ease the use of distributed processing for data reduction pipelines.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A Consistent Model of `Explosive' Financial Bubbles With Mean-Reversing Residuals

We present a self-consistent model for explosive financial bubbles, which combines a mean-reverting volatility process and a stochastic conditional return which reflects nonlinear positive feedbacks and continuous updates of the investors' beliefs and sentiments. The conditional expected returns exhibit faster-than-exponential acceleration decorated by accelerating oscillations, called "log-periodic power law." Tests on residuals show a remarkable low rate (0.2%) of false positives when applied to a GARCH benchmark. When tested on the S&P500 US index from Jan. 3, 1950 to Nov. 21, 2008, the model correctly identifies the bubbles ending in Oct. 1987, in Oct. 1997, in Aug. 1998 and the ITC bubble ending on the first quarter of 2000. Different unit-root tests confirm the high relevance of the model specification. Our model also provides a diagnostic for the duration of bubbles: applied to the period before Oct. 1987 crash, there is clear evidence that the bubble started at least 4 years earlier. We confirm the validity and universality of the volatility-confined LPPL model on seven other major bubbles that have occurred in the World in the last two decades. Using Bayesian inference, we find a very strong statistical preference for our model compared with a standard benchmark, in contradiction with Chang and Feigenbaum (2006) which used a unit-root model for residuals.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Coronal Structure and Abundances in Young Fast Rotators

AB Dor, Speedy Mic and Rst137B are in their early post-T Tauri evolutionary phase (<100Myr), at the age of fastest rotation in the life of late-type stars. They straddle the coronal saturation-supersaturation boundary first defined by young stars in open clusters. High resolution Chandra X-ray spectra have been analysed to study their coronal properties as a function of coronal activity parameters Rossby number, $L_X/L_{bol}$ and a coronal temperature index. Plasma emission measure distributions as a function of temperature show broad peaks at T~10e7K. Differences between stars suggest that as supersaturation is reached the DEM slope below the temperature of peak DEM becomes shallower, while the DEM drop-off above this temperature becomes more pronounced. A larger sample comprising our three targets and 22 active stars studied in the recent literature reveals a general increase of plasma at T>10e7 toward the saturated-supersaturated boundary but a decline beyond this among supersaturated stars. All three of the stars studied in detail here show lower coronal abundances of the low FIP elements Mg, Si and Fe, relative to the high FIP elements S, O and Ne, as compared to the solar mixture. The coronal Fe abundances of the stellar sample are inversely correlated with Lx/Lbol, declining slowly with rising Lx/Lbol, but with a much more sharp decline at Lx/Lbol>3x10e-4. For dwarfs the Fe abundance is also well-correlated with Rossby number. The coronal O/Fe ratios for dwarfs show a clear increase with decreasing Rossby number, apparently reaching saturation at [O/Fe]=0.5 at the coronal supersaturation boundary. Similar increases in O/Fe with increasing coronal temperature and $L_X/L_{bol}$ are seen.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The Evolution of Galaxy Mergers and Morphology at z<1.2 in the Extended Groth Strip

We present the quantitative rest-frame B morphological evolution and galaxy merger fractions at 0.2 < z < 1.2 as observed by the All-wavelength Extended Groth Strip International Survey (AEGIS). We use the Gini coefficent and M_20 to identify major mergers and classify galaxy morphology for a volume-limited sample of 3009 galaxies brighter than 0.4 L_B^*, assuming pure luminosity evolution of 1.3 M_B per unit redshift. We find that the merger fraction remains roughly constant at 10 +/- 2% for 0.2 < z < 1.2. The fraction of E/S0/Sa increases from 21+/- 3% at z ~ 1.1 to 44 +/- 9% at z ~ 0.3, while the fraction of Sb-Ir decreases from 64 +/- 6% at z ~ 1.1 to 47 +/- 9% at z ~ 0.3. The majority of z < 1.2 Spitzer MIPS 24 micron sources with L(IR) > 10^11 L_sun are disk galaxies, and only ~ 15% are classified as major merger candidates. Edge-on and dusty disk galaxies (Sb-Ir) are almost a third of the red sequence at z ~ 1.1, while E/S0/Sa makeup over 90% of the red sequence at z ~ 0.3. Approximately 2% of our full sample are red mergers. We conclude (1) the galaxy merger rate does not evolve strongly between 0.2 < z < 1.2; (2) the decrease in the volume-averaged star-formation rate density since z ~ 1 is a result of declining star-formation in disk galaxies rather than a disappearing population of major mergers; (3) the build-up of the red sequence at z < 1 can be explained by a doubling in the number of spheroidal galaxies since z ~ 1.2.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Bounds on the renormalized couplings in an SU(2)_L \otimes SU(2)_R symmetric Yukawa model

The vacuum stability lower bound on the mass of the Higgs boson is numerically investigated in an $SU(2)_L \otimes SU(2)_R$ symmetric Yukawa model, which describes two heavy degenerate fermion doublets in the limit of vanishing gauge couplings. Good agreement with perturbation theory is found, although the couplings are strong. The upper bound on the fermion mass and renormalized Yukawa coupling is also determined in the part of bare parameter space where reflection positivity has been proven.

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

Mass Spectrum and Bounds on the Couplings in Yukawa Models With Mirror-Fermions

The $\rm SU(2)_L\otimes SU(2)_R$ symmetric Yukawa model with mirror-fermions in the limit where the mirror-fermion is decoupled is studied both analytically and numerically. The bare scalar self-coupling $λ$ is fixed at zero and infinity. The phase structure is explored and the relevant phase transition is found to be consistent with a second order one. The fermionic mass spectrum close to that transition is discussed and a first non-perturbative estimate of the influence of fermions on the upper and lower bounds on the renormalized scalar self-coupling is given. Numerical results are confronted with perturbative predictions.

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

Numerical simulation of heavy fermions in an SU(2)_L x SU(2)_R symmetric Yukawa model

An exploratory numerical study of the influence of heavy fermion doublets on the mass of the Higgs boson is performed in the decoupling limit of a chiral $\rm SU(2)_L \otimes SU(2)_R$ symmetric Yukawa model with mirror fermions. The behaviour of fermion and boson masses is investigated at infinite bare quartic coupling on $4^3 \cdot 8$, $6^3 \cdot 12$ and $8^3 \cdot 16$ lattices. A first estimate of the upper bound on the renormalized quartic coupling as a function of the renormalized Yukawa-coupling is given.

0 citations0 reviewsNo code linkedOpen access