Tensor lattice field theory for renormalization and quantum computing
Not Available
Meurice, Yannick (ORCID:0000000209959694)↗
SEARCH · Engineering Papers
Results for “Computing continuum”
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.
Nonreciprocal Dissipation Engineering via Strong Coupling with a Continuum of Modes
Optical nonreciprocity plays a key role in almost every optical system, directing light flow and protecting optical components from backscattered light. Controllable forms of on-chip nonreciprocity are needed for the robust operation of increasingly sophisticated photonic integrated circuits (PICs) in the context of classical and quantum computation, networking, communications, and sensing. However, it has been challenging to achieve wideband, low-loss optical nonreciprocity on-chip. In this paper, we demonstrate strong coupling and Rabi-like energy exchange between photonic bands, possessing a continuum of modes, to unlock nonreciprocity and frequency translation over wide optical bandwidths in silicon. Using a traveling-wave phonon field to drive indirect interband photonic transitions, we demonstrate band hybridization that enables an intriguing form of nonreciprocal dissipation engineering. Using the converted mode to create a nonreciprocal dissipation channel, we demonstrate a frequency-neutral, low-loss (less than 1 dB) isolator with high nonreciprocal contrast (more than 14 dB) and broad operating bandwidth (more than 59 GHz). Additionally, through the implementation of complete interband conversion, we demonstrate a high extinction (more than 55 dB) optical frequency translation operation with near-unity efficiency.
Zhou, Yishu↗
A predictive discrete-continuum multiscale model of plasticity with quantified uncertainty
Multiscale models of materials, consisting of upscaling discrete simulations to continuum models, are unique in their capability to simulate complex materials behavior. The fundamental limitation in multiscale models is the presence of uncertainty in the computational predictions delivered by them. In this work, a sequential multiscale model has been developed, incorporating discrete dislocation dynamics (DDD) simulations and a strain gradient plasticity (SGP) model to predict the size effect in plastic deformations of metallic micro-pillars. The DDD simulations include uniaxial compression of micro-pillars with different sizes and over a wide range of initial dislocation densities and spatial distributions of dislocations. An SGP model is employed at the continuum level that accounts for the size-dependency of flow stress and hardening rate. Sequences of uncertainty analyses have been performed to assess the predictive capability of the multiscale model. The variance-based global sensitivity analysis determines the effect of parameter uncertainty on the SGP model prediction. The multiscale model is then constructed by calibrating the continuum model using the data furnished by the DDD simulations. A Bayesian calibration method is implemented to quantify the uncertainty due to microstructural randomness in discrete dislocation simulations (density and spatial distribution of dislocations) on the macroscopic continuum model prediction (size effect in plastic deformation). Here, the outcomes of this study indicate that the discrete-continuum multiscale model can accurately simulate the plastic deformation of micro-pillars, despite the significant uncertainty in the DDD results. Additionally, depending on the macroscopic features represented by the DDD simulations, the SGP model can reliably predict the size effect in plasticity responses of the micropillars with below 10% of error.
36 MATERIALS SCIENCE↗
Tardigrade-examples V0.1.0
Tardigrade-examples is a repository of computational workflows that exercise the Tardigrade software package. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling. These capabilities include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This repository investigates continuum upscaling of various direct numerical simulations (DNS) conducted in Abaqus finite element (FE), Ratel FE, and GEOS material point method (MPM) software. Verification of the upscaling workflow is first investigated by considering DNS of trivial stress states for homogeneous materials, results of which indicate that classical continuum behavior is recovered as expected. DNS of heterogeneous materials are then considered.
Allard, Thomas↗
Status of Lattice QCD Calculations of One and Two-Nucleon Matrix Elements
Lattice QCD calculations of basic hadronic quantities are now often performed with systematic control over the physical pion mass, continuum and infinite volume limits. This progress signals a new era in which lattice QCD will be used to compute basic properties of hadrons, nucleons and light nuclei directly from the Standard Model (SM). In particular, we will be able to use lattice QCD, combined with Effective Field Theory, to make quantitative statements about the interaction of SM matter and potential beyond the SM physics ranging from direct dark matter detection to permanent electric dipole moments. I will briefly review the status and challenges of lattice QCD calculations of select one- and two-nucleon matrix elements. These calculations are significantly more challenging than the basic quantities computed with all systematics controlled, but substantial progress is being made.
Walker-Loud, Andre↗
Numerical modeling and experimental validation of low velocity impact of woven GFRP/CFRP composites
Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. Further, the comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.
36 MATERIALS SCIENCE↗
A subdivision-stabilized B-spline mixed material point method
Subjected to external loadings, polymeric materials, e.g., biological tissues, hydrogels, and elastomers, may undergo extreme, nearly incompressible, (self-)contact deformations. For numerical modeling employing mesh-based techniques such as the finite element method (FEM), these deformations pose significant challenges due to large distortions in the deformed geometry, accuracy issues stemming from volumetric locking effects, and increased computational cost from complex contact searches. As an alternative to mesh-based methods, the material point method (MPM), a continuum-based particle technique, is gaining attention for its ability to handle extreme distortions and capture no-slip contact without added cost. For nearly incompressible material behaviors, while mixed formulations can address locking effects by treating displacements and pressure as independent fields, they can suffer from numerical instabilities close to the incompressibility limit due to the violation of the inf-sup condition, leading to inaccurate nodal pressure solutions. Here we propose an efficient and stable mixed B-spline material point method with highest achievable regularity for quasi-compressible polymeric materials. Using the two-scale relation of B-splines, we introduce a subdivision-stabilization for the two-field mixed MPM and obtain numerically stable, oscillation-free nodal solutions with equal-order interpolations with optimal regularity. Building on the Eulerian-Lagrangian nature of MPM, a previously-converged solution framework is adopted to mitigate issues related to cell-crossing and numerical fracture artifact present in standard MPM. We assess the stability and accuracy of the developed mixed MPM at large deformations for soft materials through the benchmark Cook’s membrane problem. Additionally, we test the robustness of the proposed MPM by modeling several examples, including the compression and indentation of a circular block into a quasi-compressible substrate and the twisting deformation of a rectangular block. The findings demonstrate the MPM’s capabilities for modeling practical soft material applications.
36 MATERIALS SCIENCE↗
Bound states and electromagnetic radiation of relativistically rotating cylindrical wells
Here, we compute the effect of rigid rotation on the non-relativistic bound states. The energy levels of the bound states increase with the angular velocity of rotation until at certain value of the angular velocity they are completely pushed out into the continuum which corresponds to dissociation of the bound states. When the angular velocity exceeds the critical value at which the ground state disappears into the continuum, no bound state is possible. This effect should have important consequences for the phenomenology of the quark-gluon plasma. One of the ways to study it experimentally is to observe the electromagnetic radiation emitted by a rotating bound state. We compute the corresponding intensity of electromagnetic radiation and show that it strongly depends on the angular velocity of rotation.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Simulating ice-shelf extent using damage mechanics
Abstract Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner.
54 ENVIRONMENTAL SCIENCES↗
Generalized Hall conductivities in local commuting projector models: Generalized symmetries and protected surface modes
Hall conductivities are important characterizations of phases of matter. It is known that nonzero Hall conductivities are difficult to realize in local commuting projector lattice models due to no-go theorems in (2+1)D. In this work we construct local commuting projector models in (2+1)D and (3+1)D with nonzero generalized Hall conductivities for ordinary and higher-form continuous symmetries on tensor product Hilbert space of finite local dimension. The model is given by a standard ℤ 𝑁 toric code, but the symmetries do not admit expression in terms of on-site charge operators. The symmetry do not have local charges or currents on the lattice in the absence of boundaries, but there is still a notion of Hall conductivities that coincide with the continuum field theories. We construct protected gapless boundaries of the lattice models using modified Villain formalism. The generalized Hall conductivities are computed by surface currents as well as bulk flux insertion and many-body Chern number.
Anomalies↗
Developing Platinum-Group-Metal-Free Catalysts for Oxygen Reduction Reaction in Acid: Beyond the Single Metal Site
This project is to develop M (x) -N-C catalysts with dense multiple metal center (MMC) sites to meet the DOE 2025 activity target of 0.044 mA/cm 2 at 0.9 V IR-free ., as well as other goals such as durability. We have made important contributions to both catalyst development and fundamental understandings of the active sites in M-N-C catalysts in this project. We successfully made M (x) -N-C catalysts with some multiple metal center (MMC) sites by combining ionothermal carbonization with chemical vapor deposition (CVD). These catalysts, however, are not as active as the most active single-atom Fe-N-C catalysts made by the similar CVD process, owing likely to the low site density and the presence of inorganic Fe species such as iron carbides and nanoparticles. The most significant accomplishments we achieved in this project are: (1) we unraveled the formation pathway of Fe-N 4 sites during pyrolysis step-by-step and identified the trans-metalation mechanism, in collaboration with Deborah Myers and her colleagues at Argonne National Laboratory (ANL); (2) inspired by this finding, we pioneered the CVD synthesis of M-N-C catalysts (M = Mn, Fe, and Co), in which the Fe-N-C catalyst by CVD demonstrated an ORR activity of 33 mA/cm 2 at 0.9 V in H 2 -O 2 proton exchange membrane fuel cells (PEMFCs), very close to the ultimate goal of 35 mA/cm 2 of our project. This catalyst is the first Fe-N-C catalyst that contains only D1 sites without the D2 sites; whereas D1 and D2 sites have been always identified by Mossbauer in previous Fe-N-C catalysts. This finding helps to understand what the D1 and D2 sites are and their roles in catalyzing the ORR. (3) by improving the mass transport of the carbon matrix prior to the CVD process, the revised Fe-N-C catalyst made by CVD delivered a maximum power density of 0.53 W/cm 2 in H 2 -air PEMFCs. The improvement strategy was partly inspired by the computational modeling work by Adam Weber from LBNL, the Co-PI of this project, by developing, coding, and exercising a continuum level model of transport phenomena within a PGM-free catalyst layer. The model demonstrated that local resistances combined with limited site density of the PGM-free catalyst can result in limiting currents and poor polarization performance. The model also gave design guidance for impact of ECSA and overall catalyst-layer thickness. However, both Fe-N-C and Co-N-C catalysts developed by CVD showed poor durability in PEMFCS, in comparison with the traditional M-N-C catalysts synthesized via regular pyrolysis process. Consequently, we did not achieve the proposed durability targets. Despite so, we believe the FeNC-CVD catalysts with the poor durability and D1 sites only provides an excellent platform to understand the degradation mechanism of Fe-N-C in PEMFCs, the most important challenge in the development of M-N-C catalysts.
08 HYDROGEN↗
Reconciling Work Functions and Adsorption Enthalpies for Implicit Solvent Models: A Pt (111)/Water Interface Case Study
Implicit solvent models are a computationally efficient method of representing solid/liquid interfaces prevalent in electrocatalysis, energy storage, and materials science. However, electronic structure changes induced at the metallic surface by the dielectric continuum are not fully understood. To address this, we perform DFT calculations for the Pt(111)/water interface, in order to compare Poisson–Boltzmann continuum solvation methods with ab initio molecular dynamics (AIMD) simulations of explicit solvent. We show that the implicit solvent cavity can be parametrized in terms of the electric dipole moment change at the equilibrated explicit Pt/water interface to obtain the potential of zero charge (PZC). We also compare the accuracy of aqueous enthalpies of adsorption of phenol on Pt(111) using geometry and charge density based dielectric cavitation methods. The ability to parametrize the cavity according to individual atoms, as afforded in the geometry based approach, is key to obtaining accurate enthalpy changes of adsorption under aqueous conditions. Additionally, we show that the electronic structure changes induced by explicit solvent and our proposed implicit solvent parametrization scheme yield comparable density difference profiles and d-band projected density of states. We therefore demonstrate the capability of implicit solvent approaches to capture both the energetics of adsorption processes and the main electronic effects of aqueous solvent on the metallic surface. Therefore, this work provides a scheme for computationally efficient simulations of interfacial processes for applications in areas such as heterogeneous catalysis and electrochemistry.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Quadrupole-mediated dielectric response and the charge-asymmetric solvation of ions in water
Treating water as a linearly responding dielectric continuum on molecular length scales allows very simple estimates of the solvation structure and thermodynamics for charged and polar solutes. While this approach can successfully account for basic length and energy scales of ion solvation, computer simulations indicate not only its quantitative inaccuracies but also its inability to capture some basic and important aspects of microscopic polarization response. Here, we consider one such shortcoming, a failure to distinguish the solvation thermodynamics of cations from that of otherwise-identical anions, and we pursue a simple, physically inspired modification of the dielectric continuum model to address it. The adaptation is motivated by analyzing the orientational response of an isolated water molecule whose dipole is rigidly constrained. Its free energy suggests a Hamiltonian for dipole fluctuations that accounts implicitly for the influence of higher-order multipole moments while respecting constraints of molecular geometry. Finally, we propose a field theory with the suggested form, whose nonlinear response breaks the charge symmetry of ion solvation. An approximate variational solution of this theory, with a single adjustable parameter, yields solvation free energies that agree closely with simulation results over a considerable range of solute size and charge.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Correlation effects in magic-angle twisted bilayer graphene: An auxiliary-field quantum Monte Carlo study
Magic-angle twisted bilayer graphene (MATBG) presents a fascinating platform for investigating the effects of electron interactions in topological flat bands. The Bistritzer-MacDonald (BM) model provides a simplified quantitative description of the flat bands. Introducing long-range Coulomb interactions leads to an interacting BM (IBM) Hamiltonian, a momentum-space continuum description which offers a very natural starting point for many-body studies of MATBG. Accurate and reliable many-body computations in the IBM model are challenging, however, and have been limited mostly to special fillings or smaller lattice sizes. We employ a state-of-the-art auxiliary-field quantum Monte Carlo (AFQMC) method to study the IBM model, which constrains the sign problem to enable accurate treatment of large system sizes. We determine ground-state properties and quantify errors compared to mean-field theory calculations. Our calculations identify correlated metal states and their competition with the insulating Kramers intervalley-coherent state at both half-filling and charge neutrality. Additionally, we investigate one- and three-quarter fillings, and examine the effect of many-body corrections beyond single Slater determinant solutions. We discuss the effect that details of the IBM Hamiltonian have on the results, including different forms of double-counting corrections, and the need to establish and precisely specify many-body Hamiltonians to allow more direct and quantitative comparisons with experiments in MATBG. Published by the American Physical Society 2025
Xiao, Zhi-Yu (ORCID:0000000219531579)↗
Radiative Electronic Bound States in the Continuum from Defects in Semiconductors
Continuum-buried defect states in semiconductors are generally expected to be optically inactive because of their strong coupling to continuum bands. Here, we show that such defects can instead host radiative electronic bound states in the continuum (BICs) using the silicon G center as a prototypical example. Hybrid functional first-principles calculations with a Hubbard U correction reveal that a localized defect state, initially buried below the valence band maximum (VBM) in the ground state, undergoes exchange-driven energy-level reordering under optical excitation and shifts above the VBM. This exchange-induced transition suppresses nonradiative decay and enables a robust radiative emission. By computing temperature-dependent nonradiative lifetimes and comparing them to experimental photoluminescence (PL) lifetimes, we quantitatively reproduce the observed temperature dependence of the emission. These results uncover a stabilization mechanism for continuum-embedded defect states and establish electronic BICs as a general paradigm for designing defect-based optical systems, including quantum emitters and qubits.
Hong, SeongYun↗
Impact of artificial topological changes on flow and transport through fractured media due to mesh resolution
Abstract We performed a set of numerical simulations to characterize the interplay of fracture network topology, upscaling, and mesh refinement on flow and transport properties in fractured porous media. We generated a set of generic three-dimensional discrete fracture networks at various densities, where the radii of the fractures were sampled from a truncated power-law distribution, and whose parameters were loosely based on field site characterizations. We also considered five network densities, which were defined using a dimensionless version of density based on percolation theory. Once the networks were generated, we upscaled them into a single continuum model using the upscaled discrete fracture matrix model presented by Sweeney et al. (2019). We considered steady, isothermal pressure-driven flow through each domain and then simulated conservative, decaying, and adsorbing tracers using a pulse injection into the domain. For each simulation, we calculated the effective permeability and solute breakthrough curves as quantities of interest to compare between network realizations. We found that selecting a mesh resolution such that the global topology of the upscaled mesh matches the fracture network is essential. If the upscaled mesh has a connected pathway of fracture (higher permeability) cells but the fracture network does not, then the estimates for effective permeability and solute breakthrough will be incorrect. False connections cannot be eliminated entirely, but they can be managed by choosing appropriate mesh resolution and refinement for a given network. Adopting octree meshing to obtain sufficient levels of refinement leads to fewer computational cells (up to a 90% reduction in overall cell count) when compared to using a uniform resolution grid and can result in a more accurate continuum representation of the true fracture network.
58 GEOSCIENCES↗
Data driven learning of robust nonlocal models
Nonlocal models use integral operators that embed length-scales in their definition. However, the integrands in these operators are difficult to define from the data that are typically available for a given physical system, such as laboratory mechanical property tests. In contrast, molecular dynamics (MD) does not require these integrands, but it suffers from computational limitations in the length and time scales it can address. To combine the strengths of both methods and to obtain a coarse-grained, homogenized continuum model that efficiently and accurately captures materials' behavior, we propose a learning framework to extract, from MD data, an optimal nonlocal model as a surrogate for MD displacements. Our framework guarantees that the resulting model is mathematically well-posed, physically consistent, and that it generalizes well to settings that are different from the ones used during training. The efficacy of this approach is demonstrated with several numerical tests for single layer graphene both in the case of perfect crystal and in the presence of thermal noise.
36 MATERIALS SCIENCE↗
On the Viability of Quantum Annealers to Solve Fluid Flows
This paper explores the suitability of upcoming novel computing technologies, particularly adiabatic annealing based quantum computers, to solve fluid dynamics problems that form a critical component of several science and engineering applications. For our experiments, we start with a well-studied one-dimensional simple flow problem, and provide a framework to convert such problems in continuum to a form amenable for deployment on such quantum annealers. Since the DWave annealer returns multiple states sampling the energy landscape of the problem, we explore multiple solution selection strategies to approximate the solution of the problem. We analyze the continuum solutions obtained both qualitatively and quantitatively as well as their sensitivities to the particular solution selection scheme.
42 ENGINEERING↗