Three-dimensional construction of hyperuniform, nonhyperuniform, and antihyperuniform disordered heterogeneous materials and their transport properties via spectral density functions

Rigorous theories connecting physical properties of a heterogeneous material to its microstructure offer a promising avenue to guide the computational material design and optimization. The spectral density function χ[over ̃]_{_{V}}(k), which can be obtained experimentally from scattering data, enables accurate determination of various transport and wave propagation characteristics, including the time-dependent diffusion spreadability S(t) and effective dynamic dielectric constant ε_{e} for electromagnetic wave propagation. Moreover, χ[over ̃]_{_{V}}(k) determines rigorous upper bounds on the fluid permeability K. Given the importance of χ[over ̃]_{_{V}}(k), we present here an efficient Fourier-space based computational framework to construct three-dimensional (3D) statistically isotropic two-phase heterogeneous materials corresponding to targeted spectral density functions. In particular, we employ a variety of analytical functional forms for χ[over ̃]_{_{V}}(k) that satisfy all known necessary conditions to construct disordered stealthy hyperuniform, standard hyperuniform, nonhyperuniform, and antihyperuniform two-phase heterogeneous material systems at varying phase volume fractions. We show that by tuning the correlations in the system across length scales via the targeted functions, one can generate a rich spectrum of distinct structures within each of the above classes of materials. Importantly, we present the first realization of antihyperuniform two-phase heterogeneous materials in 3D, which are characterized by autocovariance function χ_{_{V}}(r) with a power-law tail, resulting in microstructures that contain clusters of dramatically different sizes and morphologies. We also determine the diffusion spreadability S(t) and estimate the fluid permeability K associated with all of the constructed materials directly from the corresponding spectral densities. Although it is well established that the long-time asymptotic scaling behavior of S(t) only depends on the functional form of χ[over ̃]_{_{V}}(k), with the stealthy hyperuniform and antihyperuniform media, respectively, achieving the most and least efficient transport, we find that varying the length-scale parameter within each class of χ[over ̃]_{_{V}}(k) functions can also lead to orders of magnitude variation of S(t) at intermediate and long time scales. Moreover, we find that increasing the solid volume fraction ϕ_{1} and correlation length a in the constructed media generally leads to a decrease in the dimensionless fluid permeability K/a^{2}, while the antihyperuniform media possess the largest K/a^{2} among the four classes of materials with the same ϕ_{1} and a. These results indicate the feasibility of employing parameterized spectral densities for designing composites with targeted transport properties.

Anomalous suppression of large-scale density fluctuations in classical and quantum spin liquids

Classical spin liquids (CSLs) are intriguing states of matter that do not exhibit long-range magnetic order and are characterized by an extensive ground-state degeneracy. Adding quantum fluctuations, which induce dynamics between these different classical ground states, can give rise to quantum spin liquids (QSLs). QSLs are highly entangled quantum phases of matter characterized by fascinating emergent properties, such as fractionalized excitations and topological order. One such exotic quantum liquid is the [Formula: see text] QSL, which can be regarded as a resonating valence bond (RVB) state formed from superpositions of dimer coverings of an underlying lattice. In this work, we unveil a hidden large-scale structural property of archetypal CSLs and QSLs known as hyperuniformity, i.e., normalized infinite-wavelength density fluctuations are completely suppressed in these systems. In particular, we first demonstrate that classical ensembles of close-packed dimers and their corresponding quantum RVB states are perfectly hyperuniform in general. Subsequently, we focus on a ruby-lattice spin liquid that was recently realized in a Rydberg-atom quantum simulator, and show that the QSL remains effectively hyperuniform even in the presence of a finite density of spinon and vison excitations, as long as the dimer constraint is still largely preserved. Moreover, we demonstrate that metrics based on the framework of hyperuniformity can be used to distinguish the QSL from other proximate quantum phases. These metrics can help identify potential QSL candidates, which can then be further analyzed using more advanced, computationally intensive quantum numerics to confirm their status as true QSLs.

Spectroscopy of two-dimensional interacting lattice electrons using symmetry-aware neural backflow transformations

Neural networks have shown to be a powerful tool to represent the ground state of quantum many-body systems, including fermionic systems. However, efficiently integrating lattice symmetries into neural representations remains a significant challenge. In this work, we introduce a framework for embedding lattice symmetries in fermionic wavefunctions and demonstrate its ability to target both ground states and low-lying excitations. Using group-equivariant neural backflow transformations, we study the t-V model on a square lattice away from half-filling. Our symmetry-aware backflow significantly improves ground-state energies and yields accurate low-energy excitations for lattices up to 10 × 10. We also compute accurate two-point density-correlation functions and the structure factor to identify phase transitions and critical points. These findings introduce a symmetry-aware framework important for studying quantum materials and phase transitions.

Highly accurate real-space electron densities with neural networks

Variational ab initio methods in quantum chemistry stand out among other methods in providing direct access to the wave function. This allows, in principle, straightforward extraction of any other observable of interest, besides the energy, but, in practice, this extraction is often technically difficult and computationally impractical. Here, we consider the electron density as a central observable in quantum chemistry and introduce a novel method to obtain accurate densities from real-space many-electron wave functions by representing the density with a neural network that captures known asymptotic properties and is trained from the wave function by score matching and noise-contrastive estimation. We use variational quantum Monte Carlo with deep-learning Ansätze to obtain highly accurate wave functions free of basis set errors and from them, using our novel method, correspondingly accurate electron densities, which we demonstrate by calculating dipole moments, nuclear forces, contact densities, and other density-based properties.

Unified Variational Approach Description of Ground-State Phases of the Two-Dimensional Electron Gas

The two-dimensional electron gas (2DEG) is a fundamental model, which is drawing increasing interest because of recent advances in experimental and theoretical studies of 2D materials. Current understanding of the ground state of the 2DEG relies on quantum Monte Carlo calculations, based on variational comparisons of different Ansätze for different phases. We use a single variational ansatz, a general backflow-type wave function using a message-passing neural quantum state architecture, for a unified description across the entire density range. The variational optimization consistently leads to lower ground-state energies than previous best results. Transition into a Wigner crystal (WC) phase occurs automatically at r_{s}=37±1, a density lower than currently believed. Between the liquid and WC phases, the same ansatz and variational search strongly suggest the existence of intermediate states in a broad range of densities, with enhanced short-range nematic spin correlations.

Spin-symmetry-enforced solution of the many-body Schrödinger equation with a deep neural network

The integration of deep neural networks with the variational Monte Carlo (VMC) method has marked a substantial advancement in solving the Schrödinger equation. In this work we enforce spin symmetry in the neural-network-based VMC calculation using a modified optimization target. Our method is designed to solve for the ground state and multiple excited states with target spin symmetry at a low computational cost. It predicts accurate energies while maintaining the correct symmetry in strongly correlated systems, even in cases in which different spin states are nearly degenerate. Our approach also excels at spin-gap calculations, including the singlet-triplet gap in biradical systems, which is of high interest in photochemistry. Overall, this work establishes a robust framework for efficiently calculating various quantum states with specific spin symmetry in correlated systems.

Formally Exact and Practically Useful Analytic Solution of Harmonium

We provide a novel exact analytic solution of harmonium with arbitrary Coulomb interaction strength, for ground as well as all the excited states, using our recently developed method for solving Schrödinger equations. By comparing three formally exact analytic representations of the wave function including the one that utilizes biconfluent Heun function, we find that the best and practically useful representation for the ground state is given by an exact factorized form involving a noninteger power pre-exponential factor, an exponentially decaying term and a modulator function. For excited states, additional factors are needed to account for the nodal information. We show that our method is far more efficient than basis-expansion-based methods in representing the wave function. With the exact wave functions, we have also analyzed the evolution trends of the electron density and natural occupation numbers with increasing interaction strength, which gives insight into the interesting physics in the strong correlation limit.

Electronic rotons and Wigner crystallites in a two-dimensional dipole liquid

A key concept proposed by Landau to explain superfluid liquid helium is the elementary excitation of quantum particles called rotons1-8. The irregular arrangement of atoms in a liquid leads to the aperiodic dispersion of rotons, which played a pivotal role in understanding fractional quantum Hall liquids (magneto-rotons)9,10 and the supersolidity of Bose-Einstein condensates11-13. Even for a two-dimensional electron or dipole liquid, in the absence of a magnetic field, the repulsive interactions have been predicted to form a roton minimum14-19, which can be used to trace the transition to Wigner crystals20-24 and superconductivity25-27, although this has not yet been observed. Here, we report the observation of such electronic rotons in a two-dimensional dipole liquid of alkali-metal ions donating electrons to surface layers of black phosphorus. Our data reveal the striking aperiodic dispersion of rotons, which is characterized by a local minimum of energy at finite momentum. As the density of dipoles decreases so that interactions dominate over the kinetic energy, the roton gap reduces to 0, as in a crystal, signalling Wigner crystallization. Our model shows the importance of short-range order arising from repulsion between dipoles, which can be viewed as the formation of Wigner crystallites (bubbles or stripes) floating in the sea of a Fermi liquid. Our results reveal that the primary origin of electronic rotons (and the pseudogap) is strong correlations.

Local-to-non-local transition laws for stiffness-tuneable monoatomic chains preserving springs mass

Recently, non-local configurations have been proposed by adding beyond nearest neighbour couplings among elements in lattices to obtain roton-like dispersion relations and phase and group velocities with opposite signs. Even though the introduction of non-local elastic links in metamaterials has unlocked unprecedented possibilities, literature models and prototypes seem neither to provide criteria to compare local and non-local lattices nor to discuss any related rules governing the transition between the two configurations. A physically reasonable principle that monoatomic one-dimensional chains must obey to pass from single- to multi-connected systems is here proposed through a mass conservation law for elastic springs thereby introducing a suitable real dimensionless parameter [Formula: see text] to tune stiffness distribution. Therefore, the dispersion relations as a function of [Formula: see text] and of the degree of non-locality [Formula: see text] are derived analytically, demonstrating that the proposed principle can be rather interpreted as a general mechanical consistency condition to preserve proper dynamics, involving the spring-to-bead mass ratio. Finally, after discussing qualitative results and deriving some useful inequalities, numerical simulations and two-dimensional FFTs are performed for some paradigmatic examples to highlight key dynamics features exhibited by chains with finite length as the parameters [Formula: see text] and [Formula: see text] vary.This article is part of the theme issue 'Current developments in elastic and acoustic metamaterials science (Part 2)'.

Universal Hyperuniform Organization in Looped Leaf Vein Networks

The leaf vein network is a hierarchical vascular system that transports water and nutrients to the leaf cells. The thick primary veins form a branched network, while the secondary veins can develop closed loops forming a well-defined cellular structure. Through extensive analysis of a variety of distinct leaf species, we discover that the apparently disordered cellular structures of the secondary vein networks exhibit a universal hyperuniform organization and possess a hidden order on large scales. Disorder hyperuniform systems lack conventional long-range order, yet they completely suppress normalized infinite-wavelength density fluctuations like crystals. Specifically, we find that the distributions of the geometric centers associated with the vein network loops possess a vanishing static structure factor in the limit that the wave number k goes to 0, i.e., S(k)∼k^{α}, where α≈0.64±0.021, providing an example of class III hyperuniformity in biology. This hyperuniform organization leads to superior efficiency of diffusive transport, as evidenced by the much faster convergence of the time-dependent spreadability S(t) to its longtime asymptotic limit, compared to that of other uncorrelated or correlated disordered but nonhyperuniform organizations. Our results also have implications for the discovery and design of novel disordered network materials with optimal transport properties.

Hyperuniformity classes of quasiperiodic tilings via diffusion spreadability

Hyperuniform point patterns can be classified by the hyperuniformity scaling exponent α>0, that characterizes the power-law scaling behavior of the structure factor S(k) as a function of wave number k≡|k| in the vicinity of the origin, e.g., S(k)∼|k|^{α} in cases where S(k) varies continuously with k as k→0. In this paper, we show that the spreadability is an effective method for determining α for quasiperiodic systems where S(k) is discontinuous and consists of a dense set of Bragg peaks. It has been shown in [Phys. Rev. E 104, 054102 (2021)10.1103/PhysRevE.104.054102] that, for media with finite α, the long-time behavior of the excess spreadability S(∞)-S(t) can be fit to a power law of the form t^{-(d-α)/2}, where d is the space dimension, to accurately extract α for the continuous case. We first transform quasiperiodic and limit-periodic point patterns into two-phase media by mapping them onto packings of identical nonoverlapping disks, where space interior to the disks represents one phase and the space in exterior to them represents the second phase. We then compute the spectral density χ[over ̃]_{_{V}}(k) of the packings, and finally compute and fit the long-time behavior of their excess spreadabilities. Specifically, we show that the excess spreadability can be used to accurately extract α for the one-dimensional (1D) limit-periodic period-doubling chain (α=1) and the 1D quasicrystalline Fibonacci chain (α=3) to within 0.02% of the analytically known exact results. Moreover, we obtain a value of α=5.97±0.06 for the two-dimensional Penrose tiling and present plausible theoretical arguments strongly suggesting that α is exactly equal to six. We also show that, due to the self-similarity of the structures examined here, one can truncate the small-k region of the scattering information used to compute the spreadability and obtain an accurate value of α, with a small deviation from the untruncated case that decreases as the system size increases. This strongly suggests that one can obtain a good estimate of α for an infinite self-similar quasicrystal from a modestly sized finite sample. The methods described here offer a simple and general procedure to characterize accurately the large-scale translational order present in quasicrystalline and limit-periodic media in any space dimension that are self-similar. Moreover, the scattering information extracted from these two-phase media encoded in χ[over ̃]_{_{V}}(k), can be used to estimate their physical properties, such as their effective dynamic dielectric constants, effective dynamic elastic constants, and fluid permeabilities.

Vibrational properties of disordered stealthy hyperuniform 1D atomic chains

Disorder hyperuniformity is a recently discovered exotic state of many-body systems that possess a hidden order in between that of a perfect crystal and a completely disordered system. Recently, this novel disordered state has been observed in a number of quantum materials including amorphous 2D graphene and silica, which are endowed with unexpected electronic transport properties. Here, we numerically investigate 1D atomic chain models, including perfect crystalline, disordered stealthy hyperuniform (SHU) as well as randomly perturbed atom packing configurations to obtain a quantitative understanding of how the unique SHU disorder affects the vibrational properties of these low-dimensional materials. We find that the disordered SHU chains possess lower cohesive energies compared to the randomly perturbed chains, implying their potential reliability in experiments. Our inverse partition ratio (IPR) calculations indicate that the SHU chains can support fully delocalized states just like perfect crystalline chains over a wide range of frequencies, i.e.ω∈(0,100)cm-1, suggesting superior phonon transport behaviors within these frequencies, which was traditionally considered impossible in disordered systems. Interestingly, we observe the emergence of a group of highly localized states associated withω∼200cm-1, which is characterized by a significant peak in the IPR and a peak in phonon density of states at the corresponding frequency, and is potentially useful for decoupling electron and phonon degrees of freedom. These unique properties of disordered SHU chains have implications in the design and engineering of novel quantum materials for thermal and phononic applications.

Field control of quasiparticle decay in a quantum antiferromagnet

Dynamics in a quantum material is described by quantized collective motion: a quasiparticle. The single-quasiparticle description is useful for a basic understanding of the system, whereas a phenomenon beyond the simple description such as quasiparticle decay which affects the current carried by the quasiparticle is an intriguing topic. The instability of the quasiparticle is phenomenologically determined by the magnitude of the repulsive interaction between a single quasiparticle and the two-quasiparticle continuum. Although the phenomenon has been studied in several materials, thermodynamic tuning of the quasiparticle decay in a single material has not yet been investigated. Here we show, by using neutron scattering, magnetic field control of the magnon decay in a quantum antiferromagnet RbFeCl3, where the interaction between the magnon and continuum is tuned by the field. At low fields where the interaction is small, the single magnon decay process is observed. In contrast, at high fields where the interaction exceeds a critical magnitude, the magnon is pushed downwards in energy and its lifetime increases. Our study demonstrates that field control of quasiparticle decay is possible in the system where the two-quasiparticle continuum covers wide momentum-energy space, and the phenomenon of the magnon avoiding decay is ubiquitous.

Ab initio quantum chemistry with neural-network wavefunctions

Deep learning methods outperform human capabilities in pattern recognition and data processing problems and now have an increasingly important role in scientific discovery. A key application of machine learning in molecular science is to learn potential energy surfaces or force fields from ab initio solutions of the electronic Schrödinger equation using data sets obtained with density functional theory, coupled cluster or other quantum chemistry (QC) methods. In this Review, we discuss a complementary approach using machine learning to aid the direct solution of QC problems from first principles. Specifically, we focus on quantum Monte Carlo methods that use neural-network ansatzes to solve the electronic Schrödinger equation, in first and second quantization, computing ground and excited states and generalizing over multiple nuclear configurations. Although still at their infancy, these methods can already generate virtually exact solutions of the electronic Schrödinger equation for small systems and rival advanced conventional QC methods for systems with up to a few dozen electrons.

Computational design of anisotropic stealthy hyperuniform composites with engineered directional scattering properties

Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties, including large isotropic photonic band gaps, superior resistance to fracture, and nearly optimal electrical and thermal transport properties, to name but a few. Here we generalize the Fourier-space-based numerical construction procedure for designing and generating digital realizations of isotropic disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. 142, 152 (2018)1359-645410.1016/j.actamat.2017.09.053] to anisotropic microstructures with targeted spectral densities. Our generalized construction procedure explicitly incorporates the vector-dependent spectral density function χ[over ̃]_{_{V}}(k) of arbitrary form that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of anisotropic stealthy hyperuniform microstructures with χ[over ̃]_{_{V}}(k)=0 for k∈Ω, i.e., complete suppression of scattering in an "exclusion" region Ω around the origin in Fourier space. We show how different exclusion-region shapes with various discrete symmetries, including circular-disk, elliptical-disk, square, rectangular, butterfly-shaped, and lemniscate-shaped regions of varying size, affect the resulting statistically anisotropic microstructures as a function of the phase volume fraction. The latter two cases of Ω lead to directionally hyperuniform composites, which are stealthy hyperuniform only along certain directions and are nonhyperuniform along others. We find that while the circular-disk exclusion regions give rise to isotropic hyperuniform composite microstructures, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotropic exclusion regions impose strong constraints on the global symmetry of the resulting media, they can still possess structures at a local level that are nearly isotropic. Both the isotropic and anisotropic hyperuniform microstructures associated with the elliptical-disk, square, and rectangular Ω possess phase-inversion symmetry over certain range of volume fractions and a percolation threshold ϕ_{c}≈0.5. On the other hand, the directionally hyperuniform microstructures associated with the butterfly-shaped and lemniscate-shaped Ω do not possess phase-inversion symmetry and percolate along certain directions at much lower volume fractions. We also apply our general procedure to construct stealthy nonhyperuniform systems. Our construction algorithm enables one to control the statistical anisotropy of composite microstructures via the shape, size, and symmetries of Ω, which is crucial to engineering directional optical, transport, and mechanical properties of two-phase composite media.

Extreme Spatial Dispersion in Nonlocally Resonant Elastic Metamaterials

To date, the vast majority of architected materials have leveraged two physical principles to control wave behavior, namely, Bragg interference and local resonances. Here, we describe a third path: structures that accommodate a finite number of delocalized zero-energy modes, leading to anomalous dispersion cones that nucleate from extreme spatial dispersion at 0 Hz. We explain how to design such zero-energy modes in the context of elasticity and show that many of the landmark wave properties of metamaterials can also be induced at an extremely subwavelength scale by the associated anomalous cones, without suffering from the same bandwidth limitations. We then validate our theory through a combination of simulations and experiments. Finally, we present an inverse design method to produce anomalous cones at desired locations in k space.

Precise determination of pair interactions from pair statistics of many-body systems in and out of equilibrium

The determination of the pair potential v(r) that accurately yields an equilibrium state at positive temperature T with a prescribed pair correlation function g_{2}(r) or corresponding structure factor S(k) in d-dimensional Euclidean space R^{d} is an outstanding inverse statistical mechanics problem with far-reaching implications. Recently, Zhang and Torquato [Phys. Rev. E 101, 032124 (2020)2470-004510.1103/PhysRevE.101.032124] conjectured that any realizable g_{2}(r) or S(k) corresponding to a translationally invariant nonequilibrium system can be attained by a classical equilibrium ensemble involving only (up to) effective pair interactions. Testing this conjecture for nonequilibrium systems as well as for nontrivial equilibrium states requires improved inverse methodologies. We have devised an optimization algorithm to precisely determine effective pair potentials that correspond to pair statistics of general translationally invariant disordered many-body equilibrium or nonequilibrium systems at positive temperatures. This methodology utilizes a parameterized family of pointwise basis functions for the potential function whose initial form is informed by small-, intermediate- and large-distance behaviors dictated by statistical-mechanical theory. Subsequently, a nonlinear optimization technique is utilized to minimize an objective function that incorporates both the target pair correlation function g_{2}(r) and structure factor S(k) so that the small intermediate- and large-distance correlations are very accurately captured. To illustrate the versatility and power of our methodology, we accurately determine the effective pair interactions of the following four diverse target systems: (1) Lennard-Jones system in the vicinity of its critical point, (2) liquid under the Dzugutov potential, (3) nonequilibrium random sequential addition packing, and (4) a nonequilibrium hyperuniform "cloaked" uniformly randomized lattice. We found that the optimized pair potentials generate corresponding pair statistics that accurately match their corresponding targets with total L_{2}-norm errors that are an order of magnitude smaller than that of previous methods. The results of our investigation lend further support to the Zhang-Torquato conjecture. Furthermore, our algorithm will enable one to probe systems with identical pair statistics but different higher-body statistics, which will shed light on the well-known degeneracy problem of statistical mechanics.

Experimental observation of roton-like dispersion relations in metamaterials

Previously, rotons were observed in correlated quantum systems at low temperatures, including superfluid helium and Bose-Einstein condensates. Here, following a recent theoretical proposal, we report the direct experimental observation of roton-like dispersion relations in two different three-dimensional metamaterials under ambient conditions. One experiment uses transverse elastic waves in microscale metamaterials at ultrasound frequencies. The other experiment uses longitudinal air-pressure waves in macroscopic channel–based metamaterials at audible frequencies. In both experiments, we identify the roton-like minimum in the dispersion relation that is associated to a triplet of waves at a given frequency. Our work shows that designed interactions in metamaterials beyond the nearest neighbors open unprecedented experimental opportunities to tailor the lowest dispersion branch—while most previous metamaterial studies have concentrated on shaping higher dispersion branches.

Roton pair density wave in a strong-coupling kagome superconductor

The transition metal kagome lattice materials host frustrated, correlated and topological quantum states of matter^1-9. Recently, a new family of vanadium-based kagome metals, AV₃Sb₅ (A = K, Rb or Cs), with topological band structures has been discovered^10,11. These layered compounds are nonmagnetic and undergo charge density wave transitions before developing superconductivity at low temperatures^11-19. Here we report the observation of unconventional superconductivity and a pair density wave (PDW) in CsV₃Sb₅ using scanning tunnelling microscope/spectroscopy and Josephson scanning tunnelling spectroscopy. We find that CsV₃Sb₅ exhibits a V-shaped pairing gap Δ ~ 0.5 meV and is a strong-coupling superconductor (2Δ/k_BT_c ~ 5) that coexists with 4a₀ unidirectional and 2a₀ × 2a₀ charge order. Remarkably, we discover a 3Q PDW accompanied by bidirectional 4a₀/3 spatial modulations of the superconducting gap, coherence peak and gap depth in the tunnelling conductance. We term this novel quantum state a roton PDW associated with an underlying vortex-antivortex lattice that can account for the observed conductance modulations. Probing the electronic states in the vortex halo in an applied magnetic field, in strong field that suppresses superconductivity and in zero field above T_c, reveals that the PDW is a primary state responsible for an emergent pseudogap and intertwined electronic order. Our findings show striking analogies and distinctions to the phenomenology of high-T_c cuprate superconductors, and provide groundwork for understanding the microscopic origin of correlated electronic states and superconductivity in vanadium-based kagome metals.

Roton-like acoustical dispersion relations in 3D metamaterials

Roton dispersion relations have been restricted to correlated quantum systems at low temperatures, such as liquid Helium-4, thin films of Helium-3, and Bose–Einstein condensates. This unusual kind of dispersion relation provides broadband acoustical backward waves, connected to energy flow vortices due to a “return flow”, in the words of Feynman, and three different coexisting acoustical modes with the same polarization at one frequency. By building mechanisms into the unit cells of artificial materials, metamaterials allow for molding the flow of waves. So far, researchers have exploited mechanisms based on various types of local resonances, Bragg resonances, spatial and temporal symmetry breaking, topology, and nonlinearities. Here, we introduce beyond-nearest-neighbor interactions as a mechanism in elastic and airborne acoustical metamaterials. For a third-nearest-neighbor interaction that is sufficiently strong compared to the nearest-neighbor interaction, this mechanism allows us to engineer roton-like acoustical dispersion relations under ambient conditions.

Here, the authors introduce beyond-nearest-neighbour interactions as a mechanism for molding the flow of waves in acoustic metamaterials. They find that for strong third-nearest-neighbour interactions, this mechanism allows for engineering roton-like acoustical dispersion relations under ambient conditions.

Structural characterization of many-particle systems on approach to hyperuniform states

The study of hyperuniform states of matter is an emerging multidisciplinary field, impinging on topics in the physical sciences, mathematics, and biology. The focus of this work is the exploration of quantitative descriptors that herald when a many-particle system in d-dimensional Euclidean space R^{d} approaches a hyperuniform state as a function of the relevant control parameter. We establish quantitative criteria to ascertain the extent of hyperuniform and nonhyperuniform distance-scaling regimes as well as the crossover point between them in terms of the "volume" coefficient A and "surface-area" coefficient B associated with the local number variance σ^{2}(R) for a spherical window of radius R. The larger the ratio B/A, the larger the hyperuniform scaling regime, which becomes of infinite extent in the limit B/A→∞. To complement the known direct-space representation of the coefficient B in terms of the total correlation function h(r), we derive its corresponding Fourier representation in terms of the structure factor S(k), which is especially useful when scattering information is available experimentally or theoretically. We also demonstrate that the free-volume theory of the pressure of equilibrium packings of identical hard spheres that approach a strictly jammed state either along the stable crystal or metastable disordered branch dictates that such end states be exactly hyperuniform. Using the ratio B/A, as well as other diagnostic measures of hyperuniformity, including the hyperuniformity index H and the direct-correlation function length scale ξ_{c}, we study three different exactly solvable models as a function of the relevant control parameter, either density or temperature, with end states that are perfectly hyperuniform. Specifically, we analyze equilibrium systems of hard rods and "sticky" hard-sphere systems in arbitrary space dimension d as a function of density. We also examine low-temperature excited states of many-particle systems interacting with "stealthy" long-ranged pair interactions as the temperature tends to zero, where the ground states are disordered, hyperuniform, and infinitely degenerate. We demonstrate that our various diagnostic hyperuniformity measures are positively correlated with one another. The same diagnostic measures can be used to detect the degree to which imperfections in nearly hyperuniform systems cause deviations from perfect hyperuniformity. Moreover, the capacity to identify hyperuniform scaling regimes should be particularly useful in analyzing experimentally or computationally generated samples that are necessarily of finite size.

Stone-Wales defects preserve hyperuniformity in amorphous two-dimensional networks

Disordered hyperuniformity (DHU) is a recently discovered novel state of many-body systems that possesses vanishing normalized infinite-wavelength density fluctuations similar to a perfect crystal and an amorphous structure like a liquid or glass. Here, we discover a hyperuniformity-preserving topological transformation in two-dimensional (2D) network structures that involves continuous introduction of Stone-Wales (SW) defects. Specifically, the static structure factor [Formula: see text] of the resulting defected networks possesses the scaling [Formula: see text] for small wave number k, where [Formula: see text] monotonically decreases as the SW defect concentration p increases, reaches [Formula: see text] at [Formula: see text], and remains almost flat beyond this p. Our findings have important implications for amorphous 2D materials since the SW defects are well known to capture the salient feature of disorder in these materials. Verified by recently synthesized single-layer amorphous graphene, our network models reveal unique electronic transport mechanisms and mechanical behaviors associated with distinct classes of disorder in 2D materials.

Quantum liquids

This article is dedicated to Roger A Cowley and his seminal contributions to our understanding of quantum liquids, both liquid4He and3He. Roger Cowley's neutron scattering measurements of the collective and independent particle response of liquid4He were made at Chalk River Laboratories in 1965-74 chiefly with A D B (Dave) Woods. They measured the phonon-roton (P-R) mode energy, intensity and width with new precision. Particularly, they extended the measurements to higher wave vector and identified both collective and single particle response regimes. They showed that the P-R mode terminated at a finite energy as predicted by Pitaeskii rather than continuing as predicted by Feynman and Feynman and Cohen. They determined both the single P-R mode and multimode contributions to the dynamics. They made direct comparison with theory which Roger understood well. They observed the Bose-Einstein condensate (BEC) fraction in liquid4He for the first time. This appears to be the first ever observation of BEC in any Bose gas or liquid. Roger Cowley's pioneering measurements of the density excitations of liquid3He were made at the Institut Laue Langevin (ILL) in the period 1973-80. Roger, Reinhard Scherm, W G (Bill) Stirling and collaborators showed for the first time that the density response of this highly neutron absorbing liquid could indeed be observed with neutrons. They documented with others the dynamic response as a function of temperature and pressure stimulating extensive theoretical and experimental interest that continues today.

Fermionic neural-network states for ab-initio electronic structure

Neural-network quantum states have been successfully used to study a variety of lattice and continuous-space problems. Despite a great deal of general methodological developments, representing fermionic matter is however still early research activity. Here we present an extension of neural-network quantum states to model interacting fermionic problems. Borrowing techniques from quantum simulation, we directly map fermionic degrees of freedom to spin ones, and then use neural-network quantum states to perform electronic structure calculations. For several diatomic molecules in a minimal basis set, we benchmark our approach against widely used coupled cluster methods, as well as many-body variational states. On some test molecules, we systematically improve upon coupled cluster methods and Jastrow wave functions, reaching chemical accuracy or better. Finally, we discuss routes for future developments and improvements of the methods presented.

Cloaking the underlying long-range order of randomly perturbed lattices

Random, uncorrelated displacements of particles on a lattice preserve the hyperuniformity of the original lattice, that is, normalized density fluctuations vanish in the limit of infinite wavelengths. In addition to a diffuse contribution, the scattering intensity from the the resulting point pattern typically inherits the Bragg peaks (long-range order) of the original lattice. Here we demonstrate how these Bragg peaks can be hidden in the effective diffraction pattern of independent and identically distributed perturbations. All Bragg peaks vanish if and only if the sum of all probability densities of the positions of the shifted lattice points is a constant at all positions. The underlying long-range order is then "cloaked" in the sense that it cannot be reconstructed from the pair correlation function alone. On the one hand, density fluctuations increase monotonically with the strength of perturbations a, as measured by the hyperuniformity order metric Λ[over ¯]. On the other hand, the disappearance and reemergence of long-range order, depending on whether the system is cloaked as the perturbation strength increases, is manifestly captured by the τ order metric. Therefore, while the perturbation strength a may seem to be a natural choice for an order metric of perturbed lattices, the τ order metric is a superior choice. It is noteworthy that cloaked perturbed lattices allow one to easily simulate very large samples (with at least 10^{6} particles) of disordered hyperuniform point patterns without Bragg peaks.

Vortex-loop calculation of the specific heat of superfluid ^{4}He under pressure

Vortex-loop renormalization is used to compute the specific heat of superfluid ^{4}He near the lambda point at various pressures up to 26 bars. The input parameters are the pressure dependence of T_{λ} and the superfluid density, which determine the nonuniversal parameters of the vortex core energy and core size. The results for the specific heat are found to be in good agreement with experimental data, matching the expected universal pressure dependence to within about 5%. The nonuniversal critical amplitude of the specific heat is found to be in reasonable agreement, a factor of four larger than the experiments. We point out problems with recent Gross-Pitaevskii simulations that claimed the vortex-loop percolation temperature did not match the critical temperature of the superfluid phase transition.

Sound speed in Yukawa one-component plasmas across coupling regimes

A many-body system of charged particles interacting via a pairwise Yukawa potential, the so-called Yukawa one-component plasma (YOCP), is a good approximation for a variety of physical systems. Such systems are completely characterized by two parameters: the screening parameter, κ, and the nominal coupling strength, Γ. It is well known that the collective spectrum of the YOCP is governed by a longitudinal acoustic mode, both in the weakly and strongly coupled regimes. In the long-wavelength limit, the linear term in the dispersion (i.e., ω=sk) defines the sound speed s. We study the evolution of this latter quantity from the weak- through the strong-coupling regimes by analyzing the dynamic structure function S(k,ω) in the low-frequency domain. Depending on the values of Γ and κ and w=s/v_{th} (i.e., the ratio between the phase velocity of the wave and the thermal speed of the particles), we identify five domains in the (κ,Γ) parameter space in which the physical behavior of the YOCP exhibits different features. The competing physical processes are the collective Coulomb-like versus binary-collision-dominated behavior and the individual particle motion versus quasilocalization. Our principal tool of investigation is molecular dynamics (MD) computer simulation from which we obtain S(k,ω). Recent improvements in the simulation technique have allowed us to obtain a large body of high-quality data in the range Γ={0.1-10000} and κ={0.5-5}. The theoretical results based on various models are compared in order to see which one provides the most cogent physical description and the best agreement with MD data in the different domains.

Disordered hyperuniform obstacles enhance sorting of dynamically chiral microswimmers

Disordered hyperuniformity, a brand new type of arrangement with novel physical properties, provides various practical applications in extensive fields. To highlight the great potential of applying disordered hyperuniformity to active systems, a practical example is reported here by an optimal sorting of dynamically chiral microswimmers in disordered hyperuniform obstacle environments in comparison with regular or disordered ones. This optimal chirality sorting stems from a competition between advantageous microswimmer-obstacle collisions and disadvantageous trapping of microswimmers by obstacles. Based on this mechanism, optimal chirality sorting is also realized by tuning other parameters including the number density of obstacles, the strength of driven force and the noise intensity. Our findings may open a new perspective on both theoretical and experimental investigations for further applications of disordered hyperuniformity in active systems.

Particlelike Behavior of Topological Defects in Linear Wave Packets in Photonic Graphene

Topological defects, such as quantum vortices, determine the properties of quantum fluids. Their study has been at the center of activity in solid state and BEC communities. In parallel, the nontrivial behavior of linear wave packets with complex phase patterns was investigated by singular optics. Here, we study the formation, evolution, and interaction of optical vortices in wave packets at the Dirac point in photonic graphene. We show that while their exact behavior goes beyond the Dirac equation and requires a full account of the lattice properties, it can be still approximately described by an effective theory considering the phase singularities as "particles". These particles are capable of mutual interaction, with their trajectory obeying the laws of dynamics.

Backflow Transformations via Neural Networks for Quantum Many-Body Wave Functions

Obtaining an accurate ground state wave function is one of the great challenges in the quantum many-body problem. In this Letter, we propose a new class of wave functions, neural network backflow (NNB). The backflow approach, pioneered originally by Feynman and Cohen [Phys. Rev. 102, 1189 (1956)10.1103/PhysRev.102.1189], adds correlation to a mean-field ground state by transforming the single-particle orbitals in a configuration-dependent way. NNB uses a feed-forward neural network to learn the optimal transformation via variational Monte Carlo calculations. NNB directly dresses a mean-field state, can be systematically improved, and directly alters the sign structure of the wave function. It generalizes the standard backflow [L. F. Tocchio et al., Phys. Rev. B 78, 041101(R) (2008)10.1103/PhysRevB.78.041101], which we show how to explicitly represent as a NNB. We benchmark the NNB on Hubbard models at intermediate doping, finding that it significantly decreases the relative error, restores the symmetry of both observables and single-particle orbitals, and decreases the double-occupancy density. Finally, we illustrate interesting patterns in the weights and bias of the optimized neural network.

Methodology to construct large realizations of perfectly hyperuniform disordered packings

Disordered hyperuniform packings (or dispersions) are unusual amorphous two-phase materials that are endowed with exotic physical properties. Such hyperuniform systems are characterized by an anomalous suppression of volume-fraction fluctuations at infinitely long-wavelengths, compared to ordinary disordered materials. While there has been growing interest in such singular states of amorphous matter, a major obstacle has been an inability to produce large samples that are perfectly hyperuniform due to practical limitations of conventional numerical and experimental methods. To overcome these limitations, we introduce a general theoretical methodology to construct perfectly hyperuniform packings in d-dimensional Euclidean space R^{d}. Specifically, beginning with an initial general tessellation of space by disjoint cells that meets a "bounded-cell" condition, hard particles of general shape are placed inside each cell such that the local-cell particle packing fractions are identical to the global packing fraction. We prove that the constructed packings with a polydispersity in size are perfectly hyperuniform in the infinite-sample-size limit, regardless of particle shapes, positions, and numbers per cell. We use this theoretical formulation to devise an efficient and tunable algorithm to generate extremely large realizations of such packings. We employ two distinct initial tessellations: Voronoi as well as sphere tessellations. Beginning with Voronoi tessellations, we show that our algorithm can remarkably convert extremely large nonhyperuniform packings into hyperuniform ones in R^{2} and R^{3}. Implementing our theoretical methodology on sphere tessellations, we establish the hyperuniformity of the classical Hashin-Shtrikman multiscale coated-spheres structures, which are known to be two-phase media microstructures that possess optimal effective transport and elastic properties. A consequence of our work is a rigorous demonstration that packings that have identical tessellations can either be nonhyperuniform or hyperuniform by simply tuning local characteristics. It is noteworthy that our computationally designed hyperuniform two-phase systems can easily be fabricated via state-of-the-art methods, such as 2D photolithographic and 3D printing technologies. In addition, the tunability of our methodology offers a route for the discovery of novel disordered hyperuniform two-phase materials.

Probing Quantum Turbulence in ^{4}He by Quantum Evaporation Measurements

Theory of superfluid ^{4}He shows that, due to strong correlations and backflow effects, the density profile of a vortex line has the character of a density modulation and it is not a simple rarefaction region as found in clouds of cold bosonic atoms. We find that the basic features of this density modulation are represented by a wave packet of cylindrical symmetry in which rotons with a positive group velocity have a dominant role: The vortex density modulation can be viewed as a cloud of virtual excitations, mainly rotons, sustained by the phase of the vortex wave function. This suggests that in a vortex reconnection some of these rotons become real so that a vortex tangle is predicted to be a source of nonthermal rotons. The presence of such vorticity induced rotons can be verified by measurements at low temperature of quantum evaporation of ^{4}He atoms. We estimate the rate of evaporation and this turns out to be detectable by current instrumentation. Additional information on the microscopic processes in the decay of quantum turbulence will be obtained if quantum evaporation by high energy phonons should be detected.

Metal-Insulator Transition of Solid Hydrogen by the Antisymmetric Shadow Wave Function

We revisit the pressure-induced molecular-atomic metal-insulator transition of solid hydrogen by means of variational quantum Monte Carlo simulations based on the antisymmetric shadow wave function. For the purpose of facilitating the study of the electronic structure of large-scale fermionic systems, the shadow wave function formalism is extended by a series of technical advancements as implemented in our HswfQMC code. Among others, these improvements include a revised optimization method for the employed shadow wave function and an enhanced treatment of periodic systems with long-range interactions. It is found that the superior accuracy of the antisymmetric shadow wave function results in a significantly increased transition pressure with respect to previous theoretical estimates.

Nonlinear Network Description for Many-Body Quantum Systems in Continuous Space

We show that the recently introduced iterative backflow wave function can be interpreted as a general neural network in continuum space with nonlinear functions in the hidden units. Using this wave function in variational Monte Carlo simulations of liquid ^{4}He in two and three dimensions, we typically find a tenfold increase in accuracy over currently used wave functions. Furthermore, subsequent stages of the iteration procedure define a set of increasingly good wave functions, each with its own variational energy and variance of the local energy: extrapolation to zero variance gives energies in close agreement with the exact values. For two dimensional ^{4}He, we also show that the iterative backflow wave function can describe both the liquid and the solid phase with the same functional form-a feature shared with the shadow wave function, but now joined by much higher accuracy. We also achieve significant progress for liquid ^{3}He in three dimensions, improving previous variational and fixed-node energies.

Rational design of stealthy hyperuniform two-phase media with tunable order

Disordered stealthy hyperuniform materials are exotic amorphous states of matter that have attracted recent attention because of their novel structural characteristics (hidden order at large length scales) and physical properties, including desirable photonic and transport properties. It is therefore useful to devise algorithms that enable one to design a wide class of such amorphous configurations at will. In this paper, we present several algorithms enabling the systematic identification and generation of discrete (digitized) stealthy hyperuniform patterns with a tunable degree of order, paving the way towards the rational design of disordered materials endowed with novel thermodynamic and physical properties. To quantify the degree of order or disorder of the stealthy systems, we utilize the discrete version of the τ order metric, which accounts for the underlying spatial correlations that exist across all relevant length scales in a given digitized two-phase (or, equivalently, a two-spin state) system of interest. Our results impinge on a myriad of fields, ranging from physics, materials science and engineering, visual perception, and information theory to modern data science.

Excitations in the quantum liquid 4He: A review

Progress made in measuring and interpreting the elementary excitations of superfluid and normal liquid [Formula: see text] in the past 25 years is reviewed. The goal is to bring up to date the data, calculations and our understanding of the excitations since the books and reviews of the early 1990s. Only bulk liquid [Formula: see text] is considered. Reference to liquid [Formula: see text], mixtures, reduced dimensions (films and confined helium) is made where useful to enhance interpretation. The focus is on the excitations as measured by inelastic neutron scattering methods. The review covers the dynamical response of liquid [Formula: see text] from the collective excitations at low energy and long wavelength (i.e. phonon-roton modes) to the single particle excitations at high energy from which the atomic momentum distribution and Bose-Einstein condensate fraction are determined. A goal is to show the interplay of these excitations with other spectacular properties such as superfluidity and the test of fundamental calculations of quantum liquids that is possible. The role of Bose-Einstein condensation in determining the nature of the [Formula: see text] mode and particularly it's temperature dependence is emphasized. The similarity of normal liquid [Formula: see text] with other quantum and classical liquids is discussed.

Microstructure and mechanical properties of hyperuniform heterogeneous materials

A hyperuniform random heterogeneous material is one in which the local volume fraction fluctuations in an observation window decay faster than the reciprocal window volume as the window size increases. Recent studies show that this class of materials are endowed with superior physical properties such as large isotropic photonic band gaps and optimal transport properties. Here we employ a stochastic optimization procedure to systematically generate realizations of hyperuniform heterogeneous materials with controllable short-range order, which is partially quantified using the two-point correlation function S_{2}(r) associated with the phase of interest. Specifically, our procedure generalizes the widely used Yeong-Torquato reconstruction procedure by including an additional constraint for hyperuniformity, i.e., the volume integral of the autocovariance function χ(r)=S_{2}(r)-ϕ^{2} over the whole space is zero. In addition, we only require the reconstructed S_{2} to match the target function up to a certain cutoff distance γ, in order to give the system sufficient degrees of freedom to satisfy the hyperuniform condition. By systematically increasing the γ value for a given S_{2}, one can produce a spectrum of hyperuniform heterogeneous materials with varying degrees of partial short-range order compatible with the specified S_{2}. The mechanical performance including both elastic and brittle fracture behaviors of the generated hyperuniform materials is analyzed using a volume-compensated lattice-particle method. For the purpose of comparison, the corresponding nonhyperuniform materials with the same short-range order (i.e., with S_{2} constrained up to the same γ value) are also constructed and their mechanical performance is analyzed. Here we consider two specific S_{2} including the positive exponential decay function and the correlation function associated with an equilibrium hard-sphere system. For the constructed systems associated with these two specific functions, we find that although the hyperuniform materials are softer than their nonhyperuniform counterparts, the former generally possess a significantly higher brittle fracture strength than the latter. This superior mechanical behavior is attributed to the lower degree of stress concentration in the material resulting from the hyperuniform microstructure, which is crucial to crack initiation and propagation.

Structural Characterization and Statistical-Mechanical Model of Epidermal Patterns

In proliferating epithelia of mammalian skin, cells of irregular polygon-like shapes pack into complex, nearly flat two-dimensional structures that are pliable to deformations. In this work, we employ various sensitive correlation functions to quantitatively characterize structural features of evolving packings of epithelial cells across length scales in mouse skin. We find that the pair statistics in direct space (correlation function) and Fourier space (structure factor) of the cell centroids in the early stages of embryonic development show structural directional dependence (statistical anisotropy), which is a reflection of the fact that cells are stretched, which promotes uniaxial growth along the epithelial plane. In the late stages, the patterns tend toward statistically isotropic states, as cells attain global polarization and epidermal growth shifts to produce the skin's outer stratified layers. We construct a minimalist four-component statistical-mechanical model involving effective isotropic pair interactions consisting of hard-core repulsion and extra short-range soft-core repulsion beyond the hard core, whose length scale is roughly the same as the hard core. The model parameters are optimized to match the sample pair statistics in both direct and Fourier spaces. By doing this, the parameters are biologically constrained. In contrast with many vertex-based models, our statistical-mechanical model does not explicitly incorporate information about the cell shapes and interfacial energy between cells; nonetheless, our model predicts essentially the same polygonal shape distribution and size disparity of cells found in experiments, as measured by Voronoi statistics. Moreover, our simulated equilibrium liquid-like configurations are able to match other nontrivial unconstrained statistics, which is a testament to the power and novelty of the model. The array of structural descriptors that we deploy enable us to distinguish between normal, mechanically deformed, and pathological skin tissues. Our statistical-mechanical model enables one to generate tissue microstructure at will for further analysis. We also discuss ways in which our model might be extended to better understand morphogenesis (in particular the emergence of planar cell polarity), wound healing, and disease-progression processes in skin, and how it could be applied to the design of synthetic tissues.

Disordered hyperuniform heterogeneous materials

Disordered hyperuniform many-body systems are distinguishable states of matter that lie between a crystal and liquid: they are like perfect crystals in the way they suppress large-scale density fluctuations and yet are like liquids or glasses in that they are statistically isotropic with no Bragg peaks. These systems play a vital role in a number of fundamental and applied problems: glass formation, jamming, rigidity, photonic and electronic band structure, localization of waves and excitations, self-organization, fluid dynamics, quantum systems, and pure mathematics. Much of what we know theoretically about disordered hyperuniform states of matter involves many-particle systems. In this paper, we derive new rigorous criteria that disordered hyperuniform two-phase heterogeneous materials must obey and explore their consequences. Two-phase heterogeneous media are ubiquitous; examples include composites and porous media, biological media, foams, polymer blends, granular media, cellular solids, and colloids. We begin by obtaining some results that apply to hyperuniform two-phase media in which one phase is a sphere packing in d-dimensional Euclidean space [Formula: see text]. Among other results, we rigorously establish the requirements for packings of spheres of different sizes to be 'multihyperuniform'. We then consider hyperuniformity for general two-phase media in [Formula: see text]. Here we apply realizability conditions for an autocovariance function and its associated spectral density of a two-phase medium, and then incorporate hyperuniformity as a constraint in order to derive new conditions. We show that some functional forms can immediately be eliminated from consideration and identify other forms that are allowable. Specific examples and counterexamples are described. Contact is made with well-known microstructural models (e.g. overlapping spheres and checkerboards) as well as irregular phase-separation and Turing-type patterns. We also ascertain a family of autocovariance functions (or spectral densities) that are realizable by disordered hyperuniform two-phase media in any space dimension, and present select explicit constructions of realizations. These studies provide insight into the nature of disordered hyperuniformity in the context of heterogeneous materials and have implications for the design of such novel amorphous materials.