SEARCH · Engineering Papers
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.
The mean-field theory of dilute spin-glasses is studied in the limit where the average coordination number is finite. The zero-temperature phase diagram is calculated and the relationship between the spin-glass phase and the percolation transition is discussed. The present formalism is applicable also to graph optimization problems.
Li-excess disordered rocksalts (DRXs) are emerging as promising cathode materials for Li-ion batteries due to their ability to use earth-abundant transition metals. In this work, a new strategy based on partial Li deficiency engineering is introduced to optimize the overall electrochemical performance of DRX cathodes. Specifically, by using Mn-based DRX as a proof-of-concept, it is demonstrated that the introduction of cation vacancies during synthesis (e.g., Li 1.3-x Mn 2+ 0.4-x Mn 3+ x Nb 0.3 O 1.6 F 0.4 , x = 0, 0.2, and 0.4) improves both the discharge capacity and rate performance due to the more favored short-range order in the presence of Mn3+. Density functional theory calculations and Monte Carlo simulations, in combination with spectroscopic tools, reveal that introducing 10% vacancies (Li 1.1 Mn 2+ 0.2 Mn 3+ 0.2 Nb 0.3 O 1.6 F 0.4 ) enables both Mn 2+ /Mn 3+ redox and excellent Li percolation. However, a more aggressive vacancy doping (e.g., 20% vacancies in Li 0.9 Mn 3+ 0.4 Nb 0.3 O 1.6 F 0.4 ) impairs performance because it induces phase separation between an Mn-rich and a Li-rich phase.
The quantum anomalous Hall (QAH) effect is sometimes observed in twisted bilayer graphene (tBG) when it is nearly aligned with an encapsulating hexagonal boron nitride (hBN) layer. We propose that the appearance or absence of the QAH effect in individual devices could be related to commensurability between the graphene/graphene and graphene/hBN moiré patterns. We identify a series of points in the (θ GG , θ GBN ) twist-angle space at which the two moiré patterns are commensurate, allowing moiré band theory to be applied, and we show that the band Chern numbers are in this case sensitive to a rigid in-plane hBN displacement. Given this property, we argue that the QAH effect is likely only when (i) the (θ GG , θ GBN ) twist-angle-pair is close enough to a commensurate point that the two moiré patterns yield a supermoiré pattern with a sufficiently long length scale, and (ii) the supermoiré has a percolating topologically nontrivial QAH phase. For twist angles far from commensurability, the hBN layer acts as a source of disorder that can destroy the QAH effect. Our proposal can explain a number of current experimental observations. Finally, further experimental studies that can test this proposal more directly are suggested.
Not Available
Dispersal can be critical to the maintenance of ecosystems as it allows local communities to be recolonized after extinction. However, it remains unclear whether the extinction-mitigating effect of dispersal persists when the number of competing species is large. Based on a spatially explicit mathematical description of metacommunities, we show that when many species coexist, each species operates near its extinction threshold, barely surviving due to dispersal. This has general consequences for spatiotemporal abundance patterns. For short-range dispersal, species organize into fractal spatiotemporal extinction patterns characteristic of a directed percolation phase transition. As species approach their extinction threshold, biodiversity is very sensitive to perturbation, suggesting that dispersal within a metacommunity puts tight constraints on the robustness and evolution of species-rich metacommunities. Biodiversity is often attributed to a dynamic equilibrium between the immigration and extinction of species. This equilibrium forms a common basis for studying ecosystem assembly from a static reservoir of migrants—the mainland. Yet, natural ecosystems often consist of many coupled communities (i.e., metacommunities), and migration occurs between these communities. The pool of migrants then depends on what is sustained in the ecosystem, which, in turn, depends on the dynamic migrant pool. This chicken-and-egg problem of survival and dispersal is poorly understood in communities of many competing species, except for the neutral case—the “unified neutral theory of biodiversity.” Employing spatiotemporal simulations and mean-field analyses, we show that self-consistent dispersal puts rather tight constraints on the dynamic migration–extinction equilibrium. When the number of species is large, species are pushed to the edge of their global extinction, even when competition is weak. As a consequence, the overall diversity is highly sensitive to perturbations in demographic parameters, including growth and dispersal rates. When dispersal is short range, the resulting spatiotemporal abundance patterns follow broad scale-free distributions that correspond to a directed percolation phase transition. The qualitative agreement of our results for short-range and long-range dispersal suggests that this self-organization process is a general property of species-rich metacommunities. Our study shows that self-sustaining metacommunities are highly sensitive to environmental change and provides insights into how biodiversity can be rescued and maintained.
The vaporization of a droplet, interacting with its neighbors in a non-dilute spray environment is examined as well as a vaporization scaling law established on the basis of a recently developed theory of renormalized droplet. The interacting droplet consists of a centrally located droplet and its vapor bubble which is surrounded by a cloud of droplets. The distribution of the droplets and the size of the cloud are characterized by a pair-distribution function. The vaporization of a droplet is retarded by the collective thermal quenching, the vapor concentration accumulated in the outer sphere, and by the limited percolative passages for mass, momentum and energy fluxes. The retardation is scaled by the local collective interaction parameters (group combustion number of renormalized droplet, droplet spacing, renormalization number and local ambient conditions). The numerical results of a selected case study reveal that the vaporization correction factor falls from unity monotonically as the group combustion number increases, and saturation is likely to occur when the group combustion number reaches 35 to 40 with interdroplet spacing of 7.5 diameters and an environment temperature of 500 K. The scaling law suggests that dense sprays can be classified into: (1) a diffusively dense cloud characterized by uniform thermal quenching in the cloud; (2) a stratified dense cloud characterized by a radial stratification in temperature by the differential thermal quenching of the cloud; or (3) a sharply dense cloud marked by fine structure in the quasi-droplet cloud and the corresponding variation in the correction factor due to the variation in the topological structure of the cloud characterized by a pair-distribution function of quasi-droplets.
Multivalent-ion batteries offer an alternative to Li-based technologies, with the potential for greater sustainability, improved safety, and higher energy density, primarily due to their rechargeable system featuring a passivating metal anode. Although a system based on the Ca 2+ /Ca couple is particularly attractive given the low electrochemical plating potential of Ca 2+ , the remaining challenge for a viable rechargeable Ca battery is to identify Ca cathodes with fast ion transport. In this work, a high-throughput computational pipeline is adapted to (1) discover novel Ca cathodes in a largely unexplored space of empty intercalation hosts and (2) develop material design rules for Ca-ion mobility. One candidate from the screening, W 2 O 3 (PO 4 ) 2 , is confirmed to have a low Nudged Elastic Band (NEB) barrier of 168 meV within a one-dimensional (1D) ion percolation topology. This candidate is subsequently synthesized and electrochemically tested, achieving reversible Ca cycling with a capacity of 25 mA h/g. To further accelerate the screening for promising Ca intercalation electrodes, machine learning (ML) Random Forest (RF) and Extreme Gradient Boosting (XGB) classification models are created with local environment descriptors based on a large, structurally and chemically diverse dataset of minimum energy pathways, spanning over 5,000 density functional theory (DFT) site energy calculations. Accuracies of 92% are achieved, material design metrics are quantified, ML force-fields are leveraged in an accelerated iteration of the screening, and a total of 27 novel Ca cathode materials are highlighted for further investigation.
We study the cosmological transition of 5D warped compactifications, from the high-temperature black-brane phase to the low-temperature Randall-Sundrum I phase. The transition proceeds via percolation of bubbles of IR-brane nucleating from the black-brane horizon. The violent bubble dynamics can be a powerful source of observable stochastic gravitational waves. While bubble nucleation is non-perturbative in 5D gravity, it is amenable to semiclassical treatment in terms of a “bounce” configuration interpolating between the two phases. We demonstrate how such a bounce configuration can be smooth enough to maintain 5D effective field theory control, and how a simple ansatz for it places a rigorous lower-bound on the transition rate in the thin-wall regime, and gives plausible estimates more generally. When applied to the Hierarchy Problem, the minimal Goldberger-Wise stabilization of the warped throat leads to a slow transition with significant supercooling. We demonstrate that a simple generalization of the Goldberger-Wise potential modifies the IR-brane dynamics so that the transition completes more promptly. Supercooling determines the dilution of any (dark) matter abundances generated before the transition, potentially at odds with data, while the prompter transition resolves such tensions. We discuss the impact of the different possibilities on the strength of the gravitational wave signals. Via AdS/CFT duality the warped transition gives a theoretically tractable holographic description of the 4D Composite Higgs (de)confinement transition. Our generalization of the Goldberger-Wise mechanism is dual to, and concretely models, our earlier proposal in which the composite dynamics is governed by separate UV and IR RG fixed points. The smooth 5D bounce configuration we introduce complements the 4D dilaton/radion dominance derivation presented in our earlier work.
Explore the source record for details and available documents.
Abstract Point defect formation and migration in oxides governs a wide range of phenomena from corrosion kinetics and radiation damage evolution to electronic properties. In this study, we examine the thermodynamics and kinetics of anion and cation point defects using density functional theory in hematite ( α -Fe 2 O 3 ), an important iron oxide highly relevant in both corrosion of steels and water-splitting applications. These calculations indicate that the migration barriers for point defects can vary significantly with charge state, particularly for cation interstitials. Additionally, we find multiple possible migration pathways for many of the point defects in this material, related to the low symmetry of the corundum crystal structure. The possible percolation paths are examined, using the barriers to determine the magnitude and anisotropy of long-range diffusion. Our findings suggest highly anisotropic mass transport in hematite, favoring diffusion along the c -axis of the crystal. In addition, we have considered the point defect formation energetics using the largest Fe 2 O 3 supercell reported to date.
Upon intense femtosecond photo-excitation, a many-body system can undergo a phase transition through a non-equilibrium route, but understanding these pathways remains an outstanding challenge. Here, we use time-resolved second harmonic generation to investigate a photo-induced phase transition in Ca 3 Ru 2 O 7 and show that mesoscale inhomogeneity profoundly influences the transition dynamics. We observe a marked slowing down of the characteristic time τ that quantifies the transition between two structures. τ evolves non-monotonically as a function of photo-excitation fluence, rising from below 200 fs to ~1.4 ps, then falling again to below 200 fs. To account for the observed behavior, we perform a bootstrap percolation simulation that demonstrates how local structural interactions govern the transition kinetics. Finally, our work highlights the importance of percolating mesoscale inhomogeneity in the dynamics of photo-induced phase transitions and provides a model that may be useful for understanding such transitions more broadly.
There is significant technological interest in developing ever faster switching between different electronic and magnetic states of matter. Manipulating properties at terahertz rates requires accessing the intrinsic timescales of both electrons and associated phonons, which is possible with short-pulse photoexcitation. However, in many Mott insulators, the electronic transition is accompanied by the nucleation and growth of percolating domains of the changed lattice structure, leading to empirical timescales dominated by slowly coarsening dynamics. Here, in this study, we use time-resolved X-ray diffraction and reflectivity measurements to show that the photoinduced insulator-to-metal transition in an epitaxially strained Mott insulating thin film occurs without observable domain formation and coarsening effects, allowing the study of the intrinsic electronic and lattice dynamics. Above a fluence threshold, the initial electronic excitation drives a fast lattice rearrangement, which is followed by a slower electronic evolution into a metastable nonequilibrium state. Microscopic model calculations based on time-dependent dynamical mean-field theory and semiclassical lattice dynamics explain the threshold behaviour and elucidate the delayed onset of the electronic phase transition. This work highlights the importance of combined electronic and structural studies in unravelling the physics of dynamic transitions and the timescales of photoinduced processes. During a photoinduced phase transition, electronic rearrangements are usually faster than lattice ones. Time-resolved measurements now show that the insulator-to-metal transition in a thin-film Mott insulator is preceded by lattice reconfiguration.
Infiltration processes in fractured-porous media remain a crucial, yet not very well understood component of recharge and vulnerability assessment. Under partially-saturated conditions flows in fractures, percolating fracture networks, and fault zones contribute to the fastest spectrum of infiltration velocities via preferential pathways. Specifically, the partitioning dynamics at fracture intersections determine the magnitude of flow fragmentation into vertical and horizontal components and hence the bulk flow velocity and dispersion of fracture networks. In this work, we derive an approximate analytical solution for the partitioning process and validate it using smoothed particle hydrodynamics simulations. The developed transfer function is conceptually based on simulation results and laboratory experiments carried out in previous works. It allows to efficiently simulate flow through fracture networks with simple cubic structure and arbitrary number of fractures and aperture sizes via linear response theory and convolution of a given input signal. We derive a non-dimensional bulk flow velocity ($\widetilde{v}$) and dispersion coefficient ($\widetilde{D}$) to characterize fracture networks in terms of dimensionless horizontal and vertical time scales $\tau_m$ and $\tau_0$. The dispersion coefficient is shown to strongly depend on the horizontal time scale and converges towards a constant value of $0.08$ within reasonable fluid and geometrical parameter ranges, while the non-dimensional velocity exhibits a characteristic $\widetilde{v} \sim \tau_m^{-1/2}$ scaling. Given that hydraulic information is often only available at limited places within (fractured-porous) aquifer system, such as boreholes or springs, our study intends to provide a rudimentary analytical concept to potentially reconstruct internal fracture network geometries from external boundary information, e.g., the dispersive properties of discharge (groundwater level fluctuations).
Introducing particles to dynamic covalent networks is a common approach to improve their performance. However, network formation can be impacted by their size and functionality. The influence can be predicted by common theories for small molecular precursors, but it is unclear whether they are applicable to precursors bearing numerous reactive groups and micrometer-scale dimensions. In this work, an experimental study was undertaken using dynamic covalent networks formed by the Diels–Alder reaction between furan and maleimide groups. The gelation behavior of the Diels–Alder networks was studied using rheometry to track their network formation at 40 °C with varying maleimide-functionalized microsphere loading. The highly functional microspheres can interact with the furan precursor, aiding in the formation of the Diels–Alder networks. A 5 wt% microsphere sample can reduce the gelation time by 23% and facilitate network formation in an unbalanced stoichiometry near the critical composition to form a percolating network.
Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.
The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024
Mechanical forces generated by myosin II molecular motors drive diverse cellular processes, most notably shape change, division and locomotion. These forces may be transmitted over long range through the cytoskeletal medium – a disordered, viscoelastic network of biopolymers. The resulting cell size scale force chains can in principle mediate mechanical interactions between distant actomyosin units, leading to self-organized structural order in the cell cytoskeleton. Inspired by such force transmission through elastic structures in the cytoskeleton, we consider a percolated fiber lattice network, where fibers are represented as linear elastic elements that can both bend and stretch, and the contractile activity of myosin motors is represented by force dipoles. Then, by using a variety of metrics, we show how two such contractile force dipoles interact with each other through their mutual mechanical deformations of the elastic fiber network. As a prelude to two-dipole interactions, we quantify how forces propagate through the network from a single anisotropic force dipole by analyzing clusters of nodes connected by highly strained bonds, as well as through the decay rate of strain energy with distance from a force dipole. We show that predominant fiber bending screens out force propagation, resulting in reduced and strongly network configuration-dependent dipole interactions. On the other hand, stretching-dominated networks support longer-ranged inter-dipole interactions that recapitulate the predictions of linear elasticity theory. By characterizing the differences between tensile and compressive force propagation in the fiber network, we show how inter-dipole interaction depends on the dipoles’ mutual separation and orientation. The resulting elastic interaction energy may mediate a force between multiple distant dipoles, leading to their self-organization into ordered configurations. In conclusion, this provides a potential pathway for active mechanical force-driven structural order in elastic biopolymer networks.
A comprehensive framework of characterizing complex self-assembled structures with a set of orientational order parameters is presented. It is especially relevant in the context of using anisotropic building blocks with various symmetries. Two classes of tensor order parameters are associated with polyhedral nematic and bond orientational order. For the latter, a variation of classical bond order parameters (BOPs) is introduced, which takes advantage of the symmetry of constituent particles and/or expected crystalline phases. These symmetrized BOPs can be averaged over an entire system or assigned locally to an individual bond. By combining that with the bond percolation procedure, one is able to identify coherent domains within a self-assembled structure. As a demonstration of the proposed framework, we apply it to a simulated hybrid system that combines isotropic and patchy particles with octahedral symmetry. In conclusion, not only does the methodology allow one to identify individual crystalline domains but also it detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.