Search NASA⌕ Search

SEARCH · Search NASA

Results for “transfer matrix calculations”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Interference-Limited Absorption in Dense Molecular Nanolayers Near Reflecting Surfaces

We investigate linear resonant absorption by a dense ensemble of molecules confined to a sub-wavelength layer in two geometries: (i) a free-standing film in a homogeneous space and (ii) the same film placed at a controlled distance from a reflecting surface. In both cases, increasing the effective light–matter coupling (via molecular density/oscillator strength) produces a nonmonotonic response: absorption rises to an optimum and then decreases as the film becomes increasingly radiatively bright and reflective. Finite-difference time-domain simulations and analytical transfer-matrix calculations agree quantitatively and yield compact ridge conditions for the optimum. We interpret the trends using a scattering/port picture: the isolated film is a symmetric two-port system (reflection and transmission), which bounds single-sided resonant absorption to ≤50% in the ultrathin limit (reflecting transition saturation), whereas adding a mirror suppresses transmission and converts the structure into an effectively one-port absorber. In the mirror-backed geometry, interference can cancel reflection, and unity absorption is obtained at critical coupling, when radiative leakage is balanced by intrinsic molecular loss. These results clarify fundamental limits and design rules for collective absorption in dense molecular layers near dielectric or metallic boundaries.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Spontaneous Formation of Single-Crystalline Spherulites in a Chiral 2D Hybrid Perovskite

In two-dimensional (2D) chiral metal-halide perovskites (MHPs), chiral organic spacers induce structural chirality and chiroptical properties in the metal-halide sublattice. This structural chirality enables reversible crystalline-glass phase transitions in (S-NEA) 2 PbBr 4 , a prototypical chiral 2D MHP where NEA + represents 1-(1-naphthyl)ethylammonium. Here, in this study, we investigate two distinct spherulite states of (S-NEA) 2 PbBr 4 , exhibiting either radial-like or stripe-like banded patterns depending on the annealing conditions of the amorphous film. Despite similarities in optical absorption and photoluminescence, the stripe-like, banded spherulite exhibits higher crystallinity and improved optical transparency compared to those of radial-like spherulite. X-ray nanoprobe measurements reveal tilting-angle modulations in the octahedral plane of stripe-like spherulites, correlating with the film’s surface geometry. Transfer matrix calculations indicate that the optical contrast in stripe-like patterns, seen in bright-field optical microscopy, arises from optical interference effects, differing from the contrast mechanism observed in polymer spherulites. Ultrafast carrier dynamics experiments suggest that the stripe-like spherulites resemble single crystals more closely than radial-like spherulites, while electrical conductivity measurements show enhanced charge carrier transport in stripe-like spherulites. These findings offer insights into MHP spherulite states with a single composition but different morphologies, previously observed only in polymers, highlighting their potential for optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Lanczos algorithm for lattice QCD matrix elements

Recent work [M. L. Wagman, Lanczos, the transfer matrix, and the signal-to-noise problem, .] found that an analysis formalism based on the Lanczos algorithm allows energy levels to be extracted from Euclidean correlation functions with faster ground-state convergence than effective masses, convergent estimators for multiple states from a single correlator, and two-sided error bounds. After filtering out spurious eigenvalues and using outlier-robust estimators within a nested bootstrap framework, Lanczos estimators behave more like multistate fit results than effective masses—but without involving statistical fitting. We extend this formalism to the determination of matrix elements from three-point correlation functions and provide a physical picture of “spurious-state filtering” involving restriction to a Hermitian subspace. We demonstrate similar advantages for matrix elements as for spectroscopy through example applications to noiseless mock-data and (bare) forward matrix elements of the strange scalar current between both ground and excited states with the quantum numbers of the nucleon.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Disordered hydrogen adsorption at the three-fold site of W(110)

Hydrogen often forms disordered phases on metals; however, deciphering the atomic-scale behavior requires experimental and modeling techniques that directly probe the short-range order between neighboring hydrogen adsorbates. Here, to address this challenge, we applied direct recoil spectroscopy (DRS) to investigate hydrogen adsorption on the W(110) surface. We show that the recoiled hydrogen flux measured during DRS is sensitive to the short range order of the adsorbed hydrogen. Ordered hydrogen neighbors are much more efficient at dechanneling incident ions along low-index surface channels, strongly suppressing the recoiled hydrogen flux along these directions. By modeling the DRS measurements with molecular dynamics for ordered and disordered hydrogen phases, we find that at a hydrogen surface coverage of approximately Θ = 0.6, the hydrogen adsorbates are bound to three-fold sites in a disordered phase at room temperature. This finding is consistent with transfer-matrix scaling model calculations that suggest the transition to an ordered hydrogen phase occurs at a greater surface coverage of Θ = 0.83.

08 HYDROGEN↗
DOE OSTI2025

Generalized parton distributions from lattice QCD with asymmetric momentum transfer: Tensor case

The calculation of generalized parton distributions (GPDs) in lattice QCD was traditionally done by calculating matrix elements in the symmetric frame. Recent advancements have significantly reduced computational costs by calculating these matrix elements in the asymmetric frame, allowing us to choose the momentum transfer to be in either the initial or final states only. The theoretical methodology requires a new parametrization of the matrix element to obtain Lorentz-invariant amplitudes, which are then related to the GPDs. The formulation and implementation of this approachaveh already been established for the unpolarized and helicity GPDs. Building upon this idea, we extend this formulation to the four leading-twist quark transversity GPDs ($𝐻_𝑇$, $𝐸_𝑇$, $\tilde{𝐻}_𝑇$, $\tilde{𝐸}_𝑇$). We also present numerical results for zero skewness using an 𝑁 𝑓 = 2 +1 +1 ensemble of twisted mass fermions with a clover improvement. The light quark masses employed in these calculations correspond to a pion mass of about 260 MeV. Furthermore, we include a comparison between the symmetric and asymmetric frame calculations to demonstrate frame independence of the Lorentz-invariant amplitudes. Analysis of the matrix elements in the asymmetric frame is performed at several values of the momentum transfer squared, −𝑡, ranging from 0.17 to 2.29 GeV 2 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Terahertz Dirac Hyperbolic Metamaterial

Hyperbolic metamaterials (HMMs) are engineered materials with a hyperbolic isofrequency surface, enabling a range of interesting phenomena and applications including negative refraction, enhanced sensing, and subdiffraction imaging, focusing, and waveguiding. Existing HMMs primarily work in the visible and infrared spectral range due to the inherent properties of their constituent materials. Here, we demonstrate a THz-range Dirac HMM using topological insulators as the building blocks. Here, we find that the structure houses up to three high-wavevector volume plasmon polariton (VPP) modes, consistent with transfer matrix modeling and effective medium theory calculations. The VPPs have mode indices greater than 100, significantly larger than observed for VPP modes in HMMs made from metals or doped semiconductors while maintaining comparable quality factors. We attribute these properties to the two-dimensional Dirac nature of the electrons occupying the topological insulator surface states. Because these are van der Waals materials, these structures can be grown at a wafer-scale on a variety of substrates, allowing them to be integrated with existing THz structures and enabling next-generation THz optical devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Three-dimensional imaging of pion using lattice QCD: generalized parton distributions

In this work, we report a lattice calculation of x-dependent valence pion generalized parton distributions (GPDs) at zero skewness with multiple values of the momentum transfer −t. The calculations are based on an N f = 2 + 1 gauge ensemble of highly improved staggered quarks with Wilson-Clover valence fermion. The lattice spacing is 0.04 fm, and the pion valence mass is tuned to be 300 MeV. We determine the Lorentz-invariant amplitudes of the quasi-GPD matrix elements for both symmetric and asymmetric momenta transfers with similar values and show the equivalence of both frames. Then, focusing on the asymmetric frame, we utilize a hybrid scheme to renormalize the quasi-GPD matrix elements obtained from the lattice calculations. After the Fourier transforms, the quasi-GPDs are then matched to the light-cone GPDs within the framework of large momentum effective theory with improved matching, including the next-to-next-to-leading order perturbative corrections, and leading renormalon and renormalization group resummations. We also present the 3-dimensional image of the pion in impact-parameter space through the Fourier transform of the momentum transfer −t.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Moments of axial-vector GPD from lattice QCD: quark helicity, orbital angular momentum, and spin-orbit correlation

In this work, we present a lattice QCD calculation of the Mellin moments of the twist-2 axial-vector generalized parton distribution (GPD), $\overset{\sim }{H}\left(x,\xi, t\right)$ , at zero skewness, ξ, with multiple values of the momentum transfer, t. Our analysis employs the short-distance factorization framework on ratio-scheme renormalized quasi-GPD matrix elements. The calculations are based on an N f = 2 + 1 + 1 twisted mass fermions ensemble with clover improvement, a lattice spacing of a = 0.093 fm, and a pion mass of m π = 260 MeV. We consider both the iso-vector and iso-scalar cases, utilizing next-to-leading-order perturbative matching while omitting the disconnected contributions and gluon mixing in the iso-scalar case. For the first time, we determine the Mellin moments of $\overset{\sim }{H}$ up to the fifth order. From these moments, we discuss the quark helicity and orbital angular momentum contributions to the nucleon spin, as well as the spin-orbit correlations of the quarks. Additionally, we perform a Fourier transform over the momentum transfer, which allows us to explore the spin structure in the impact-parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Multi-stage load partitioning in additively manufactured Al6061+TiC nanocomposite characterized by in-situ neutron diffraction

Metal matrix composites (MMCs), combining metal matrix and ceramic particles, exhibit high mechanical strength and stiffness compared to conventional metallic materials. During fusion-based additive manufacturing (AM), MMCs undergo a rapid heating and cooling thermal history, which affects the interfacial bonding, thermal misfit stress and mechanical load transfer behavior between the metal matrix and particles, thus impacting the bulk mechanical performance. Here, we investigate the deformation dynamics and phase-specific load transfer in fusion-based AMed Al6061 +TiC MMCs through in-situ neutron diffraction. By calculating the phase-specific lattice strains, a multi-stage load transfer and deformation behavior during uniaxial compression is revealed: (1) elastic deformation in Stage I for both phases; (2) sudden stress rebalance between two phases in Stage II, evidenced by a sudden decrease in Al lattice strain and increase in TiC; (3) active load carrying by both phases in Stage III, with both phases experiencing an increase in lattice strains. Further, the deformation mechanism of each stage is deducted by correlating the evolutions of the interphase stresses and peaks’ broadening and intensity of the Al matrix. The local plastic deformation of the Al matrix near the phase interface is triggered and leads to stress rebalancing in Stage II. The global plastic deformation subsequently propagates throughout the Al matrix in Stage III, and the composite is further hardened through both dislocation multiplication and load transfer. The findings offer valuable insights into deformation and load-sharing behavior in AMed MMCs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An FPGA-based hardware accelerator supporting sensitive sequence homology filtering with profile hidden Markov models

Abstract Background Sequence alignment lies at the heart of genome sequence annotation. While the BLAST suite of alignment tools has long held an important role in alignment-based sequence database search, greater sensitivity is achieved through the use of profile hidden Markov models (pHMMs). Here, we describe an FPGA hardware accelerator, called HAVAC, that targets a key bottleneck step (SSV) in the analysis pipeline of the popular pHMM alignment tool, HMMER. Results The HAVAC kernel calculates the SSV matrix at 1739 GCUPS on a $$\sim$$ ∼ $3000 Xilinx Alveo U50 FPGA accelerator card, $$\sim$$ ∼ 227× faster than the optimized SSV implementation in nhmmer . Accounting for PCI-e data transfer data processing, HAVAC is 65× faster than nhmmer’s SSV with one thread and 35× faster than nhmmer with four threads, and uses $$\sim$$ ∼ 31% the energy of a traditional high end Intel CPU. Conclusions HAVAC demonstrates the potential offered by FPGA hardware accelerators to produce dramatic speed gains in sequence annotation and related bioinformatics applications. Because these computations are performed on a co-processor, the host CPU remains free to simultaneously compute other aspects of the analysis pipeline.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Photon-odderon interference in exclusive χ c charmonium production at the Electron-Ion Collider

Exclusive C = + 1 scalar, axial-vector, and tensor quarkonium production in high-energy electron-proton scattering requires a C -odd t -channel exchange of a photon or a three-gluon ladder. We derive the expressions for the corresponding amplitudes. The relative phase of the photon vs three-gluon exchange amplitudes is determined by the sign of the light-front matrix element of the eikonal color current operator d a b c J + a J + b J + c at moderate x , and is not affected by small- x QCD evolution. Model calculations predict constructive interference, which is particularly strong for momentum transfer | t | ∼ 1 GeV 2 where the cross section for χ c J production exceeds that for pure photon exchange by up to a factor of 4. Exclusive χ c J electroproduction at the high luminosity Electron-Ion Collider occurs with well measurable rates and measurements of these processes should find evidence for the perturbative odderon exchange. We also compute the total electroproduction cross section as a function of energy and provide first estimates of the number of χ c J events per month at the Electron-Ion Collider design luminosity. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Skewness-dependent moments of the pion GPD from nonlocal quark-bilinear correlators

We present lattice QCD calculations of the odd Mellin moments of pion valence-quark generalized parton distribution up to fifth order ⟨𝑥 4 ⟩ and for the skewness range [−0.33, 0] using operator product expansion of bilocal quark-bilinear operators. The calculations are performed on an ensemble with lattice spacing 𝑎 = 0.04 fm and valence pion mass 300 MeV, employing boosted pion states with momenta up to 2.428 GeV and momentum transfers reaching 2.748 GeV 2 . We employ ratio-scheme renormalization and next-to-leading logarithmic resummed perturbative matching. At zero skewness, our results are consistent with previous lattice studies. By combining matrix elements at multiple values of skewness and momentum transfer, skewness-dependent moments are obtained through simultaneous polynomiality-constrained fits.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Unifying Combinatorial and Graphical Methods in Artificial Intelligence

Recently, a new graph Laplacian, called the inner product Laplacian, was introduced which generalizes many existing Laplacians, including the normalized and combinatorial Laplacian and their weighted variants. The key observation behind the inner product Laplacian is that by defining appropriate inner product spaces on the vertices and edges, the standard Laplacians can be recovered as Hodge Laplacians over the simplicial complex formed by the edges and vertices. These inner product spaces form a natural way to incorporate non-combinatorial information into the definition of a domain-specific Laplacian. In particular, in contrast to current domain-specific weighting schemes which rely solely on edge weights, information regarding the similarity of non-adjacent vertices and arbitrary pairs of edges can be effectively incorporated into the Laplacian. In order to illustrate this approach we consider the problem of calculating the potential energy of an atomistic configuration using Graph Neural Networks. In comparison with start-of-the-art approaches, such as SchNet, our approach replaces a learned (via auto-encoder) representation of the atom types with an inner product space on atoms based on scientific knowledge (e.g., electronegativity). We will illustrate how this approach captures key chemical properties of the molecules and compare the energy calculations with state-of-the-art neural network approaches. However, to compute the resulting Laplacian involves a mixture of sparse and dense matrix computation and yields a dense matrix as the basis for the graph convolution. This dense convolutional kernel necessitates moving away from the standard message passing framework for graph neural networks and increases the computational cost of applying the kernel. In order to mitigate these costs we investigate means of leveraging the mixed sparse and dense computations to reduce the overall computational cost and how these approaches can be automatically transferred to energy efficient hardware (e.g., field programmable gate arrays (FPGAs)).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Mechanistic Insights into Visible-Light-Induced ATRA Reactions Powered by the Symbiotic Relationship between Cu(II)/Cu(I) Phenanthroline Complexes

As photoredox catalysis continues to yield promising chemical transformations, there is an increased need to understand how specific photocatalysts function to improve reaction efficiencies while expanding their scope. Copper-phenanthroline-based photocatalysts such as Cu II (dap)Cl 2 (dap = 2,9-di(p-anisyl)-1,10-phenanthroline) and [Cu I (dap) 2 ]Cl were both found to be equally capable of olefin activation through electrophilic atom transfer radical addition (ATRA) reactions. Although these molecular catalysts have proven successful, many intermediates suggested in the proposed catalytic cycles have never been detected. One undetermined aspect in this chemistry is related to how one equivalent of Cu II (dap)Cl 2 generates half of an equivalent of [Cu I (dap) 2 ] + during the photocatalytic sequence. To this end, we initially used more synthetically accessible model systems, namely, [Cu I (dpp) 2 ]Cl and Cu II (dpp)Cl 2 (dpp = 2,9-diphenyl-1,10-phenanthroline), to glean detailed mechanistic insights into this unusual symbiotic relationship. We directly detected several intermediates involved in the ATRA photocatalytic cycle using these model chromophores in conjunction with electronic spectroscopy, infrared matrix-assisted laser desorption electrospray ionization (IR-MALDESI) mass spectrometry, electronic structure calculations, EPR spin-trap experiments, and 1 H NMR spectroscopy. We found that the unique ligand lability and coordinating properties of acetonitrile enable both the in situ oxidation of [Cu I (dpp) 2 ] + by tosyl chloride into Cu II (dpp)Cl 2 and the visible-light-induced homolysis of the Cu II –Cl bond, which initiates the conversion to the Cu I species [Cu I (dpp) 2 ][Cu I Cl 2 ]. Here, the combined findings from the present study of the catalytic cycle demonstrate that the symbiotic relationship between Cu II (dpp)Cl 2 and [Cu I (dpp) 2 ] + , as well as between Cu II (dap)Cl 2 and [Cu I (dap) 2 ] + , is the critical factor enabling the ATRA photoreaction by departing from either photocatalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2025

𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ Resonance Form Factors from 𝐵→ 𝜋⁢𝜋⁢ℓ⁢$\bar{v}$ in Lattice QCD

The decay 𝐵 → 𝜌⁢ℓ⁢$\bar{v}$ is an attractive process for determining the magnitude of the smallest Cabibbo-Kobayashi-Maskawa matrix element, |𝑉 𝑢⁢𝑏 |, and can provide new insights into the origin of the long-standing exclusive-inclusive discrepancy in determinations of this standard-model parameter. This requires a nonperturbative QCD calculation of the 𝐵 → 𝜌 form factors 𝑉, 𝐴 0 , 𝐴 1 , and 𝐴 12 . The unstable nature of the 𝜌 resonance has prevented precise lattice QCD calculations of these form factors to date. Here, we present the first lattice QCD calculation of the 𝐵 → 𝜌 form factors in which the 𝜌 is treated properly as a resonance in 𝑃-wave 𝜋⁢𝜋 scattering. To this end, we use the Lellouch-Lüscher finite-volume formalism to compute the 𝐵 → 𝜋⁢𝜋 form factors as a function of both momentum transfer and 𝜋⁢𝜋 invariant mass, and then analytically continue to the 𝜌 resonance pole. This calculation is performed with 2 + 1 dynamical quark flavors at a pion mass of approximately 320 MeV, and demonstrates a clear path toward results at the physical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗