skip to main content

This content will become publicly available on August 21, 2024

Title: Development of a new quantum trajectory molecular dynamics framework
An extension to the wave packet description of quantum plasmas is presented, where the wave packet can be elongated in arbitrary directions. A generalized Ewald summation is constructed for the wave packet models accounting for long-range Coulomb interactions and fermionic effects are approximated by purpose-built Pauli potentials, self-consistent with the wave packets used. We demonstrate its numerical implementation with good parallel support and close to linear scaling in particle number, used for comparisons with the more common wave packet employing isotropic states. Ground state and thermal properties are compared between the models with differences occurring primarily in the electronic subsystem. Especially, the electrical conductivity of dense hydrogen is investigated where a 15% increase in DC conductivity can be seen in our wave packet model compared with other models. This article is part of the theme issue ‘Dynamic and transient processes in warm dense matter’.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    A suite of general circulation models is used to investigate the surface magnetic perturbations due to the ionospheric currents driven by an eastward‐propagating ultrafast Kelvin wave (UFKW) packet with periods between 2 and 4 days and zonal wave number. The simulated daytime UFKW‐driven meridional magnetic perturbations dBn (∼±5 nT) (or zonal currents) between about 5° and 20° magnetic latitude in each hemisphere are opposite in sign to those equatorward of±5° and produced by the equatorial electrojet (EEJ), with the directions on any given day determined by the phase of the UFKW as it propagates eastward with respect to the sunlit ionosphere. Since the nominal daytimeSqzonal current between∼±30° is uniformly eastward flowing, the present results are consistent with the hypothesis that the EEJ is part of a local current vortex with oppositely directed currents near the equator versus those between 5° and 20° at low latitudes. UFKWs are a special wave type wherein meridional winds are relatively small, which leads to our finding that the EEJ dBn constitutes a simple quantitative proxy forE‐region UFKW neutral winds near the 107‐km peak height of the Hall conductivity, including the variable wave period of the UFKW packet. Numerical experiments are also performed to understand the longitude distribution of actual ground magnetometer measurements that are needed to reliably extract the UFKW dBn signal and hence the neutral winds, both of which are closely linked to plasma drifts and electron densities in the equatorialFregion. Using actual magnetometer data it is moreover shown that the UFKW dBn signal is easily measurable. Therefore measurements of EEJ dBn can potentially be used to infer UFKW activity for scientific investigations focusing on coupling between the tropical troposphere and the ionosphere‐thermosphere.

    more » « less
  2. Abstract

    Numerical experiments are performed using a suite of general circulation models that enable the interaction between a Kelvin wave packet and the ionosphere‐thermosphere (IT) to be elucidated. Focus is on an eastward‐propagating ultra‐fast Kelvin wave (UFKW) packet with periods between 2 and 4 days and zonal wavenumbers=−1 during day of year (DOY) 266–281, 2009. Dissipative processes modify the classic UFKW dynamics (equatorially trapped, small meridional wind component) in three ways: (1) molecular diffusion acts to spread the UFKW zonal (u) and meridional (v) wind fields meridionally, pole to pole, asuandv, respectively, decay and grow with increasing height; (2) due to molecular diffusion, the UFKW spectrum at longer periods and with shorter vertical wavelengths preferentially dissipates with height; and (3) interaction with the diurnally varying IT introduces a westward‐propagatings=+2 component to the wind field that significantly modifies its longitude‐UT structure to include a diurnal modulation. The F‐region ionosphere also responds withs=+2, which originates from the influence of diurnally varying E‐region conductivity onE×Bdrifts. Additional spectral peaks invand ionospheric parameters arise due to longitude variations in the magnetic field. Maximum excursions in NmF2 (as compared with those from a simulation without UFKW forcing) achieve values as large as ±50% but more commonly occur in the range of ±20–30%. The combination of positive and negative responses, and their relative magnitudes, depends on the phasing of the UFKW as it moves zonally relative to the Sun‐synchronous diurnal variation of the ionosphere, in addition to its changing amplitude between DOY 266 and 282. Modifications of order 10 ms−1and −7% to zonal‐mean zonal winds and NmF2, respectively, also result from dissipation of the UFKW packet.

    more » « less
  3. We have developed a memory and operation-count efficient 2.5D inversion algorithm of electrical resistivity (ER) data that can handle fine discretization domains imposed by other geophysical (e.g, ground penetrating radar or seismic) data. Due to numerical stability criteria and available computational memory, joint inversion of different types of geophysical data can impose different grid discretization constraints on the model parameters. Our algorithm enables the ER data sensitivities to be directly joined with other geophysical data without the need of interpolating or coarsening the discretization. We have used the adjoint method directly in the discretized Maxwell’s steady state equation to compute the data sensitivity to the conductivity. In doing so, we make no finite-difference approximation on the Jacobian of the data and avoid the need to store large and dense matrices. Rather, we exploit matrix-vector multiplication of sparse matrices and find successful convergence using gradient descent for our inversion routine without having to resort to the Hessian of the objective function. By assuming a 2.5D subsurface, we are able to linearly reduce memory requirements when compared to a 3D gradient descent inversion, and by a power of two when compared to storing a 2D Hessian. Moreover, our method linearly outperforms operation counts when compared with 3D Gauss-Newton conjugate-gradient schemes, which scales cubically in our favor with respect to the thickness of the 3D domain. We physically appraise the domain of the recovered conductivity using a cutoff of the electric current density present in our survey. We evaluate two case studies to assess the validity of our algorithm. First, on a 2.5D synthetic example, and then on field data acquired in a controlled alluvial aquifer, where we were able to match the recovered conductivity to borehole observations. 
    more » « less
  4. A bstract We study damping signatures at the Jiangmen Underground Neutrino Observatory (JUNO), a medium-baseline reactor neutrino oscillation experiment. These damping signatures are motivated by various new physics models, including quantum decoherence, ν 3 decay, neutrino absorption, and wave packet decoherence. The phenomenological effects of these models can be characterized by exponential damping factors at the probability level. We assess how well JUNO can constrain these damping parameters and how to disentangle these different damping signatures at JUNO. Compared to current experimental limits, JUNO can significantly improve the limits on τ 3 / m 3 in the ν 3 decay model, the width of the neutrino wave packet σ x , and the intrinsic relative dispersion of neutrino momentum σ rel . 
    more » « less
  5. The temperature dependence of the thermal conductivity is linked to the nature of the energy transport at a frequency ω , which is quantified by thermal diffusivity d ( ω ). Here we study d ( ω ) for a poorly annealed glass and a highly stable glass prepared using the swap Monte Carlo algorithm. To calculate d ( ω ), we excite wave packets and find that the energy moves diffusively for high frequencies up to a maximum frequency, beyond which the energy stays localized. At intermediate frequencies, we find a linear increase of the square of the width of the wave packet with time, which allows for a robust calculation of d ( ω ), but the wave packet is no longer well described by a Gaussian as for high frequencies. In this intermediate regime, there is a transition from a nearly frequency independent thermal diffusivity at high frequencies to d ( ω ) ∼ ω −4 at low frequencies. For low frequencies the sound waves are responsible for energy transport and the energy moves ballistically. The low frequency behavior can be predicted using sound attenuation coefficients. 
    more » « less