Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Random walks”

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 361 records · Page 20

MONACO: accurate biological network alignment through optimal neighborhood matching between focal nodes

Motivation: Alignment of protein–protein interaction networks can be used for the unsupervised prediction of functional modules, such as protein complexes and signaling pathways, that are conserved across different species. To date, various algorithms have been proposed for biological network alignment, many of which attempt to incorporate topological similarity between the networks into the alignment process with the goal of constructing accurate and biologically meaningful alignments. Especially, random walk models have been shown to be effective for quantifying the global topological relatedness between nodes that belong to different networks by diffusing node-level similarity along the interaction edges. However, these schemes are not ideal for capturing the local topological similarity between nodes. Results: Here, we propose MONACO, a novel and versatile network alignment algorithm that finds highly accurate pairwise and multiple network alignments through the iterative optimal matching of ‘local’ neighborhoods around focal nodes. Extensive performance assessment based on real networks as well as synthetic networks, for which the ground truth is known, demonstrates that MONACO clearly and consistently outperforms all other state-of-the-art network alignment algorithms that we have tested, in terms of accuracy, coherence and topological quality of the aligned network regions. Furthermore, despite the sharply enhanced alignment accuracy, MONACO remains computationally efficient and it scales well with increasing size and number of networks.

97 MATHEMATICS AND COMPUTING↗

CHRR: coordinate hit-and-run with rounding for uniform sampling of constraint-based models

In constraint-based metabolic modelling, physical and biochemical constraints define a polyhedral convex set of feasible flux vectors. Uniform sampling of this set provides an unbiased characterization of the metabolic capabilities of a biochemical network. However, reliable uniform sampling of genome-scale biochemical networks is challenging due to their high dimensionality and inherent anisotropy. Here, we present an implementation of a new sampling algorithm, coordinate hit-and-run with rounding (CHRR). This algorithm is based on the provably efficient hit-and-run random walk and crucially uses a preprocessing step to round the anisotropic flux set. CHRR provably converges to a uniform stationary sampling distribution. We apply it to metabolic networks of increasing dimensionality. We show that it converges several times faster than a popular artificial centering hit-and-run algorithm, enabling reliable and tractable sampling of genome-scale biochemical networks.

59 BASIC BIOLOGICAL SCIENCES↗

Effects of linear trends on estimation of noise in GNSS position time-series

A thorough understanding of time-dependent noise in Global Navigation Satellite System (GNSS) position time-series is necessary for computing uncertainties in any signals found in the data. However, estimation of time-correlated noise is a challenging task and is complicated by the difficulty in separating noise from signal, the features of greatest interest in the time-series. In this study, we investigate how linear trends affect the estimation of noise in daily GNSS position time-series. We use synthetic time-series to study the relationship between linear trends and estimates of time-correlated noise for the six most commonly cited noise models. We find that the effects of added linear trends, or conversely de-trending, vary depending on the noise model. The commonly adopted model of random walk (RW), flicker noise (FN) and white noise (WN) is the most severely affected by de-trending, with estimates of low-amplitude RW most severely biased. FN plus WN is least affected by adding or removing trends. Non-integer power-law noise estimates are also less affected by de-trending, but are very sensitive to the addition of trend when the spectral index is less than one. We derive an analytical relationship between linear trends and the estimated RW variance for the special case of pure RW noise. Finally, overall, we find that to ascertain the correct noise model for GNSS position time-series and to estimate the correct noise parameters, it is important to have independent constraints on the actual trends in the data.

58 GEOSCIENCES↗

Blazar variability power spectra from radio up to TeV photon energies: Mrk 421 and PKS 2155−304

ABSTRACT We present the results of the power spectral density (PSD) analysis for the blazars Mrk 421 and PKS 2155−304, using good-quality, densely sampled light curves at multiple frequencies, covering 17 decades of the electromagnetic spectrum, and variability time-scales from weeks up to a decade. The data were collected from publicly available archives of observatories at radio from Owens Valley Radio Observatory, optical and infrared (B, V, R, I, J, H, and Kbands), X-rays from the Swift and the Rossi X-ray Timing Explorer, high and very high energy (VHE) γ-rays from the Fermi and Very Energetic Radiation Imaging Telescope Array System as well as the High Energy Stereoscopic System. Our results are: (1) the power-law form of the variability power spectra at radio, infrared, and optical frequencies have slopes ∼1.8, indicative of random-walk-type noise processes; (2) the power-law form of the variability power spectra at higher frequencies, from X-rays to VHE γ-rays, however, have slopes ∼1.2, suggesting a flicker noise-type process; and (3) there is significantly more variability power at X-rays, high and VHE γ-rays on time-scales ≲ 100 d, as compared to lower energies. Our results do not easily fit into a simple model, in which a single compact emission zone is dominating the radiative output of the blazars across all the time-scales probed in our analysis. Instead, we argue that the frequency-dependent shape of the variability power spectra points out a more complex picture, with highly inhomogeneous outflow producing non-thermal emission over an extended, stratified volume.

Goyal, Arti↗

Optical variability of quasars with 20-yr photometric light curves

ABSTRACT We study the optical gri photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during ∼1998−2020 with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average ∼200 nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10–15 yr baselines and ≲ 100 epochs. We find that the average damping time-scale τDRW continues to rise with increased baseline, reaching a median value of ∼750 d (g band) in the rest frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained τDRW (less than 20 per cent of the baseline), we confirm a weak wavelength dependence of τDRW∝λ0.51 ± 0.20. We further quantify optical variability of these quasars over days to decades time-scales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping time-scales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on time-scales less than ∼ a month for these luminous quasars, and this second break point correlates with the longer DRW damping time-scale.

79 ASTRONOMY AND ASTROPHYSICS↗

OzDES Reverberation Mapping Programme: Mg ii lags and R − L relation

ABSTRACT The correlation between the broad line region radius and continuum luminosity (R–L relation) of active galactic nuclei (AGNs) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At z ∼ 1–2, where AGN activity peaks, the R–L relation is constrained by the reverberation mapping (RM) lags of the Mg ii line. We present 25 Mg ii lags from the Australian Dark Energy Survey RM project based on 6 yr of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the rescaled, resampled versions of the observed light curves of local, well-measured AGN. Our sample significantly increases the number of Mg ii lags and extends the R–L relation to higher redshifts and luminosities. The relative iron line strength $\mathcal {R}_{\rm Fe}$ has little impact on the R–L relation. The best-fitting Mg iiR–L relation has a slope α = 0.39 ± 0.08 with an intrinsic scatter $\sigma _{\rm rl} = 0.15^{+0.03}_{-0.02}$ . The slope is consistent with previous measurements and shallower than the H β R–L relation. The intrinsic scatter of the new R–L relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H β R–L relation. Our new R–L relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon.

79 ASTRONOMY AND ASTROPHYSICS↗

Searching for quasar candidates with periodic variations from the Zwicky Transient Facility: results and implications

ABSTRACT We conduct a systematic search for quasars with periodic variations from the archival photometric data of the Zwicky Transient Facility by cross-matching with the quasar catalogues of the Sloan Digital Sky Survey and Véron-Cetty and Véron. We first select out 184 primitive periodic candidates using the generalized Lomb–Scargle periodogram and autocorrelation function and then estimate their statistical significance of periodicity based on two red-noise models, i.e. damped random walk (DRW) and single power-law (SPL) models. As such, we finally identify 106 (DRW) and 86 (SPL) candidates with the most significant periodic variations out of 143 700 quasars. We further compare DRW and SPL models using Bayes factors, which indicate a relative preference of the SPL model for our primitive sample. We thus adopt the candidates identified with SPL as the final sample and summarize its basic properties. We extend the light curves of the selected candidates by supplying other archival survey data to verify their periodicity. However, only three candidates (with 6–8 cycles of periods) meet the selection criteria. This result clearly implies that, instead of being strictly periodic, the variability must be quasi-periodic or caused by stochastic red-noise. This exerts a challenge to the existing search approaches and calls for developing new effective methods.

Chen, Yong-Jie (ORCID:0000000342807673)↗

RoboPol: first season rotations of optical polarization plane in blazars

Here, we present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma-rays is investigated using the data set obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations during non-rotating periods have significantly larger amplitude and faster variations of polarization angle than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability ≤1.5 × 10 -2 ) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely (~5 × 10 -5 ) that none of our rotations is physically connected with an increase in gamma-ray activity.

79 ASTRONOMY AND ASTROPHYSICS↗

The ASAS-SN catalogue of variable stars – IV. Periodic variables in the APOGEE survey

ABSTRACT We explore the synergy between photometric and spectroscopic surveys by searching for periodic variable stars among the targets observed by the Apache Point Observatory Galactic Evolution Experiment (APOGEE) using photometry from the All-Sky Automated Survey for Supernovae (ASAS-SN). We identified 1924 periodic variables among more than $258\, 000$ APOGEE targets; 465 are new discoveries. We homogeneously classified 430 eclipsing and ellipsoidal binaries, 139 classical pulsators (Cepheids, RR Lyrae, and δ Scuti), 719 long-period variables (pulsating red giants), and 636 rotational variables. The search was performed using both visual inspection and machine learning techniques. The light curves were also modelled with the damped random walk stochastic process. We find that the median [Fe/H] of variable objects is lower by 0.3 dex than that of the overall APOGEE sample. Eclipsing binaries and ellipsoidal variables are shifted to a lower median [Fe/H] by 0.2 dex. Eclipsing binaries and rotational variables exhibit significantly broader spectral lines than the rest of the sample. We make ASAS-SN light curves for all the APOGEE stars publicly available and provide parameters for the variable objects.

Pawlak, Michał↗

Measurement of the longitudinal flow-plane decorrelation using multiplane cumulants in $\sqrt{s_{NN}}$ = 200 GeV Au+Au, Ru+Ru, and Zr+Zr collisions

Measurements of the variation of anisotropic flow-plane angles (𝛹 𝑛 ) with rapidity, commonly known as the flow-plane decorrelation, provide important insights into the initial conditions of the matter produced in heavy-ion collisions. Here, in this paper, using data collected by the STAR experiment, we report the first measurement of the four-plane correlator observable 𝑇 𝑛 ⁡{𝑏⁢𝑎;𝑑⁢𝑐}=⟨⟨sin⁡[𝑛⁢(𝛹$^𝑏_𝑛$ − 𝛹$^𝑎_𝑛$)]⁢sin⁡[𝑛⁢(𝛹$^𝑑_𝑛$ − 𝛹$^𝑐_𝑛$)]⟩⟩, where superscripts 𝑎, 𝑏, 𝑐, and 𝑑 denote sequential pseudorapidity (𝜂) regions with 𝑎 corresponding to the most backward region, 𝑏 and 𝑐 close to mid-rapidity with 𝜂 𝑏 < 0 and 𝜂 𝑐 > 0, and 𝑑 being the most forward. The measurement is performed for the elliptic and triangular flow (i.e.~𝑛 = 2 and 3) in Au+Au and isobar (Ru+Ru, Zr+Zr) collisions at $\sqrt{s_{NN}}$=200~GeV. The goal of calculating the correlation of the flow-plane angle variations from backward to mid-central, and from mid-central to forward regions, is to probe the systematic variation of flow angle over a wide 𝜂 range. In mid-central collisions (10−30% centrality), we find 𝑇 2 ⁡{𝑏⁢𝑎;𝑑⁢𝑐} =−0.004 ±0.001⁢(s⁢t⁢a⁢t) ±0.002⁢(s⁢y⁢s⁢t) independent of the collision system. Such a small value of 𝑇 2 favors a random-walk'' variation of the flow-plane angles, where the rapidity correlation length is smaller than the entire region under study. These measurements provide new information on the decorrelation patterns in the system and offer a quantitative estimate of possible systematic variations in anisotropic flow angles such astwist’’ between forward and backward regions. This opens new opportunities for understanding the three-dimensional structure and the time evolution of the quark-gluon plasma created in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Trapping interlayer excitons in van der Waals heterostructures by potential arrays

Transition metal dichalcogenide heterostructures can host interlayer excitons (IXs), which consist of electrons and holes spatially separated in different layers. IXs possess permanent dipoles and have proven to offer a wealth of novel physics. Here we develop a discrete, random-walk model which includes annihilation and repulsion interactions among IXs. Using this model, we simulate the trapping of IXs in traps of different depths, densities, and shapes. Our results show that dipole-dipole interactions play an important role in regulating IX trapping. The effects of dipole interactions can be mitigated with small, deep traps which are realizable with atomic defects and moire potentials.

2-dimensional systems↗

Distinct universality classes of diffusive transport from full counting statistics

The hydrodynamic transport of local conserved densities furnishes an effective coarse-grained description of the dynamics of a many-body quantum system. However, the full quantum dynamics contains much more structure beyond the simplified hydrodynamic description. Here we show that systems with the same hydrodynamics can nevertheless belong to distinct dynamical universality classes, as revealed by new classes of experimental observables accessible in synthetic quantum systems, which can, for instance, measure simultaneous site-resolved snapshots of all of the particles in a system. Specifically, we study the full counting statistics of spin transport, whose first moment is related to linear-response transport, but the higher moments go beyond. We present an analytic theory of the full counting statistics of spin transport in various integrable and nonintegrable anisotropic one-dimensional spin models, including the XXZ spin chain. We find that spin transport, while diffusive on average, is governed by a distinct non-Gaussian dynamical universality class in the models considered. We consider a setup in which the left and right half of the chain are initially created at different magnetization densities, and consider the probability distribution of the magnetization transferred between the two half-chains. We derive a closed-form expression for the probability distribution of the magnetization transfer, in terms of random walks on the half-line. We show that this distribution strongly violates the large-deviation form expected for diffusive chaotic systems, and explain the physical origin of this violation. Here, we discuss the crossovers that occur as the initial state is brought closer to global equilibrium. Our predictions can directly be tested in experiments using quantum gas microscopes or superconducting qubit arrays.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

C 60 -induced Devil's Staircase transformation on a Pb/Si(111) wetting layer

Density functional theory is used to study structural energetics of Pb vacancy cluster formation on C 60 /Pb/Si(111) to explain the unusually fast and error-free transformations between the “Devil's Staircase” (DS) phases on the Pb/Si(111) wetting layer at low temperature (~110K). The formation energies of vacancy clusters are calculated in C 60 /Pb/Si(111) as Pb atoms are progressively ejected from the initial dense Pb wetting layer. Vacancy clusters larger than five Pb atoms are found to be stable with seven being the most stable, while vacancy clusters smaller than five are highly unstable, which agrees well with the observed ejection rate of ~5 Pb atoms per C 60 . Furthermore, the high energy cost (~0.8 eV) for the small vacancy clusters to form indicates convincingly that the unusually fast transformation observed experimentally between the DS phases, upon C 60 adsorption at low temperature, cannot be the result of single-atom random walk diffusion but of correlated multi-atom processes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Hydrodynamics of operator spreading and quasiparticle diffusion in interacting integrable systems

In this work, we address the hydrodynamics of operator spreading in interacting integrable lattice models. In these models, operators spread through the ballistic propagation of quasiparticles, with an operator front whose velocity is locally set by the fastest quasiparticle velocity. In interacting integrable systems, this velocity depends on the density of the other quasiparticles, so equilibrium density fluctuations cause the front to follow a biased random walk, and therefore to broaden diffusively. Ballistic front propagation and diffusive front broadening are also generically present in nonintegrable systems in one dimension; thus, although the mechanisms for operator spreading are distinct in the two cases, these coarse-grained measures of the operator front do not distinguish between the two cases. We present an expression for the front-broadening rate; we explicitly derive this for a particular integrable model (the “Floquet-Fredrickson-Andersen” model), and argue on kinetic grounds that it should apply generally. In conclusion, our results elucidate the microscopic mechanism for diffusive corrections to ballistic transport in interacting integrable models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Primary fission fragment mass yields across the chart of nuclides

In this work, we have calculated a complete set of primary fission fragment mass yields, Y(A), for heavy nuclei across the chart of nuclides, including those of particular relevance to the rapid neutron capture process (r process) of nucleosynthesis. We assume that the nuclear shape dynamics are strongly damped, which allows for a description of the fission process via Brownian shape motion across nuclear potential-energy surfaces. The macroscopic energy of the potential was obtained with the Finite-Range Liquid-Drop Model (FRLDM), while the microscopic terms were extracted from the single-particle level spectra in the fissioning system by the Strutinsky procedure for the shell energies and the BCS treatment for the pairing energies. For each nucleus considered, the fission fragment mass yield, Y(A), is obtained from 50 000 to 500 000 random walks on the appropriate potential-energy surface. The full mass and charge yield, Y(Z,A), is then calculated by invoking the Wahl systematics. With this method, we have calculated a comprehensive set of fission-fragment yields from over 3800 nuclides bounded by 80 ≤ Z ≤ 130 and A ≤ 330; these yields are provided as an ASCII formatted database in the Supplemental Material. We compare our yields to known data and discuss general trends that emerge in low-energy fission yields across the chart of nuclides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Formation of transfermium elements in reactions with Pb 208

Within the Langevin framework, we investigate the dynamics of the fusion process for production of transfermium elements in reactions of Ca 48 , Ti 50 , Cr 54 , and Fe 58 with Pb 208 . After the reacting nuclei have made contact, the early dynamical stage is dominated by the dissipation of the initial radial kinetic energy, while the subsequent shape evolution is diffusive. The probability for surmounting the inner barrier and forming a compound system is obtained by simulating the evolution as a Metropolis random walk in a five-dimensional potential-energy landscape. Good agreement with the available data is obtained, especially for the maximal formation probability. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Calculated fission-fragment yield systematics in the region 74 <=Z <= 94 and 90 <=N <= 150

Background: In the seminal experiment by Schmidt et al . [Nucl. Phys. A 665 , 221 (2000)] in which fission-fragment charge distributions were obtained for 70 nuclides, asymmetric distributions were seen above nucleon number A ≈ 226 and symmetric ones below. Because asymmetric fission had often loosely been explained as a preference for the nucleus to always exploit the extra binding of fragments near ¹³²Sn it was assumed that all systems below A ≈ 226 would fission symmetrically because available isotopes do not have a proton-to-neutron Z/N ratio that allows division into fragments near ¹³²Sn. But the finding by Andreyev et al . [Phys. Rev. Lett. 105 , 252502 (2010)] did not conform to this expectation because the compound system ¹⁸⁰Hg was shown to fission asymmetrically. It was suggested that this was a new type of asymmetric fission, because no strong shell effects occur for any possible fragment division. Purpose: We calculate a reference database for fission-fragment mass yields for a large region of the nuclear chart comprising 987 nuclides. A particular aim is to establish whether ¹⁸⁰Hg is part of a contiguous region of asymmetric fission, and if so, its extent, or if not, in contrast to the actinides, there are scattered smaller groups of nuclei that fission asymmetrically in this area of the nuclear chart. Methods: We use the by now well benchmarked Brownian shape-motion method and perform random walks on the previously calculated five-dimensional potential-energy surfaces. The calculated shell corrections are damped out with energy according to a prescription developed earlier. Results: We have obtained a theoretical reference database of fission-fragment mass yields for 987 nuclides. These results show an extended region of asymmetric fission with approximate extension 74 ≤ Z ≤ 85 and 100 ≤ N ≤ 120. The calculated yields are highly variable. We show 20 representative plots of these variable features and summarize the main aspects of our results in terms of “nuclear-chart” plots showing calculated degrees of asymmetry versus N and Z. Conclusions: Experimental data in this region are rare: only ten or so yield distributions have been measured, some with very limited statistics. We agree with several measurements with higher statistics. Regions where there might be differences between our calculated results and measurements lie near the calculated transition line between symmetric and asymmetric fission. To draw more definite conclusions about the accuracy of the present implementation of the Brownian shape-motion approach in this region experimental data, with reliable statistics, for a fair number of suitably located additional nuclides are clearly needed. Because the nuclear potential-energy structure is so different in this region compared to the actinide region, additional experimental data together with fission theory studies that incorporate additional, dynamical aspects should provide much new insight.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. Here, the method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015 E h of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10 –5 E h without an importance sampling or Jastrow factor.

74 ATOMIC AND MOLECULAR PHYSICS↗