Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase field method”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Complexity of many-body interactions in transition metals via machine-learned force fields from the TM23 data set

Abstract This work examines challenges associated with the accuracy of machine-learned force fields (MLFFs) for bulk solid and liquid phases ofd-block elements. In exhaustive detail, we contrast the performance of force, energy, and stress predictions across the transition metals for two leading MLFF models: a kernel-based atomic cluster expansion method implemented using sparse Gaussian processes (FLARE), and an equivariant message-passing neural network (NequIP). Early transition metals present higher relative errors and are more difficult to learn relative to late platinum- and coinage-group elements, and this trend persists across model architectures. Trends in complexity of interatomic interactions for different metals are revealed via comparison of the performance of representations with different many-body order and angular resolution. Using arguments based on perturbation theory on the occupied and unoccupieddstates near the Fermi level, we determine that the large, sharpddensity of states both above and below the Fermi level in early transition metals leads to a more complex, harder-to-learn potential energy surface for these metals. Increasing the fictitious electronic temperature (smearing) modifies the angular sensitivity of forces and makes the early transition metal forces easier to learn. This work illustrates challenges in capturing intricate properties of metallic bonding with current leading MLFFs and provides a reference data set for transition metals, aimed at benchmarking the accuracy and improving the development of emerging machine-learned approximations.

Chemistry↗
DOE OSTI2022

First-principles calculation of gate-tunable ferromagnetism in magic-angle twisted bilayer graphene under pressure

Magic-angle twisted bilayer graphene (MATBG) is notable as a highly tunable platform for investigating strongly correlated phenomena such as unconventional superconductivity and quantum spin liquids, due to easy control of doping level through gating and sensitive dependence of the magic angle on hydrostatic pressure. Experimental observations of correlated insulating states, unconventional superconductivity and ferromagnetism in MATBG indicate that this system exhibits rich exotic phases. In this work, using density functional theory calculations in conjunction with the effective screening medium method, we find the MATBG under pressure at a twisting angle of 2.88 °and simulate how its electronic states evolve when doping level and electric field perpendicular to plane are tuned by gating. Our calculations show that, at doping levels between two electrons and four holes per moiré unit cell, a ferromagnetic (FM) solution with spin density localized at AA stacking sites is lower in energy than the nonmagnetic solution. The magnetic moment of this FM state decreases with both electron and hole doping and vanishes at four electrons/holes doped per moiré unit cell. Hybridization between the flat bands at the Fermi level and the surrounding dispersive bands can take place at finite doping. On increasing the out-of-plane electric field at zero doping, a transition from the FM state to the nonmagnetic one is seen. An investigation of impurity effects shows that both absorption of O 2 molecules and occurrence of Stone–Wales impurities suppress the FM state, and the mechanisms are understood from our calculations. We also analyze the interlayer bonding character due to flat bands via Wannier functions. Finally, we report trivial band topology of the flat bands in the FM state at a certain doping level.

Physics↗
DOE OSTI2023

Isotope effect and critical magnetic fields of superconducting YH 6 : A Migdal-Eliashberg theory approach

The emergence of near-ambient temperature superconductivity under pressure in metal hydride systems has motivated a desire to further understand such remarkable properties, specifically critical magnetic fields. Here, YH 6 is suggested to be a departure from conventional superconductivity, due to apparent anomalous behavior. Using density functional calculations in conjunction with Migdal-Eliashberg theory we show that in YH 6 the critical temperature and the isotope effect under pressure, as well as the high critical fields, are consistent with strong-coupling conventional superconductivity, a property anticipated to extend to other related systems. Furthermore, strong-coupling corrections occur to the expected BCS values for the isotope effect coefficient (α), Ginzburg-Landau parameter [κ 1 ⁡(T)], London penetration depth [λ L (T)], electromagnetic coherence length [ξ⁡(T)], and the energy gap (Δ 0 ).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Group-equivariant autoencoder for identifying spontaneously broken symmetries

We introduce the group-equivariant autoencoder (GE autoencoder), a deep neural network (DNN) method that locates phase boundaries by determining which symmetries of the Hamiltonian have spontaneously broken at each temperature. We use group theory to deduce which symmetries of the system remain intact in all phases, and then use this information to constrain the parameters of the GE autoencoder such that the encoder learns an order parameter invariant to these “never-broken” symmetries. This procedure produces a dramatic reduction in the number of free parameters such that the GE-autoencoder size is independent of the system size. We include symmetry regularization terms in the loss function of the GE autoencoder so that the learned order parameter is also equivariant to the remaining symmetries of the system. By examining the group representation by which the learned order parameter transforms, we are then able to extract information about the associated spontaneous symmetry breaking. We test the GE autoencoder on the 2D classical ferromagnetic and antiferromagnetic Ising models, finding that the GE autoencoder (1) accurately determines which symmetries have spontaneously broken at each temperature; (2) estimates the critical temperature in the thermodynamic limit with greater accuracy, robustness, and time efficiency than a symmetry-agnostic baseline autoencoder; and (3) detects the presence of an external symmetry-breaking magnetic field with greater sensitivity than the baseline method. Lastly, we describe various key implementation details, including a quadratic-programming-based method for extracting the critical temperature estimate from trained autoencoders and calculations of the DNN initialization and learning rate settings required for fair model comparisons.

42 ENGINEERING↗
DOE OSTI2022

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Adaptive sampling for accelerating neutron diffraction-based strain mapping *

Abstract Neutron diffraction is a useful technique for mapping residual strains in dense metal objects. The technique works by placing an object in the path of a neutron beam, measuring the diffracted signals and inferring the local lattice strain values from the measurement. In order to map the strains across the entire object, the object is stepped one position at a time in the path of the neutron beam, typically in raster order, and at each position a strain value is estimated. Typical dwell times at neutron diffraction instruments result in an overall measurement that can take several hours to map an object that is several tens of centimeters in each dimension at a resolution of a few millimeters, during which the end users do not have an estimate of the global strain features and are at risk of incomplete information in case of instruments outages. In this paper, we propose an object adaptive sampling strategy to measure the significant points first. We start with a small initial uniform set of measurement points across the object to be mapped, compute the strain in those positions and use a machine learning technique to predict the next position to measure in the object. Specifically, we use a Bayesian optimization based on a Gaussian process regression method to infer the underlying strain field from a sparse set of measurements and predict the next most informative positions to measure based on estimates of the mean and variance in the strain fields estimated from the previously measured points. We demonstrate our real-time measure-infer-predict workflow on additively manufactured steel parts—demonstrating that we can get an accurate strain estimate even with 30%–40% of the typical number of measurements—leading the path to faster strain mapping with useful real-time feedback. We emphasize that the proposed method is general and can be used for fast mapping of other material properties such as phase fractions from time-consuming point-wise neutron measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Impact of the plastic deformation microstructure in metals on the kinetics of recrystallization: A phase-field study

The sensitivity of recrystallization kinetics in metals to the heterogeneity of microstructure and deformation history is a widely accepted experimental fact. However, most of the available recrystallization models employ either a mean field approach or use grain-averaged parameters, and thus neglecting the mesoscopic heterogeneity induced by prior deformation. Here in the present study, we investigate the impact of deformation-induced dislocation (subgrain) structure on the kinetics of recrystallization in metals using the phase-field approach. The primary focus here is upon the role of dislocation cell boundaries. The free energy formulation of the phase-field model accounts for the heterogeneity of the microstructure by assigning localized energy to the resulting dislocation microstructure realizations generated from experimental data. These microstructure realizations are created using the universal scaling laws for the spacing and the misorientation angles of both the geometrically necessary and incidental dislocation boundaries. The resulting free energy is used into an Allen-Cahn based model of recrystallization kinetics, which are solved using the finite element method. The solutions thus obtained shed light on the critical role of the spatial heterogeneity of deformation in the non-smooth growth of recrystallization nuclei and on the final grain structure. The results showed that, in agreement with experiment, the morphology of recrystallization front exhibits protrusions and retrusions. By resolving the subgrain structure, the presented algorithm paves the way for developing predictive kinetic models that fully account for the deformed state of recrystallizing metals.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Surrogate Model Integration with MOOSE XFEM for Creep Crack Growth

Ferritic-martensitic steels are key structural materials for advanced reactors but experience time-dependent deformation and damage under prolonged high temperature and irradiation, leading to creep-driven crack initiation and growth. High-fidelity models—crystal plasticity with irradiation mechanisms, phase-field for microstructural evolution, and continuum-damage viscoplasticity—capture the underlying physics but are too computationally intensive for broad design-space exploration and uncertainty quantification. This milestone advances a scalable alternative by integrating a microstructure-sensitive surrogate creep model into the Multiphysics Object-Oriented Simulation Environment (MOOSE) finite element framework and extending it to fracture via the extended finite element method (XFEM). The surrogate model, developed with collaborators at Sandia and Los Alamos National Laboratories, maps relevant microstructural descriptors to the viscoplastic response of HT9. We embed this surrogate within a coupled deformation-damage workflow in MOOSE/XFEM to simulate creep-driven crack initiation and propagation. Implementation enhancements include updates to the material interface, a plastic correction phase involving microstructure evolution, and fracture criteria to ensure numerical robustness and compatibility with the surrogate structure. Demonstrations on canonical creep benchmarks spanning uniaxial and multiaxial states show that the surrogate reproduces key trends of high-fidelity models while substantially reducing computational cost. The resulting capability bridges physics fidelity and performance, providing a practical path to a predictive, microstructure-aware assessment of creep and fracture in reactor materials.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

A systematic study and framework of fringe projection profilometry with improved measurement performance for in-situ LPBF process monitoring

Fringe Projection Profilometry (FPP) is a cost-effective and non-invasive technology that has been shown to measure finer features. Here, in this work, we developed an in-situ FPP method to measure the dynamic topography of powder bed and printed layer during Laser Powder Bed Fusion (LPBF) additive manufacturing (AM) process. A systematic study towards developing a comprehensive framework of LPBF-specific FPP is demonstrated to enhance and evaluate the performance of applying FPP for in-situ LPBF monitoring, including 1) a modified sensor model with localized correction; 2) improved phase unwrapping with FFT filtering 3) quantitative uncertainty analysis; and 4) experimental validation with ex-situ characterization. The developed LPBF-specific FPP system and methods are implemented on a commercial LPBF-AM machine, achieving better accuracy, more robustness, and increased field of view while maintaining sufficient measurement range and decent resolution, in contrast to literature methods. The established FPP framework will facilitate the development of closed-loop control strategies for advancing LPBF based AM.

42 ENGINEERING↗
DOE OSTI2024

A Frequency Domain Methodology for Quantitative Evaluation of Diffuse Wavefield With Applications to Seismic Imaging

Abstract Ambient Noise Imaging (ANI) of subsurface structures relies on seismic interferometry of diffuse seismic wavefields. However, the lack of effective methods to quantify and identify highly diffuse waves hampers applications of ANI, particularly in evaluating seismic attenuation and monitoring structural changes with high temporal resolution. Conventional ANI approaches require data normalization, which effectively suppresses the non‐diffuse component with large amplitude but also results in significant loss of amplitude and phase information in the continuous seismic records. In this study, we propose a frequency domain method to quantitatively evaluate the degree of diffuseness of seismic wavefields by analyzing their statistical characteristics of modal amplitudes for stationarity and randomness. Tests on synthetic waveform and field nodal records show that the proposed method can effectively distinguish between diffuse and non‐diffuse waveforms for either single‐ or three‐component data. As an application, we identify a 60‐s‐long diffuse coda of a local M 2.2 earthquake recorded by a dense nodal array on the San Jacinto Fault Zone, and successfully extract high‐quality dispersion curve andQ‐value without performing data normalization. These results are consistent with those obtained by conventional methods that assess the correlation between coherency and the Green's function, and by modeling ballistic waves generated by road traffic. Our proposed method can advance the imaging of subsurface velocity and attenuation structures as well as monitoring temporal changes for scientific studies and engineering applications.

Geochemistry & Geophysics↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2024

Extracting Higher Central Charge from a Single Wave Function

A (2+1)D topologically ordered phase may or may not have a gappable edge, even if its chiral central charge c_ is vanishing. Recently, it was discovered that a quantity regarded as a "higher" version of chiral central charge gives a further obstruction beyond c_ to gapping out the edge. Here, in this Letter, we show that the higher central charges can be characterized by the expectation value of the partial rotation operator acting on the wave function of the topologically ordered state. This allows us to extract the higher central charge from a single wave function, which can be evaluated on a quantum computer. Our characterization of the higher central charge is analytically derived from the modular properties of edge conformal field theory, as well as the numerical results with the ν = 1/2 bosonic Laughlin state and the non-Abelian gapped phase of the Kitaev honeycomb model, which corresponds to U(1) 2 and Ising topological order, respectively. The Letter establishes a numerical method to obtain a set of obstructions to the gappable edge of (2+1)D bosonic topological order beyond c_, which enables us to completely determine if a (2+1)D bosonic Abelian topological order has a gappable edge or not. We also point out that the expectation values of the partial rotation on a single wave function put a constraint on the low-energy spectrum of the bulk-boundary system of (2+1)D bosonic topological order, reminiscent of the Lieb-Schultz-Mattis-type theorems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Development of a Beam Halo Monitor

Muons, Inc. proposed improved Residual-Gas luminescent Beam Profile Monitors for Intense Beams by adding novel features including a mirascope, a light collecting device devised in the 16th century by an inventor of the telescope. A modern mirascope consists of 2 paraboloids. Light generated by the beam is collected by the paraboloids and creates a beam image. By placing a lens near the image it is possible to collect almost 100% of the emitted light, up to ten times better than existing fluorescence profilometers. Intense (multi-MW) proton beams pose a challenge for beam profile measurements. Solid-based beam intercepting instrumentation produces unallowable levels of radiation at these high powers. An alternative is to use a zero-or-low mass device such as a residual-gas profile monitor, either through collection of ionization or fluorescence. Fluorescence has an intrinsic advantage over the usual collection of ionized electrons in that it is not affected by the strong and variable electric and magnetic fields created by the beam that can distort profiles and make calibration difficult. For this reason, we proposed the fluorescence technique using the mirascope and other methods to enhance light collection and versatility by modifying the residual gas and its spatial distribution. In Phase I, we designed a mirascope and its electronics for commissioning and beam tests in a proposed Phase II. The design was based on optical simulations done in Phase I, and a physics run plan was developed for optimizing and testing. The Phase II experimental site and program have been planned and our goal is to produce a universal low cost Mirascope Residual-Gas Luminescent Beam Profile Monitor that can replace all other profile monitors currently used in accelerators.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Quantum simulation of $\phi^4$ theories in qudit systems

We discuss the implementation of quantum algorithms for lattice $\Phi^4$ theory on circuit quantum electrodynamics (cQED) system. The field is represented on qudits in a discretized field amplitude basis. The main advantage of qudit systems is that its multi-level characteristic allows the field interaction to be implemented only with diagonal single-qudit gates. Considering the set of universal gates formed by the single-qudit phase gate and the displacement gate, we address initial state preparation and single-qudit gate synthesis with variational methods.

Kurkcuoglu, Doga Murat↗
DOE OSTI2024

Identification of Kelvin-Helmholtz generated vortices in magnetised fluids

The Kelvin-Helmholtz Instability (KHI), arising from velocity shear across the magnetopause, plays a significant role in the viscous-like transfer of mass, momentum, and energy from the shocked solar wind into the magnetosphere. While the KHI leads to growth of surface waves and vortices, suitable detection methods for these applicable to magnetohydrodynamics (MHD) are currently lacking. A novel method is derived based on the well-established λ-family of hydrodynamic vortex identification techniques, which define a vortex as a local minimum in an adapted pressure field. The J × B Lorentz force is incorporated into this method by using an effective total pressure in MHD, including both magnetic pressure and a pressure-like part of the magnetic tension derived from a Helmholtz decomposition. The λ MHD method is shown to comprise of four physical effects: vortical momentum, density gradients, fluid compressibility, and the rotational part of the magnetic tension. A local three-dimensional MHD simulation representative of near-flank magnetopause conditions (plasma β’s 0.5 – 5 and convective Mach numbers M f ∼ 0.4) under northward interplanetary magnetic field (IMF) is used to validate λ MHD . Analysis shows it correlates well with hydrodynamic vortex definitions, though the level of correlation decreases with vortex evolution. Overall, vortical momentum dominates λ MHD at all times. During the linear growth phase, density gradients act to oppose vortex formation. By the highly nonlinear stage, the formation of small-scale structures leads to a rising importance of the magnetic tension. Compressibility was found to be insignificant throughout. Finally, a demonstration of this method adapted to tetrahedral spacecraft observations is performed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Unveiling the Finite Temperature Physics of Hydrogen Chains via Auxiliary Field Quantum Monte Carlo

The ability to accurately predict the finite temperature properties and phase diagrams of realistic quantum solids is central to uncovering new phases and engineering materials with novel properties ripe for device applications. Nonetheless, there remain comparatively few many-body techniques capable of elucidating the finite temperature physics of solids from first principles. Here, we take a significant step toward developing such a technique by generalizing our previous, exact fully ab initio finite temperature Auxiliary Field Quantum Monte Carlo (FT-AFQMC) method to model periodic solids and employing it to uncover the finite temperature physics of periodic hydrogen chains. Our chains’ unit cells consist of 10 hydrogen atoms modeled in a minimal basis, and we sample 5 k-points from the first Brillouin zone to arrive at a supercell consisting of 50 orbitals and 50 electrons. Based upon our calculations of these chains’ many-body energies, free energies, entropies, heat capacities, double and natural occupancies, and charge and spin correlation functions, we outline their metal–insulator and magnetic ordering as a function of both H–H bond distance and temperature. At low temperatures approaching the ground state, we observe both metal–insulator and ferromagnetic–antiferromagnetic crossovers at bond lengths between 0.5 and 0.75 Å. We then demonstrate how this low-temperature ordering evolves into a metallic phase with decreasing magnetic order at higher temperatures. In order to contextualize our results, we compare the features we observe to those previously seen in one-dimensional, half-filled Hubbard models at finite temperature and in ground state hydrogen chains. Interestingly, we identify signatures of the Pomeranchuk effect in hydrogen chains for the first time and show that spin and charge excitations that typically arise at distinct temperatures in the Hubbard model are indistinguishably coupled in these systems. Beyond qualitatively revealing the many-body phase behavior of hydrogen chains in a numerically exact manner without invoking the phaseless approximation, our efforts shed light on the further theoretical developments that will be required to construct the phase diagrams of the more complex transition metal, lanthanide, and actinide solids of longstanding interest to physicists.

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

Quantum paramagnetism in a non-Kramers rare-earth oxide: Monoclinic Pr 2 Ti 2 O 7

Little is so far known about the magnetism of the A 2 B 2 O 7 monoclinic layered perovskites that replace the spin-ice supporting pyrochlore structure for r A /r B > 1.78. We show that high quality monoclinic Pr 2 Ti 2 O 7 single crystals with a three-dimensional network of non-Kramers Pr 3+ ions that interact through edge-sharing superexchange interactions, form a singlet ground-state quantum paramagnet that does not undergo any magnetic phase transitions down to, at least, 1.8 K. The chemical phase stability, structure, and magnetic properties of the layered perovskite Pr 2 Ti 2 O 7 were investigated using x-ray diffraction, transmission electron microscopy, and magnetization measurements. Synthesis of polycrystalline samples with the nominal compositions of Pr 2 Ti 2+x O 7 (–0.16 ≤ x ≤ 0.16 ) showed that deviations from the Pr 2 Ti 2 O 7 stoichiometry lead to secondary phases of related structures including the perovskite phase Pr 2/3 TiO 3 and the orthorhombic phases Pr 4 Ti 9 O 24 and Pr 2 TiO 5 . No indications of site disordering (stuffing and antistuffing) or vacancy defects were observed in the Pr 2 Ti 2 O 7 majority phase. A procedure for growth of high-structural-quality stoichiometric single crystals of Pr 2 Ti 2 O 7 by the traveling solvent floating zone method is reported. Thermomagnetic measurements of single-crystalline Pr 2 Ti 2 O 7 reveal an isolated singlet ground state that we associate with the low-symmetry crystal electric-field environments that split the (2J + 1 = 9)-fold degenerate spin-orbital multiplets of the four differently coordinated Pr 3+ ions into 36 isolated singlets resulting in an anisotropic temperature-independent van Vleck susceptibility at low T. Here, a small isotropic Curie term is associated with 0.96(2)% noninteracting Pr 4+ impurities.

36 MATERIALS SCIENCE↗