Engineering Papers⌕ Search

Engineering topics

Meurice, Y.

Publications and source records attributed to Meurice, Y..

16 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Truncation effects in the charge representation of the O(2) model

The O(2) model in Euclidean space-time is the zero-gauge-coupling limit of the compact scalar quantum electrodynamics. In this work, we obtain a dual representation of it called the charge representation. We study the quantum phase transition in the charge representation with a truncation to “spin $\textit{S}$," where the quantum numbers have an absolute value less than or equal to $\textit{S}$. The charge representation preserves the gapless-to-gapped phase transition even for the smallest spin truncation $\textit{S}$ = 1. The phase transition for $\textit{S}$ = 1 is an infinite-order Gaussian transition with the same critical exponents $\textit{δ}$ and $\textit{η}$ as the Berezinskii-Kosterlitz-Thouless (BKT) transition, while there are true BKT transitions for $\textit{S}$ ≥ 2 . The essential singularity in the correlation length for $\textit{S}$ = 1 is different from that for $\textit{S}$ ≥ 2 . The exponential convergence of the phase-transition point is studied in both Lagrangian and Hamiltonian formulations. We discuss the effects of replacing the truncated $\hat{U}^±$ = exp ($± i\hat{θ}$) operators by the spin ladder operators $\hat{S}^±$ in the Hamiltonian. The marginal operators vanish at the Gaussian transition point for $\textit{S}$ = 1, which allows us to extract the $\textit{η}$ exponent with high accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantum Joule expansion of one-dimensional systems

Here, we investigate the Joule expansion of nonintegrable quantum systems that contain bosons or spinless fermions in one-dimensional lattices. A barrier initially confines the particles to be in half of the system in a thermal state described by the canonical ensemble and is removed at time t = 0. We investigate the properties of the time-evolved density matrix, the diagonal ensemble density matrix, and the corresponding canonical ensemble density matrix with an effective temperature determined by the total energy conservation using exact diagonalization. The weights for the diagonal ensemble and the canonical ensemble match well for high initial temperatures that correspond to negative effective final temperatures after the expansion. At long times after the barrier is removed, the time-evolved Rényi entropy of subsystems bigger than half can equilibrate to the thermal entropy with exponentially small fluctuations. The time-evolved reduced density matrix at long times can be approximated by a thermal density matrix for small subsystems. Few-body observables, like the momentum distribution function, can be approximated by a thermal expectation of the canonical ensemble with strongly suppressed fluctuations. The negative effective temperatures for finite systems go to non-negative temperatures in the thermodynamic limit for bosons, but there is a true thermodynamic effect for fermions, which is confirmed by finite-temperature density matrix renormalization group calculations. We propose the Joule expansion as a way to dynamically create negative temperature states for fermion systems with repulsive interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

B s → K ℓ ν decay from lattice QCD

We use lattice QCD to calculate the form factors f+(q 2 ) and f 0 (q 2 ) for the semileptonic decay Bs→Kℓν. Our calculation uses six MILC asqtad 2+1 flavor gauge-field ensembles with three lattice spacings. At the smallest and largest lattice spacing the light-quark sea mass is set to 1/10 the strange-quark mass. At the intermediate lattice spacing, we use four values for the light-quark sea mass ranging from 1/5 to 1/20 of the strange-quark mass. We use the asqtad improved staggered action for the light valence quarks, and the clover action with the Fermilab interpolation for the heavy valence bottom quark. We use SU(2) hard-kaon heavy-meson rooted staggered chiral perturbation theory to take the chiral-continuum limit. A functional z expansion is used to extend the form factors to the full kinematic range. We present predictions for the differential decay rate for both Bs→Kμν and Bs→Kτν. We also present results for the forward-backward asymmetry, the lepton polarization asymmetry, ratios of the scalar and vector form factors for the decays Bs→Kℓν and Bs→Dsℓν. Our results, together with future experimental measurements, can be used to determine the magnitude of the Cabibbo-Kobayashi-Maskawa matrix element |$V_{ub}$|.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Linear sigma model for multiflavor gauge theories

We consider a linear sigma model describing $2N_f^2$ bosons ($\sigma$, a 0 , $\eta$' and π) as an approximate effective theory for a $SU(3)$ local gauge theory with $N_f$ Dirac fermions in the fundamental representation. The model has a renormalizable U(N f ) L Ⓧ invariant part, which has an approximate $O(2N_f^2)$ symmetry, and two additional terms, one describing the effects of a $SU(N_f)_V$ invariant mass term and the other the effects of the axial anomaly. We calculate the spectrum for arbitrary $N_f$. Using preliminary and published lattice results from the LatKMI collaboration, we found combinations of the masses that vary slowly with the explicit chiral symmetry breaking and $N_f$. This suggests that the anomaly term plays a leading role in the mass spectrum and that simple formulas such as $M_\sigma^2\simeq (2/N_f-C_\sigma)M_{\eta'}^2$ should apply in the chiral limit. Lattice measurements of $M_{\eta'}^2$ and of approximate constants such as $C_\sigma$ could help locating the boundary of the conformal window. Here, we show that our calculation can be adapted for arbitrary representations of the gauge group and in particular to the minimal model with two sextets, where similar patterns are likely to apply.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Semileptonic $B$-Meson Decays to Light Pseudoscalar Mesons on the HISQ Ensembles

We report the status of an ongoing lattice-QCD calculation of form factors for exclusive semileptonic decays of $B$ mesons with both charged currents ($B\to\pi\ell\nu$, $B_s\to K\ell\nu$) and neutral currents ($B\to\pi\ell^+\ell^-$, $B\to K\ell^+\ell^-$). The results are important for constraining or revealing physics beyond the Standard Model. This work uses MILC's (2+1+1)-flavor ensembles with the HISQ action for the sea and light valence quarks and the clover action in the Fermilab interpretation for the $b$ quark. Simulations are carried out at three lattice spacings down to $0.088$ fm, with both physical and unphysical sea-quark masses. We present preliminary results for correlation-function fits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Probing the conformal Calabrese-Cardy scaling with cold atoms

We demonstrate that current experiments using cold bosonic atoms trapped in one-dimensional optical lattices and designed to measure the second-order Rényi entanglement entropy S 2 can be used to verify detailed predictions of conformal field theory (CFT) and estimate the central charge c. We discuss the adiabatic preparation of the ground state at half filling and small hopping parameter J/U, where we expect a CFT with c = 1. We provide two complementary methods to estimate and subtract the classical entropy generated by the experimental preparation and imaging processes. We compare numerical calculations for the classical O(2) model with a chemical potential on a (1 + 1)-dimensional lattice, and the quantum Bose-Hubbard Hamiltonian implemented in the experiments. S 2 is very similar for the two models and follows closely the Calabrese-Cardy scaling, (c/8) ln (N s ) , for N s sites with open boundary conditions, provided that the large subleading corrections are taken into account.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Gauge-invariant implementation of the Abelian-Higgs model on optical lattices

We present a gauge-invariant effective action for the Abelian Higgs model (scalar electrodynamics) with a chemical potential $\mu$ on a 1+1 dimensional lattice. This formulation provides an expansion in the hopping parameter $\kappa$ which we test with Monte Carlo simulations for a broad range of the inverse gauge coupling $\beta_{pl}=1/g^2$ and small values of the scalar self-coupling $\lambda$. In the opposite limit of infinitely large $\lambda$, the partition function can be written as a traced product of local tensors which allows us to write exact blocking formulas. Gauss's law is automatically satisfied and the introduction of $\mu$ has consequences only if we have an external electric field, $g^2=0$ or an explicit gauge symmetry breaking. The time continuum limit of the blocked transfer matrix can be obtained numerically and, for $g^2=0$ and a spin-1 truncation, the small volume energy spectrum is identical to the low energy spectrum of a two-species Bose-Hubbard model in the limit of large on-site repulsion. We extend this procedure for finite $\beta_{pl}$ and derive a spin-1 approximation of the Hamiltonian. It involves new terms corresponding to transitions among the two species in the Bose-Hubbard model. Finally, we propose an optical lattice implementations involving a ladder structure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Progress towards quantum simulating the classical O ( 2 ) mode

We connect explicitly the classical $O(2)$ model in 1+1 dimensions, a model sharing important features with $U(1)$ lattice gauge theory, to physical models potentially implementable on optical lattices and evolving at physical time. Using the tensor renormalization group formulation, we take the time continuum limit and check that finite dimensional projections used in recent proposals for quantum simulators provide controllable approximations of the original model. We propose two-species Bose-Hubbard models corresponding to these finite dimensional projections at strong coupling and discuss their possible implementations on optical lattices using a $^{87}$Rb and $^{41}$K Bose-Bose mixture.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Controlling sign problems in spin models using tensor renormalization

We consider the sign problem for classical spin models at complex $\beta =1/g_0^2$ on $L\times L$ lattices. We show that the tensor renormalization group method allows reliable calculations for larger Im$\beta$ than the reweighting Monte Carlo method. For the Ising model with complex $\beta$ we compare our results with the exact Onsager-Kaufman solution at finite volume. The Fisher zeros can be determined precisely with the TRG method. We check the convergence of the TRG method for the O(2) model on $L\times L$ lattices when the number of states $D_s$ increases. We show that the finite size scaling of the calculated Fisher zeros agrees very well with the Kosterlitz-Thouless transition assumption and predict the locations for larger volume. The location of these zeros agree with Monte Carlo reweighting calculation for small volume. The application of the method for the O(2) model with a chemical potential is briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗