Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lattice Integrable Models”

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 325 records · Page 18

Many-body chaos in the antiferromagnetic quantum critical metal

Here, we compute the scrambling rate at the antiferromagnetic (AFM) quantum critical point, using the fixed point theory found by Schlief, Lunts, and Lee. At this strongly coupled fixed point, there is an emergent control parameter $w \ll 1$ that is a ratio of natural parameters of the theory. The strong coupling is unequally felt by the two degrees of freedom: the bosonic AFM collective mode is heavily dressed by interactions with the electrons, while the electron is only marginally renormalized. We find that the scrambling rates act as a measure of the “degree of integrability” of each sector of the theory: the Lyapunov exponent for the boson $λ^{(B)}_L \sim \mathcal{O} (\sqrt{w})k_BT/ℏ$ is significantly larger than the fermion one $λ^{(F)}_L \sim \mathcal{O}(w^2)k_BT/ℏ$, where $\textit{T}$ is the temperature. Although the interaction strength in the theory is of order unity, the larger Lyapunov exponent is still parametrically smaller than the universal upper bound of $λ_L = 2πk_BT/ℏ$. We also compute the spatial spread of chaos by the boson operator, whose low-energy propagator is highly nonlocal. We find that this nonlocality leads to a scrambled region that grows exponentially fast at intermediate distances, giving an infinite “butterfly velocity” of the chaos front, a result that has also been found in lattice models with long-range interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamic spin localization and γ -magnets

Here, we explore an unusual type of quantum matter that can be realized by qubits having different physical origins. Interactions in this matter are described by essentially different coupling operators for all qubits. We show that the simplest such models, which can be realized with localized states in Dirac materials, satisfy integrability conditions that we use to describe pseudospin dynamics in a linearly time-dependent magnetic field. Generalizing to an arbitrary number of qubits, we construct a spin Hamiltonian, which we call the γ -magnet. This system does not conserve polarization of any spin and the net spin polarization. Nevertheless, for arbitrarily strong interactions, nonadiabatic dynamics, and any initial eigenstate, we find that quantum interference suppresses spin flips. This behavior resembles many-body localization but occurs in phase space of many spins rather than real space. This effect may not have a counterpart in classical physics and can be a signature of a new type of spin ordering, which is different from both disordered spin glasses and ordered phases of spin lattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Parametric Deformation of Discrete Geometry for Aerodynamic Shape Design

We present a versatile discrete geometry manipulation platform for aerospace vehicle shape optimization. The platform is based on the geometry kernel of an open-source modeling tool called Blender and offers access to four parametric deformation techniques: lattice, cage-based, skeletal, and direct manipulation. Custom deformation methods are implemented as plugins, and the kernel is controlled through a scripting interface. Surface sensitivities are provided to support gradient-based optimization. The platform architecture allows the use of geometry pipelines, where multiple modelers are used in sequence, enabling manipulation difficult or impossible to achieve with a constructive modeler or deformer alone. We implement an intuitive custom deformation method in which a set of surface points serve as the design variables and user-specified constraints are intrinsically satisfied. We test our geometry platform on several design examples using an aerodynamic design framework based on Cartesian grids. We examine inverse airfoil design and shape matching and perform lift-constrained drag minimization on an airfoil with thickness constraints. A transport wing-fuselage integration problem demonstrates the approach in 3D. In a final example, our platform is pipelined with a constructive modeler to parabolically sweep a wingtip while applying a 1-G loading deformation across the wingspan. This work is an important first step towards the larger goal of leveraging the investment of the graphics industry to improve the state-of-the-art in aerospace geometry tools.

Anderson, George R.↗

Airframe Noise Simulations of a Full-Scale Large Civil Transport in Landing Configuration

This paper summarizes the results obtained from an extensive computational campaign to accurately predict full-scale landing gear noise for large civil transports. A highly accurate digital model of a full-scale Boeing 777-300ER aircraft with as-flown nose and main landing gear components was developed for use in the simulations. Two aircraft configurations were selected: nose and main landing gear deployed with wing high-lift devices retracted, and nose and main landing gear deployed with wing high-lift devices deflected. The two configurations were simulated without and with toboggan fairings installed on the main gear to represent the principal configurations evaluated during the 2005 QTD2 flight test. All simulations were performed with the lattice Boltzmann solver PowerFLOW® to resolve and capture the highly complex, unsteady flow field in the immediate vicinity of the aircraft. The far-field noise sig-nature of the aircraft was computed via a Ffowcs-Williams and Hawkings integral approach, with flow quantities on a permeable surface enclosing the source regions used as input. Synthetic pressure records at ground array microphone locations used during the QDT2 test were employed to generate narrowband acoustic maps and integrated far-field noise spectra. With high-lift devices retracted, the predicted spectra showed that landing gear noise is equivalent to total airframe noise, with no other airframe sources appearing within 10 dB of gear peak levels. Application of a toboggan fairing to the main gear produced modest noise reductions of 1-2 dB across the resolved frequency range. With high-lift devices deflected, undercarriage noise was within 3-4 dB of the total airframe noise, thus comprising nearly half of the total airframe noise. For this configuration, the toboggan fairing did not produce a reduction in noise, corroborating trends previously observed in QTD2 flight test data.

airframe noise↗

Characterization of thermally heat-treated polyacrylonitrile carbon fibers

This study investigates the graphitization process of polyacrylonitrile (PAN) carbon fibers by subjecting commercial fibers to thermal heat treatment at temperatures ranging from 1400 to 2100 °C in 100 °C increments, using either argon or nitrogen gas atmospheres. Changes in crystallinity, surface morphology, and lattice parameters were analyzed for two commercial carbon fibers using X-ray diffraction, scanning electron microscopy, and Raman spectroscopy. Results indicated minimal changes in surface morphology with increasing heat-treatment temperature; however, crystallinity significantly increased. Crystallinity changes were more strongly dependent on temperature rather than gas atmosphere or fiber type. At intermediate heat-treatment temperatures (1600–1800 °C), fibers treated in argon showed a slight preference for graphitization. The highest level of graphitization was measured at 2100 °C. Crystallite size increased as the intensity ratio of the D1 to G Raman peaks increased, reaching a peak around ~1800 °C, after which the ratio started to decrease. This behavior aligns with Ferrari's three-stage model of carbon crystallization and is consistent with both the Marie-Mering degree of graphitization and Brubaker's Integrated Absolute Differential models, all of which describe the transformation from an amorphous to a more graphitic structure. At the higher heat-treatment temperatures, the changes between atmospheres and fiber types were measured to converge to similar levels of graphitization. In conclusion, this study evaluates the progressive change in commercial grade carbon fibers when heat-treated.

Characterization↗

Foundations of variational discrete action theory

Variational wave functions and Green's functions are two important paradigms for solving quantum Hamiltonians, each having their own advantages. Here we detail the variational discrete action theory (VDAT), which exploits the advantages of both paradigms in order to approximately solve the ground state of quantum Hamiltonians. VDAT consists of two central components: the sequential product density matrix (SPD) ansatz and a discrete action associated with the SPD. The SPD is a variational ansatz inspired by the Trotter decomposition and characterized by an integer $\mathscr{N}$, recovering many well-known variational wave functions, in addition to the exact solution for $\mathscr{N}$ = ∞. The discrete action describes all dynamical information of an effective integer time evolution with respect to the SPD. We generalize the path integral to our integer time formalism, which converts a dynamic correlation function in integer time to a static correlation function in a compound space. We also generalize the usual many-body Green's function formalism to integer time, which results in analogous but distinct mathematical structures, yielding integer time versions of the generating functional, Dyson equation, and Bethe-Salpeter equation. We prove that the SPD can be exactly evaluated in the multiband Anderson impurity model (AIM) by summing a finite number of diagrams. For the multiband Hubbard model, we prove that the self-consistent canonical discrete action approximation (SCDA), which is the integer time analog of the dynamical mean-field theory, exactly evaluates the SPD for d = ∞. VDAT within the SCDA provides an efficient yet reliable method for capturing the local physics of quantum lattice models, which will have broad applications for strongly correlated electron materials. More generally, VDAT should find applications in various many-body problems in physics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Elucidating the Structural and Electronic Effects of Ni and Mn Cationic Incorporation on CoOOH for Efficient Benzyl Alcohol Electrooxidation

Transition-metal oxyhydroxides such as CoOOH are promising low-cost electrocatalysts for the selective electrooxidation of organic molecules, yet the influence of ubiquitous transition-metal impurities on their performance and durability remains poorly understood. Here, we experimentally probed the individual and synergistic electrochemical and structural effects of Ni and Mn incorporations into model CoOOH electrocatalysts toward an efficient benzyl alcohol oxidation reaction (BAOR). Comprehensive electrochemical, microscopic, and spectroscopic analyses reveal that Ni incorporation enhances charge-transfer kinetics and overall activity through the formation of catalytically active Ni 3+ sites, whereas Mn exhibited a more complex but interesting role. At the early stages of operation, Mn 4+ acts as a stabilizing surface layer that mitigates catalyst degradation but partially blocks Co sites before they undergo gradual leaching. The concurrent incorporation of both Ni and Mn yields a trimetallic 2NMC@NF electrocatalyst that integrates the activity benefits of Ni with the stability conferred by Mn, achieving 92.9% benzyl alcohol conversion and 91.4% Faradaic efficiency after 24 h at 1.5 V vs RHE. These findings elucidate how trace Ni and Mn impurities, often introduced from electrolytes or external sources, can modulate the lattice and electronic structure of CoOOH, offering a design strategy for enhancing both activity and long-term stability in electrocatalytic organic oxidation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides (Final Report)

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

30 DIRECT ENERGY CONVERSION↗

Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides. Final Report

The overarching goal of this research program, as originally delineated in the proposal “Tuning anisotropic bonding via chemistry and pressure in layered pnictides and chalcogenides” was to develop a predictive, chemistry-driven understanding of the impact of the phonon behavior on thermal properties of bulk layered materials. At finite temperatures, atomic vibrations (phonons) strongly impact the thermodynamics, thermal and electrical transport, and phase-switching properties of functional materials. In particular, soft phonon modes and strongly anharmonic potentials can have spectacular consequences, including structural phase transitions (for example in ferroelectrics and phase-change memory materials), metal-insulator transitions, and extreme thermal resistance preventing heat propagation. Bulk materials with highly-anisotropic bonding may provide unique strategies to induce soft-phonon modes and lattice instabilities. Recently, increasingly detailed investigations of the lattice dynamics in layered materials have been made possible by the advent of first-principles phonon calculations and advanced characterization techniques based on neutron and X-ray scattering. However, due to the lack of studies in which composition and bonding character are systematically varied, there are still fundamental questions regarding the impacts of anisotropic bonding and anharmonicity on lattice stability and thermal transport. One of the major goals of this research program is therefore to address this gap by coherently tuning bonding anisotropy and anharmonicity across families of related compounds. Such approaches have revealed new strategies for exploiting structural anisotropy in quasi-1D and 2D bulk materials to obtain tailored functional properties. This project systematically explored the lattice dynamics, phase stability, and transport properties in bulk layered materials by using both composition and applied pressure to tune the degree of bonding anisotropy and anharmonicity. To accomplish this work, we combined i) single-crystal growth of key material systems with tunable anisotropy, ii) in-situ high-temperature/high-pressure characterization of structure and phonons to probe bonding anisotropy and anharmonicity, including state-of-the-art inelastic X-ray scattering (IXS) and inelastic neutron scattering (INS), and iii) first-principles simulations leveraging large-scale computing to identify the fundamental origins of the observed effects, by relating atomic structure and dynamics to electronic orbital interactions. Finally, we modeled and verified the impact of the phonon behavior on thermal transport to identify new strategies for a-priori design of thermal conductivity. Our integrated collaborative approach helped to systematically unravel the effects of anisotropy and bonding anharmonicity on phonon transport, thermodynamics, and thermal properties of complex anisotropic materials.

36 MATERIALS SCIENCE↗

Magnetic Susceptibility and Neutron Scattering of Graphene in Antiferromagnetic State: a Tight-Binding Approach

We address here a tight-binding model study of the frequency-dependent antiferromagnetic spin susceptibility for the graphene systems. The Hamiltonian consists of electron hopping up to the third-nearest-neighbors, substrate and impurity effects in presence of electron-electron interactions at A and B sub-lattices. To calculate susceptibility, we evaluate the two-particle electron Green’s function by using Zubarev’s Green’s function technique. The frequency-dependent antiferromagnetic susceptibility of the system is computed numerically by taking 1000 X 1000 grid points of the electron momentum. The susceptibility displays a sharp peak at the neutron momentum transfer energy at low energies and another higher-energy peak associated with the substrate-induced gap. The evolution of these two peaks are investigated by varying neutron wave vector, Coulomb correlation energy, substrate-induced gap, electron hopping integrals and A- and B-site electron-doping concentrations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Origin of mean-field behavior in an elastic Ising model

Simple elastic models of spin-crossover compounds are known empirically to exhibit classical critical behavior. Here, we demonstrate how the long-range interactions responsible for this behavior arise naturally upon integrating out mechanical fluctuations of such a model. A mean-field theory applied to the resulting effective Hamiltonian quantitatively accounts for both thermodynamics and kinetics observed in computer simulations, including a barrier to magnetization reversal that grows extensively with system size. For nanocrystals, which break translational symmetry, a straightforward extension of mean-field theory yields similarly accurate results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semileptonic form factors for $$B\rightarrow D^*\ell \nu $$ at nonzero recoil from $$2+1$$-flavor lattice QCD: Fermilab Lattice and MILC Collaborations

Abstract We present the first unquenched lattice-QCD calculation of the form factors for the decay $$B\rightarrow D^*\ell \nu $$ B → D ∗ ℓ ν at nonzero recoil. Our analysis includes 15 MILC ensembles with $$N_f=2+1$$ N f = 2 + 1 flavors of asqtad sea quarks, with a strange quark mass close to its physical mass. The lattice spacings range from $$a\approx 0.15$$ a ≈ 0.15 fm down to 0.045 fm, while the ratio between the light- and the strange-quark masses ranges from 0.05 to 0.4. The valence b and c quarks are treated using the Wilson-clover action with the Fermilab interpretation, whereas the light sector employs asqtad staggered fermions. We extrapolate our results to the physical point in the continuum limit using rooted staggered heavy-light meson chiral perturbation theory. Then we apply a model-independent parametrization to extend the form factors to the full kinematic range. With this parametrization we perform a joint lattice-QCD/experiment fit using several experimental datasets to determine the CKM matrix element $$|V_{cb}|$$ | V cb | . We obtain $$\left| V_{cb}\right| = (38.40 \pm 0.68_{\text {th}} \pm 0.34_{\text {exp}} \pm 0.18_{\text {EM}})\times 10^{-3}$$ V cb = ( 38.40 ± 0 . 68 th ± 0 . 34 exp ± 0 . 18 EM ) × 10 - 3 . The first error is theoretical, the second comes from experiment and the last one includes electromagnetic and electroweak uncertainties, with an overall $$\chi ^2\text {/dof} = 126/84$$ χ 2 /dof = 126 / 84 , which illustrates the tensions between the experimental data sets, and between theory and experiment. This result is in agreement with previous exclusive determinations, but the tension with the inclusive determination remains. Finally, we integrate the differential decay rate obtained solely from lattice data to predict $$R(D^*) = 0.265 \pm 0.013$$ R ( D ∗ ) = 0.265 ± 0.013 , which confirms the current tension between theory and experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Mutual modulation between surface chemistry and bulk microstructure within secondary particles of nickel-rich layered oxides

Abstract Surface lattice reconstruction is commonly observed in nickel-rich layered oxide battery cathode materials, causing unsatisfactory high-voltage cycling performance. However, the interplay of the surface chemistry and the bulk microstructure remains largely unexplored due to the intrinsic structural complexity and the lack of integrated diagnostic tools for a thorough investigation at complementary length scales. Herein, by combining nano-resolution X-ray probes in both soft and hard X-ray regimes, we demonstrate correlative surface chemical mapping and bulk microstructure imaging over a single charged LiNi 0.8 Mn 0.1 Co 0.1 O 2 (NMC811) secondary particle. We reveal that the sub-particle regions with more micro cracks are associated with more severe surface degradation. A mechanism of mutual modulation between the surface chemistry and the bulk microstructure is formulated based on our experimental observations and finite element modeling. Such a surface-to-bulk reaction coupling effect is fundamentally important for the design of the next generation battery cathode materials.

25 ENERGY STORAGE↗

Interzeolite Transformation through Cross-Nucleation: A Molecular Mechanism for Seed-Assisted Synthesis

Polymorph selection and efficient crystallization are central goals in zeolite synthesis. Crystalline seeds are used for both purposes. While it has been proposed that zeolite seeds induce interzeolite transformation by dissolving into structural units that promote nucleation of the daughter crystal, the seed’s structural elements do not always match those of the target zeolite. This discrepancy raises the question of how the seed promotes the daughter phase. Here, we present the first molecularly resolved investigation of seed-assisted zeolite synthesis. Using molecular simulations, we reproduce the experimental finding that a parent zeolite can promote the nucleation of a daughter zeolite even when it lacks common composite building units (CBUs) or crystal planes. Modeling the seed-assisted synthesis of an AFI-type zeolite using zeolite CHA, our simulations indicate that stand-alone CBUs from the parent seed do not facilitate daughter crystal formation. However, introducing the intact seed significantly reduces the synthesis time, supporting that seed integrity is key to increased efficiency. This reduction arises from the cross-nucleation of the AFI-type zeolite on the CHA (001) face. We find that parent and daughter zeolites are connected by an interfacial transition layer with an order distinct from that of both zeolites. Simulations reveal that cross-nucleation occurs over a broad range of synthesis conditions. We argue that cross-nucleation would be most favorable for zeolite pairs that share crystalline planes such as those forming intergrowths. In conclusion, our findings suggest that the prevalence of intergrowths with a common lattice plane in zeolite synthesis is likely a kinetic effect of accelerated cross-nucleation.

Crystallization↗

Closeout of CRADA JSA 2012S004: Chapter 5, Integrated Control System, of the document of the ESS Conceptual Design Report, publicly available at https://europeanspallationsource.se/accelerator-documents

The integrated control system (ICS) is responsible for the whole ESS machine and facility: accelerator, target, neutron scattering instruments and conventional facilities. This unified approach keeps the costs of development, maintenance and support relatively low. ESS has selected a standardised, field-proven controls framework, the Experimental Physics and Industrial Control System (EPICS), which was originally developed jointly by Argonne and Los Alamos National Laboratories. Complementing this selection are best practices and experience from similar facilities regarding platform standardisation, control system development and device integration and commissioning. The components of ICS include the control system core, the control boxes, the BLED database management system, and the human machine interface. The control system core is a set of systems and tools that make it possible for the control system to provide required data, information and services to engineers, operators, physicists and the facility itself. The core components are the timing system that makes possible clock synchronisation across the facility, the machine protection system (MPS) and the personnel protection system (PPS) that prevent damage to the machine and personnel, and a set of control system services. Control boxes are servers that control a collection of equipment (for example a radio frequency cavity). The integrated control system will include many control boxes that can be assigned to one supplier, such as an internal team, a collaborating institute or a commercial vendor. This approach facilitates a clear division of responsibilities and makes integration much easier. A control box is composed of a standardised hardware platform, components, development tools and services. On the top level, it interfaces with the core control system components (timing, MPS, PPS) and with the human-machine interface. At the bottom, it interfaces with the equipment and parts of the facility through a set of analog and digital signals, real-time control loops and other communication buses. The ICS central data management system is named BLED (beam line element databases). BLED is a set of databases, tools and services that is used to store, manage and access data. It holds vital control system configuration and physics-related (lattice) information about the accelerator, target and instruments. It facilitates control system configuration by bringing together direct input-output controller (IOC) con guration and real-time data from proton and neutron beam line models. BLED also simplifies development and speeds up the code-test-debug cycle. The set of tools that access BLED will be tailored to the needs of different categories of users, such as ESS staff physicists, engineers, and operators; external partner laboratories; and visiting experimental instrument users. The human-machine interface is vital to providing a high-quality experience to ICS users. It encompasses a wide array of devices and software tools, from control room screens to engineer terminal windows; from beam physics data tools to post-mortem data analysis tools. It serves users with a wide range of skills from widely varied backgrounds. The Controls Group is developing a set of user profiles to accommodate this diverse range of use-cases and users.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

Direct growth of single-crystalline III–V semiconductors on amorphous substrates

The III–V compound semiconductors exhibit superb electronic and optoelectronic properties. Traditionally, closely lattice-matched epitaxial substrates have been required for the growth of high-quality single-crystal III–V thin films and patterned microstructures. To remove this materials constraint, here we introduce a growth mode that enables direct writing of single-crystalline III–V’s on amorphous substrates, thus further expanding their utility for various applications. The process utilizes templated liquid-phase crystal growth that results in user-tunable, patterned micro and nanostructures of single-crystalline III–V’s of up to tens of micrometres in lateral dimensions. InP is chosen as a model material system owing to its technological importance. The patterned InP single crystals are configured as high-performance transistors and photodetectors directly on amorphous SiO 2 growth substrates, with performance matching state-of-the-art epitaxially grown devices. In conclusion, the work presents an important advance towards universal integration of III–V’s on application-specific substrates by direct growth.

36 MATERIALS SCIENCE↗

Chiral matrix model of the semi-QGP in QCD

Previously, a matrix model of the region near the transition temperature, in the “semi”quark gluon plasma, was developed for the theory of SU(3) gluons without quarks. In this paper we develop a chiral matrix model applicable to QCD by including dynamical quarks with 2+1 flavors. This requires adding a nonet of scalar fields, with both parities, and coupling these to quarks through a Yukawa coupling, y. Treating the scalar fields in mean field approximation, the effective Lagrangian is computed by integrating out quarks to one loop order. As is standard, the potential for the scalar fields is chosen to be symmetric under the flavor symmetry of SU (3) L × SU(3) R × Z (3) A , except for a term linear in the current quark mass, m qk . In addition, at a nonzero temperature T it is necessary to add a new term, ~ m qk T 2 . The parameters of the gluon part of the matrix model are identical to those for the pure glue theory without quarks. The parameters in the chiral matrix model are fixed by the values, at zero temperature, of the pion decay constant and the masses of the pions, kaons, η , and η' . The temperature for the chiral crossover at T$χ$ = 155 MeV is determined by adjusting the Yukawa coupling y . We find reasonable agreement with the results of numerical simulations on the lattice for the pressure and related quantities. In the chiral limit, besides the divergence in the chiral susceptibility there is also a milder divergence in the susceptibility between the Polyakov loop and the chiral order parameter, with critical exponent β $-$ 1 . We compute derivatives with respect to a quark chemical potential to determine the susceptibilities for baryon number, the $χ$ 2n . Especially sensitive tests are provided by $χ$ 4 $-$ $χ$ 2 and by $χ$ 6 , which changes in sign about T$χ$ . In conclusion, the behavior of the susceptibilities in the chiral matrix model strongly suggests that as the temperature increases from T$χ$ , that the transition to deconfinement is significantly quicker than indicated by the measurements of the (renormalized) Polyakov loop on the lattice.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗