Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic processes.”

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

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

At least 703 records · Page 39

Anomalous decoherence and absence of thermalization in a photonic many-body system

The intention of this work is twofold: first, to present a most simple system capable of simulating the intrinsic bosonic Josephson effect with photons and, second, to study various outcomes deriving from inherent or external decoherence. A qubit induces an effective coupling between two externally pumped cavity modes. Without cavity losses and in the dispersive regime, intrinsic Josephson oscillations of photons between the two modes occurs. In this case, contrary to regular Markovian decoherence, the qubit purity shows a Gaussian decay and recurrence of its coherence. Due to intrinsic nonlinearities, both the Josephson oscillations as well as the qubit properties display a rich collapse-revival structure, where, however, the complexity of the qubit evolution is in some sense stronger. The qubit as a meter of the photon dynamics is considered, and it is shown that qubit dephasing, originating, for example, from nondemolition measurements, results in an exponential destruction of the oscillations which manifests the collectiveness of the Josephson effect. Nonselective qubit measurements, on the other hand, render a Zeno effect seen in a slowing down of the Josephson oscillations. Contrary to dephasing, cavity dissipation results in a Gaussian decay of the scaled Josephson oscillations. Finally, following Ponomarev et al. [Phys. Rev. Lett. 106, 010405 (2011)], we analyze aspects of thermalization. In particular, despite similarities with the generic model studied by Ponomarev et al., our system does not seem to thermalize.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Least squares parameter estimation methods for material decomposition with energy discriminating detectors

Purpose: Energy resolving detectors provide more than one spectral measurement in one image acquisition. The purpose of this study is to investigate, with simulation, the ability to decompose four materials using energy discriminating detectors and least squares minimization techniques. Methods: Three least squares parameter estimation decomposition techniques were investigated for four-material breast imaging tasks in the image domain. The first technique treats the voxel as if it consisted of fractions of all the materials. The second method assumes that a voxel primarily contains one material and divides the decomposition process into segmentation and quantification tasks. The third is similar to the second method but a calibration was used. The simulated computed tomography (CT) system consisted of an 80 kVp spectrum and a CdZnTe (CZT) detector that could resolve the x-ray spectrum into five energy bins. A postmortem breast specimen was imaged with flat panel CT to provide a model for the digital phantoms. Hydroxyapatite (HA) (50, 150, 250, 350, 450, and 550 mg/ml) and iodine (4, 12, 20, 28, 36, and 44 mg/ml) contrast elements were embedded into the glandular region of the phantoms. Calibration phantoms consisted of a 30/70 glandular-to-adipose tissue ratio with embedded HA (100, 200, 300, 400, and 500 mg/ml) and iodine (5, 15, 25, 35, and 45 mg/ml). The x-ray transport process was simulated where the Beer-Lambert law, Poisson process, and CZT absorption efficiency were applied. Qualitative and quantitative evaluations of the decomposition techniques were performed and compared. The effect of breast size was also investigated. Results: The first technique decomposed iodine adequately but failed for other materials. The second method separated the materials but was unable to quantify the materials. With the addition of a calibration, the third technique provided good separation and quantification of hydroxyapatite, iodine, glandular, and adipose tissues. Quantification with this technique was accurate with errors of 9.83% and 6.61% for HA and iodine, respectively. Calibration at one point (one breast size) showed increased errors as the mismatch in breast diameters between calibration and measurement increased. A four-point calibration successfully decomposed breast diameter spanning the entire range from 8 to 20 cm. For a 14 cm breast, errors were reduced from 5.44% to 1.75% and from 6.17% to 3.27% with the multipoint calibration for HA and iodine, respectively. Conclusions: The results of the simulation study showed that a CT system based on CZT detectors in conjunction with least squares minimization technique can be used to decompose four materials. The calibrated least squares parameter estimation decomposition technique performed the best, separating and accurately quantifying the concentrations of hydroxyapatite and iodine.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

ELECTRON ACCELERATION IN PULSAR-WIND TERMINATION SHOCKS: AN APPLICATION TO THE CRAB NEBULA GAMMA-RAY FLARES

The γ -ray flares from the Crab Nebula observed by AGILE and Fermi -LAT reaching GeV energies and lasting several days challenge the standard models for particle acceleration in pulsar-wind nebulae because the radiating electrons have energies exceeding the classical radiation-reaction limit for synchrotron. Previous modeling has suggested that the synchrotron limit can be exceeded if the electrons experience electrostatic acceleration, but the resulting spectra do not agree very well with the data. As a result, there are still some important unanswered questions about the detailed particle acceleration and emission processes occurring during the flares. We revisit the problem using a new analytical approach based on an electron transport equation that includes terms describing electrostatic acceleration, stochastic wave-particle acceleration, shock acceleration, synchrotron losses, and particle escape. An exact solution is obtained for the electron distribution, which is used to compute the associated γ -ray synchrotron spectrum. We find that in our model the γ -ray flares are mainly powered by electrostatic acceleration, but the contributions from stochastic and shock acceleration play an important role in producing the observed spectral shapes. Our model can reproduce the spectra of all the Fermi -LAT and AGILE flares from the Crab Nebula, using magnetic field strengths in agreement with the multi-wavelength observational constraints. We also compute the spectrum and duration of the synchrotron afterglow created by the accelerated electrons, after they escape into the region on the downstream side of the pulsar-wind termination shock. The afterglow is expected to fade over a maximum period of about three weeks after the γ -ray flare.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear intrinsic variables and state reconstruction in multiscale simulations

Finding informative low-dimensional descriptions of high-dimensional simulation data (like the ones arising in molecular dynamics or kinetic Monte Carlo simulations of physical and chemical processes) is crucial to understanding physical phenomena, and can also dramatically assist in accelerating the simulations themselves. In this paper, we discuss and illustrate the use of nonlinear intrinsic variables (NIV) in the mining of high-dimensional multiscale simulation data. In particular, we focus on the way NIV allows us to functionally merge different simulation ensembles, and different partial observations of these ensembles, as well as to infer variables not explicitly measured. The approach relies on certain simple features of the underlying process variability to filter out measurement noise and systematically recover a unique reference coordinate frame. We illustrate the approach through two distinct sets of atomistic simulations: a stochastic simulation of an enzyme reaction network exhibiting both fast and slow time scales, and a molecular dynamics simulation of alanine dipeptide in explicit water.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Continuous-time random walk as a guide to fractional Schroedinger equation

We argue that the continuous-time random walk approach may be a useful guide to extend the Schroedinger equation in order to incorporate nonlocal effects, avoiding the inconsistencies raised by Jeng et al. [J. Math. Phys. 51, 062102 (2010)]. As an application, we work out a free particle in a half space, obtaining the time dependent solution by considering an arbitrary initial condition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THE MASS DISTRIBUTION OF STELLAR-MASS BLACK HOLES

We perform a Bayesian analysis of the mass distribution of stellar-mass black holes using the observed masses of 15 low-mass X-ray binary systems undergoing Roche lobe overflow and 5 high-mass, wind-fed X-ray binary systems. Using Markov Chain Monte Carlo calculations, we model the mass distribution both parametrically-as a power law, exponential, Gaussian, combination of two Gaussians, or log-normal distribution-and non-parametrically-as histograms with varying numbers of bins. We provide confidence bounds on the shape of the mass distribution in the context of each model and compare the models with each other by calculating their relative Bayesian evidence as supported by the measurements, taking into account the number of degrees of freedom of each model. The mass distribution of the low-mass systems is best fit by a power law, while the distribution of the combined sample is best fit by the exponential model. This difference indicates that the low-mass subsample is not consistent with being drawn from the distribution of the combined population. We examine the existence of a 'gap' between the most massive neutron stars and the least massive black holes by considering the value, M{sub 1%}, of the 1% quantile from each black hole mass distribution as the lower bound of black hole masses. Our analysis generates posterior distributions for M{sub 1%}; the best model (the power law) fitted to the low-mass systems has a distribution of lower bounds with M{sub 1%}>4.3 M{sub sun} with 90% confidence, while the best model (the exponential) fitted to all 20 systems has M{sub 1%}>4.5 M{sub sun} with 90% confidence. We conclude that our sample of black hole masses provides strong evidence of a gap between the maximum neutron star mass and the lower bound on black hole masses. Our results on the low-mass sample are in qualitative agreement with those of Ozel et al., although our broad model selection analysis more reliably reveals the best-fit quantitative description of the underlying mass distribution. The results on the combined sample of low- and high-mass systems are in qualitative agreement with Fryer and Kalogera, although the presence of a mass gap remains theoretically unexplained.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-Markovian environments and entanglement preservation

Using the Shabani-Lidar post-Markovian master equation, we derive non-Markovian generalizations of important quantum decohering operations on single qubits. When environmental memory effects are being taken into account, both single-qubit coherence and two-qubit entanglement may be preserved over a longer period of time, in contrast to the corresponding situations where these are totally neglected. We argue that recognizing the fact that every environment is inherently non-Markovian could be the key to the resolution of the issue of entanglement sudden death.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Non-Markovian finite-temperature two-time correlation functions of system operators of a pure-dephasing model

We evaluate the non-Markovian finite-temperature two-time correlation functions (CF's) of system operators of a pure-dephasing spin-boson model in two different ways, one by the direct exact operator technique and the other by the recently derived evolution equations, valid to second order in the system-environment interaction Hamiltonian. This pure-dephasing spin-boson model that is exactly solvable has been extensively studied as a simple decoherence model. However, its exact non-Markovian finite-temperature two-time system operator CF's, to our knowledge, have not been presented in the literature. This may be mainly due to the fact, illustrated in this article, that in contrast to the Markovian case, the time evolution of the reduced density matrix of the system (or the reduced quantum master equation) alone is not sufficient to calculate the two-time system operator CF's of non-Markovian open systems. The two-time CF's obtained using the recently derived evolution equations in the weak system-environment coupling case for this non-Markovian pure-dephasing model happen to be the same as those obtained from the exact evaluation. However, these results significantly differ from the non-Markovian two-time CF's obtained by wrongly directly applying the quantum regression theorem (QRT), a useful procedure to calculate the two-time CF's for weak-coupling Markovian open systems. This demonstrates clearly that the recently derived evolution equations generalize correctly the QRT to non-Markovian finite-temperature cases. It is believed that these evolution equations will have applications in many different branches of physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nakajima-Zwanzig versus time-convolutionless master equation for the non-Markovian dynamics of a two-level system

We consider the exact reduced dynamics of a two-level system coupled to a bosonic reservoir, further obtaining the exact time-convolutionless and Nakajima-Zwanzig non-Markovian equations of motion. The system considered includes the damped and undamped Jaynes-Cummings model. The result is obtained by exploiting an expression of quantum maps in terms of matrices and a simple relation between the time evolution map and the time-convolutionless generator as well as the Nakajima-Zwanzig memory kernel. This nonperturbative treatment shows that each operator contribution in Lindblad form appearing in the exact time-convolutionless master equation is multiplied by a different time-dependent function. Similarly, in the Nakajima-Zwanzig master equation each such contribution is convoluted with a different memory kernel. It appears that, depending on the state of the environment, the operator structures of the two sets of equations of motion can exhibit important differences.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlocal fluctuations and control of dimer entanglement dynamics

We report on the dissipative dynamics of an entangled, bipartite interacting system. We show how to induce and control the so-called early-stage disentanglement (and 'delayed' entanglement generation) dynamics by means of a driving laser field. We demonstrate that some of the features currently associated with pure non-Markovian effects in such entanglement behavior can actually take place in Markovian environments if background-noise quantum electrodynamics fluctuations are considered. We illustrate this for the case of a dimer interacting molecular system for which emission rates, interaction strength, and radiative corrections have been previously measured. We also show that even in the absence of collective decay mechanisms and qubit-qubit interactions, the entanglement still exhibits collapse-revival behavior. Our results indicate that zero-point energy fluctuations should be taken into account when formulating precise entanglement dynamics statements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Separating gravitational wave signals from instrument artifacts

Central to the gravitational wave detection problem is the challenge of separating features in the data produced by astrophysical sources from features produced by the detector. Matched filtering provides an optimal solution for Gaussian noise, but in practice, transient noise excursions or ''glitches'' complicate the analysis. Detector diagnostics and coincidence tests can be used to veto many glitches which may otherwise be misinterpreted as gravitational wave signals. The glitches that remain can lead to long tails in the matched filter search statistics and drive up the detection threshold. Here we describe a Bayesian approach that incorporates a more realistic model for the instrument noise allowing for fluctuating noise levels that vary independently across frequency bands, and deterministic glitch fitting using wavelets as glitch templates, the number of which is determined by a transdimensional Markov chain Monte Carlo algorithm. We demonstrate the method's effectiveness on simulated data containing low amplitude gravitational wave signals from inspiraling binary black-hole systems, and simulated nonstationary and non-Gaussian noise comprised of a Gaussian component with the standard LIGO/Virgo spectrum, and injected glitches of various amplitude, prevalence, and variety. Glitch fitting allows us to detect significantly weaker signals than standard techniques.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic Theory for Distribution Functions of Wave-Particle Interactions in Plasmas

The evolution of a charged particle distribution function under the influence of coherent electromagnetic waves in a plasma is determined from kinetic theory. For coherent waves, the dynamical phase space of particles is an inhomogeneous mix of chaotic and regular orbits. The persistence of long time correlations between the particle motion and the phase of the waves invalidates any simplifying Markovian or statistical assumptions--the basis for usual quasilinear theories. The generalized formalism in this Letter leads to a hierarchy of evolution equations for the reduced distribution function. The evolution operators, in contrast to the quasilinear theories, are time dependent and nonsingular and include the rich phase space dynamics of particles interacting with coherent waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A stochastic model and Monte Carlo algorithm for fluctuation-induced H{sub 2} formation on the surface of interstellar dust grains

A stochastic algorithm for simulation of fluctuation-induced kinetics of H{sub 2} formation on grain surfaces is suggested as a generalization of the technique developed in our recent studies [1] where this method was developed to describe the annihilation of spatially separate electrons and holes in a disordered semiconductor. The stochastic model is based on the spatially inhomogeneous, nonlinear integro-differential Smoluchowski equations with random source term. In this paper we derive the general system of Smoluchowski type equations for the formation of H{sub 2} from two hydrogen atoms on the surface of interstellar dust grains with physisorption and chemisorption sites. We focus in this study on the spatial distribution, and numerically investigate the segregation in the case of a source with a continuous generation in time and randomly distributed in space. The stochastic particle method presented is based on a probabilistic interpretation of the underlying process as a stochastic Markov process of interacting particle system in discrete but randomly progressed time instances. The segregation is analyzed through the correlation analysis of the vector random field of concentrations which appears to be isotropic in space and stationary in time.

79 ASTRONOMY AND ASTROPHYSICS↗

Processes of energy exchange in systems of nonidentical particles with inhomogeneous sources of heat

Processes of energy exchange in dissipative systems of nonidentical interacting particles (having different sizes, charges, etc.) with an inhomogeneous distribution of sources of heat and/or any other sources of stochastic kinetic energy are considered. A theoretical model is proposed for analyzing the energy balance in such systems. Analytical relations describing the redistribution of the “kinetic temperature” between interacting particles of the system are obtained within this model. These relations are tested by the numerical simulation of the problem for Yukawa systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Policy Iteration Algorithm for Average Continuous Control of Piecewise Deterministic Markov Processes

The main goal of this paper is to apply the so-called policy iteration algorithm (PIA) for the long run average continuous control problem of piecewise deterministic Markov processes (PDMP's) taking values in a general Borel space and with compact action space depending on the state variable. In order to do that we first derive some important properties for a pseudo-Poisson equation associated to the problem. In the sequence it is shown that the convergence of the PIA to a solution satisfying the optimality equation holds under some classical hypotheses and that this optimal solution yields to an optimal control strategy for the average control problem for the continuous-time PDMP in a feedback form.

97 MATHEMATICS AND COMPUTING↗

Generation and decay of entanglement in exciton/polariton trapped systems

The entanglement evolution of interacting trapped solid-state bosons, i.e. excitons or polaritons, is theoretically studied. From an initially separable state the quasi-particle interactions generate boson entanglement which is degraded as dissipation takes place. Based on exact solutions of a Markovian Master equation at T = 0 K, we discuss how dissipation affects the quantum entanglement dynamics of coupled solid-state boson systems. The single trap Wigner function is also calculated and its relationship to the trapped boson entanglement is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗