Tensor lattice field theory for renormalization and quantum computing
Not Available
Meurice, Yannick (ORCID:0000000209959694)↗
SEARCH · Engineering Papers
Results for “continuum modeling”
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.
Theory of subcycle time-resolved photoemission: Application to terahertz photodressing in graphene
Here, motivated by recent experimental progress we revisit the theory of pump–probe time- and angle-resolved photoemission spectroscopy (trARPES), which is one of the most powerful techniques to trace transient pump-driven modifications of the electronic properties. The pump-induced dynamics can be described in different gauges for the light–matter interaction. Standard minimal coupling leads to the velocity gauge, defined by linear coupling to the vector potential. In the context of tight-binding (TB) models, the Peierls substitution is the commonly employed scheme for single-band models. Multi-orbital extensions – including the coupling of the dipole moments to the electric field – have been introduced and tested recently. In this work, we derive the theory of time-resolved photoemission within both gauges from the perspective of nonequilibrium Green’s functions. This approach naturally incorporates the photoelectron continuum, which allows for a direct calculation of the observable photocurrent. Following this route we introduce gauge-invariant expressions for the time-resolved photoemission signal. The theory is applied to graphene pumped with short terahertz pulses, which we treat within a first-principles TB model. We investigate the gauge invariance and discuss typical effects observed in subcycle time-resolved photoemission. Our formalism is an ideal starting point for realistic trARPES simulations including scattering effects.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
A multiphase flow model of water droplets dielectrophoretic-induced air dehumidification phenomena
Air humidity in indoor spaces plays a critical role in human comfort and health. Dehumidification systems are used for building humidity controls, but they can take significant energy consumption, especially in geographic locations with high outdoor humidity and warm climates. Consequently, there is a growing demand for innovative dehumidification processes that consume minimal energy. Dielectrophoretic air dehumidification represents one such promising approach. However, it has not garnered significant attention due to the absence of engineering models and simulation tools capable of evaluating its performance and limitations at large-scale airflows. A new numerical multiphase CFD model, which is also experimentally validated, is developed in a customized Reacting Foam solver based on OpenFOAM® version 9. The newly developed model seeks to decrease substantial energy consumption and lower costs by leveraging the dielectrophoretic phenomenon to regulate moisture levels in the air. The solver integrates a hybrid Eulerian-Lagrangian framework to track the droplet's trajectory and growth rate while solving the continuum equations for the moist air. An electrospray produces electrically charged droplets, which grow during their in-flight trajectories as water vapor condenses onto their surfaces. The role of electrostatic forces in promoting vapor condensation within a high-gradient electrical field is investigated, and the dielectrophoretic vapor nucleation process on charged water droplets is discussed. The CFD model was validated against results from the literature and from proof-of-concept experiments conducted by the authors, which showed a 2 % air dehumidification with a single electrospray and airflow rate of 5 cubic feet per minute. The simulation results indicated that augmenting the number of electrically charged spray droplets increased the dehumidification of the air to 25 %. The initial mean droplet diameter, the orientation of the injector and relative humidity significantly influence the assessment of dehumidification. As a result, scaling up this approach to larger airflow volumes is identified as a potential future research direction.
42 ENGINEERING↗
An Evaluation of Advanced Tools for Distributed Wind Turbine Performance Estimation
We evaluate various classes of distributed wind turbine performance tools across two sites in the United States. The class of tools ranges from the simple mass conservation model to the coupled Reynolds-averaged Navier-Stokes model, all initiated by the WIND Toolkit data set. The resource estimation at the site is evaluated against measured data at the mast location. Taking a sample 100-kW wind turbine and constant losses, we evaluate Openwind, Continuum, and WindNinja tools and document annual energy production (AEP) and time-series statistics associated with the performance estimation of the wind turbine. Using a methodology that is consistent and unbiased across the three options currently available in the industry, we elaborate results at the two sample locations and discuss the probable sources of discrepancy in the AEP estimates. Two main sources of the discrepancy come from the input WIND Toolkit data and the spatial modeling techniques of the tools to capture atmospheric physics. The discussion includes additional values that these tools may bring into the energy assessment process to enhance the owners' confidence over the distributed wind power systems.
17 WIND ENERGY↗
Electromagnetic coupling in tight-binding models for strongly correlated light and matter
In this work, we discuss the construction of low-energy tight-binding Hamiltonians for condensed-matter systems with a strong coupling to the quantum electromagnetic field. Such Hamiltonians can be obtained by projecting the continuum theory on a given set of Wannier orbitals. However, different representations of the continuum theory lead to different low-energy formulations because different representations may entangle light and matter, transforming orbitals into light-matter hybrid states before the projection. In particular, a multicenter Power-Zienau-Woolley transformation yields a dipolar Hamiltonian which incorporates the light-matter coupling via both Peierls phases and a polarization density. We compare this dipolar gauge Hamiltonian and the straightforward Coulomb gauge Hamiltonian for a one-dimensional solid to describe subcycle light-driven electronic motion in the semiclassical limit and a coupling of the solid to a quantized cavity mode which renormalizes the band-structure into electron-polariton bands. Both descriptions yield the same result when many bands are taken into account but the dipolar Hamiltonian is more accurate when the model is restricted to few electronic bands, while the Coulomb Hamiltonian requires fewer electromagnetic modes.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Surrogate-constructed scalable-circuits adaptive variational quantum eigensolver in the Schwinger model
Inspired by recent advancements in simulating periodic systems on quantum computers, we develop an approach to further advance the simulation of these systems, named (SC) 2 -ADAPT-VQE. Our approach extends the scalable-circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable “surrogate constructed” method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC) 2 methodology. Our approach allows for a wider set of classical computations on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute the properties of the Schwinger model—quantum electrodynamics for a single, massive fermion in 1 +1 dimensions—and show that our method can be used to accurately extrapolate to the continuum limit.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Surrogate Constructed Scalable Circuits ADAPT-VQE in the Schwinger model
Inspired by recent advancements of simulating periodic systems on quantum computers, we develop a new approach, (SC)$^2$-ADAPT-VQE, to further advance the simulation of these systems. Our approach extends the scalable circuits ADAPT-VQE framework, which builds an ansatz from a pool of coordinate-invariant operators defined for arbitrarily large, though not arbitrarily small, volumes. Our method uses a classically tractable ``Surrogate Constructed'' method to remove irrelevant operators from the pool, reducing the minimum size for which the scalable circuits are defined. Bringing together the scalable circuits and the surrogate constructed approaches forms the core of the (SC)$^2$ methodology. Our approach allows for a wider set of classical computations, on small volumes, which can be used for a more robust extrapolation protocol. While developed in the context of lattice models, the surrogate construction portion is applicable to a wide variety of problems where information about the relative importance of operators in the pool is available. As an example, we use it to compute properties of the Schwinger model - quantum electrodynamics for a single, massive fermion in $1+1$ dimensions - and show that our method can be used to accurately extrapolate to the continuum limit.
Gustafson, Erik [RIACS, Mtn. View] (ORCID:00000001↗
Quantification of the Impact of Acidified Brine on Fracture-Matrix Transport in a Naturally Fractured Shale Using in Situ Imaging and Modeling
Understanding flow, transport, chemical reactions, and hydromechanical processes in fractured geologic materials is key for optimizing a range of subsurface processes including carbon dioxide and hydrogen storage, unconventional energy resource extraction, and geothermal energy recovery. Flow and transport processes in naturally fractured shale rocks have been challenging to characterize due to experimental complexity and the multiscale nature of quantifying continuum scale descriptions of mass exchange between micrometer-scale fractures and nanometer-scale pores. In this study, we use positron emission tomography (PET) to image the transport of a conservative tracer in a naturally fractured Wolfcamp shale core before and after the core was exposed to low pH brine conditions. Image-based experimental observations are interpreted by fitting an analytical transport model to fracture-containing voxels in the core. Results of this analysis indicate subtle increases in matrix diffusivity and a slightly more uniform fracture velocity distribution following exposure to low pH conditions. These observations are compared with a multicomponent one-dimensional reactive transport model that indicates the capacity for a 10% increase in porosity at the fracture-matrix interface as a result of the low pH brine exposure. This porosity change is the result of the dissolution of carbonate minerals in the shale matrix to low pH conditions. Further, this image-based workflow represents a new approach for quantifying spatially resolved fracture-matrix transport processes and provides a foundation for future work to better understand the role of coupled transport, reaction, and mechanical processes in naturally fractured rocks.
58 GEOSCIENCES↗
Ising on $\mathbb{S}^2$-- The Affine Conjecture
We review the recent construction \cite{brower2024isingmodelmathbbs2} of the 2d Ising model on a triangulated sphere $\mathbb{S}^2$. Surprisingly, this led to a precise map of the lattice couplings to the target geometry in order to reach the conform field theory (CFT) in the continuum limit. For the integrable 2d Ising CFT, the map was found analytically \cite{Brower_2023}. Here we conjecture how this might be generalized. The discrete geometry is implemented by the piecewise flat triangulation introduced by Regge in 1960 for the Einstein Hilbert action \cite{Regge1961GeneralRW}. Then following our Ising example, we posit the existence of a smooth map of lattice couplings in affine parameters consistent with quantum correlators. A sequence of theoretical investigations and numerical simulations are recommended to test this conjecture. They begin with non-integrable CFT's -- the 2d $\phi^4$ theory on $\mathbb{S}^2$; the 3d Ising model on $\mathbb{S}^3$ and $\mathbb{R} \times \mathbb{S}^2$; QED3 on $\mathbb{R} \times \mathbb{S}^{2} $ as an intermediate step to 4d non-Abelian lattice gauge theory on $\mathbb{R} \times \mathbb{S}^3$.
Brower, Richard C. [Boston U.] (ORCID:000000028185↗
A Three-Dimensional, Thermodynamically and Variationally Consistent, Fully Coupled, Electro-Chemo-Thermo-Mechanical Model of Solid-State Batteries
Here we present a theoretical framework for high fidelity modeling of three-dimensional electro-chemo-thermo-mechanical coupled systems such as solid-state batteries. The framework provides a useful means of utilizing molecular level information to predict the overall performance of the system at the continuum level. The central statement of the theoretical framework is a thermodynamically consistent inf-sup problem. The total potential of the coupled system consists of the Helmholtz free energy, the electrical kinetic potential, the chemical potential, the Fourier potential, the chemical reaction potential, and the external power expenditure. All governing equations derive from variational principles. We also show that the resulting initial-boundary-value problem (IBVP) is consistent with existing, widely adopted equations for modeling batteries. Moreover, the IBVP can be naturally discretized using the (variational) finite element method to simulate practical solid-state batteries with arbitrary geometries. Several examples illustrate multiphysics coupling between pairs of the physical field equations, culminating in a comprehensive three-dimensional model of charge-discharge in a solid state battery.
25 ENERGY STORAGE↗
Computational Methods for Modeling Electrospray Microdroplet Chemistry for Improved Quantitative Mass Spectrometry
This project aimed at enhancing the quantitative analysis capabilities of electrospray ionization mass spectrometry (ESI-MS) by developing advanced computational methods. The primary focus was to integrate continuum and molecular dynamics simulations to study the behavior of microdroplets in the ESI process, from formation to evaporation. Through this research, we sought to bridge significant length and time scales to provide a comprehensive understanding of how analyte concentrations evolve from bulk solutions into gas-phase ions. This understanding is crucial for addressing challenges such as ionization efficiency, solvent effects, and ion suppression, which currently limit the accuracy of quantitative ESI-MS.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability
Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.
79 ASTRONOMY AND ASTROPHYSICS↗
Identifying Band Inversions in Topological Materials Using Diffusion Monte Carlo
Topological insulators are characterized by insulating bulk states and robust metallic surface states. Band inversion is a hallmark of topological insulators. At time-reversal invariant points in the Brillouin zone, spin–orbit coupling (SOC) induces a swapping of orbital character at the bulk band edges. Reliably detecting band inversion in solid-state systems with many-body methods would aid in identifying possible candidates for spintronics and quantum computing applications and improve our understanding of the physics behind topologically nontrivial systems. Density functional theory (DFT) methods are a well-established means of investigating these interesting materials due to their favorable balance of computational cost and accuracy but often struggle to accurately model the electron–electron correlations present in the many materials containing heavier elements. In this work, we develop a novel method to detect band inversion within continuum quantum Monte Carlo (QMC) methods that can accurately treat the electron correlation and spin–orbit coupling that are crucial to the physics of topological insulators. Our approach applies a momentum-space-resolved atomic population analysis throughout the first Brillouin zone utilizing the Löwdin method and the one-body reduced density matrix produced with diffusion Monte Carlo (DMC). We integrate this method into QMCPACK, an open source ab initio QMC package, so that these ground-state methods can be used to complement experimental studies and validate prior DFT work on predicting the band structures of correlated topological insulators. Here, we demonstrate this new technique on the topological insulator bismuth telluride, which displays band inversion between its Bi-p and Te-p states at the Γ-point. We show an increase in charge on the bismuth-p orbital and a decrease in charge on the tellurium-p orbital when comparing band structures with and without SOC. Additionally, we use our method to compare the degree of band inversion present in monolayer Bi 2 Te 3 , which has no interlayer van der Waals interactions, to that seen in the bilayer and bulk. The method presented here will enable future many-body studies of band inversion that can shed light on the delicate interplay between correlation and topology in correlated topological materials.
Band structure↗
Implicit full-F simulations of neoclassical ion transport
The development of implicit time integration capabilities for axisymmetric full-F continuum simulations of ion neoclassical transport is reported. The approach involves the implicit treatment of the gyrokinetic Vlasov equation coupled to the nonlinear Fokker–Planck collision model in the long-wavelength limit approximation. To facilitate implicit simulations, advanced preconditioning of individual physics operators is developed, and a global multi-physics preconditioner is constructed by adopting an operator splitting methodology. The algorithm is implemented in the finite-volume code COGENT and is applied to study neoclassical transport properties for both the main ion species and the lithium impurity species in the closed-field-line region of the LTX- β tokamak. The implicit COGENT simulations elucidate the role of non-local transport effects, while demonstrating substantial speedup over the corresponding explicit approach.
Dorf, Mikhail [Lawrence Livermore National Laborat↗
A cautionary tale of attenuation in star-forming regions
ABSTRACT The attenuation of light from star-forming galaxies is correlated with a multitude of physical parameters including star formation rate, metallicity and total dust content. This variation in attenuation is even more evident on kiloparsec scales, which is the relevant size for many current spectroscopic integral field unit surveys. To understand the cause of this variation, we present and analyse Swift/UVOT near-UV (NUV) images and SDSS/MaNGA emission-line maps of 29 nearby (z < 0.084) star-forming galaxies. We resolve kiloparsec-sized star-forming regions within the galaxies and compare their optical nebular attenuation (i.e. the Balmer emission line optical depth, $\tau ^{l}_{B}\equiv \tau _{\textrm {H}\beta }-\tau _{\textrm {H}\alpha }$) and NUV stellar continuum attenuation (via the NUV power-law index, β) to the attenuation law described by Battisti et al. We show the data agree with that model, albeit with significant scatter. We explore the dependence of the scatter of the β–$\tau ^{l}_{B}$ measurements from the star-forming regions on different physical parameters, including distance from the nucleus, star formation rate and total dust content. Finally, we compare the measured $\tau ^{l}_{B}$ and β values for the individual star-forming regions with those of the integrated galaxy light. We find a strong variation in β between the kiloparsec scale and the larger galaxy scale that is not seen in $\tau ^{l}_{B}$. We conclude that the sightline dependence of UV attenuation and the reddening of β due to the light from older stellar populations could contribute to the scatter in the β–$\tau ^{l}_{B}$ relation.
Molina, Mallory↗
Evaluating quenching in cosmological simulations of galaxy formation with spectral covariance in the optical window
ABSTRACT Cosmological hydrodynamical simulations provide valuable insights on galaxy evolution when coupled with observational data. Comparisons with real galaxies are typically performed via scaling relations of the observables. Here, we follow an alternative approach based on the spectral covariance in a model-independent way. We build upon previous work by Sharbaf et al. that studied the covariance of high-quality SDSS (Sloan Digital Sky Survey) continuum-subtracted spectra in a relatively narrow range of velocity dispersion ($\sigma \in [100,150]$ km s$^{-1}$). Here, the same analysis is applied to synthetic data from the eagle and IllustrisTNG100 simulations, to assess the ability of these runs to mimic real galaxies. The real and simulated spectra are consistent regarding spectral covariance, although with subtle differences that can inform the implementation of subgrid physics. Spectral fitting done a posteriori on stacks segregated with respect to latent space reveals that the first principal component (PC1) is predominantly influenced by the stellar age distribution, with an underlying age–metallicity degeneracy. Good agreement is found regarding star formation prescriptions but there is disagreement with active galactic nucleus (AGN) feedback, that also affects the subset of quiescent galaxies. We show a substantial difference in the implementation of the AGN subgrid prescriptions, regarding central black hole seeding, that could lead to the mismatch. Differences are manifest between these two simulations in the star formation histories stacked with respect to latent space. We emphasize that this methodology only relies on the spectral covariance to assess whether simulations provide a true representation of galaxy formation.
Sharbaf, Z. (ORCID:0009000450545946)↗
Generalized susceptibilities and the properties of charm degrees of freedom across the QCD crossover temperature
We study the generalized charm susceptibilities in 2+1 flavor QCD on the lattice at several lattice spacings. We show that, below the chiral crossover, these susceptibilities are well described by the hadron resonance gas (HRG) model if charmed hadrons not listed in tables of the Particle Data Group are included. However, the HRG description abruptly breaks down just above the chiral crossover. To understand this, we use a model for the charm pressure in which it is expressed as the sum of partial pressures from charmed baryons, charmed mesons, and charm quarks. We present continuum estimates of these partial pressures and find that, while the partial pressures of charmed mesons and baryons drop below their respective HRG predictions, the charm quark pressure becomes nonzero above the chiral crossover.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds
Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.
97 MATHEMATICS AND COMPUTING↗