Activation of indistinguishability-based quantum coherence for enhanced metrological applications with particle statistics imprint

Quantum coherence, an essential feature of quantum mechanics allowing quantum superposition of states, is a resource for quantum information processing. Coherence emerges in a fundamentally different way for nonidentical and identical particles. For the latter, a unique contribution exists linked to indistinguishability that cannot occur for nonidentical particles. Here we experimentally demonstrate this additional contribution to quantum coherence with an optical setup, showing that its amount directly depends on the degree of indistinguishability and exploiting it in a quantum phase discrimination protocol. Furthermore, the designed setup allows for simulating fermionic particles with photons, thus assessing the role of exchange statistics in coherence generation and utilization. Our experiment proves that independent indistinguishable particles can offer a controllable resource of coherence and entanglement for quantum-enhanced metrology.

A quantum system can reside in coherent superpositions of states, which have a role in the interpretation of quantum mechanics (1–4), lead to nonclassicality (5, 6), and imply the intrinsically probabilistic nature of predictions in the quantum realm (7, 8). Besides this fundamental role, quantum coherence is also at the basis of quantum algorithms (9–14) and, from a modern information-theoretic perspective, constitutes a paradigmatic basis-dependent quantum resource (15–17), providing a quantifiable advantage in certain quantum information protocols.For a single quantum particle, coherence manifests itself when the particle is found in a superposition of a reference basis, for instance, the computational basis of the Hilbert space. Formally, any quantum state whose density matrix contains nonzero diagonal elements when expressed in the reference basis is said to display quantum coherence (16). This is the definition of quantum coherence employed in our work. For multiparticle compound systems, the physics underlying the emergence of quantum coherence is richer and strictly connected to the nature of the particles, with fundamental differences for nonidentical and identical particles. A particularly intriguing observation is that the states of identical particle systems can manifest coherence even when no particle resides in superposition states, provided that the wave functions of the particles overlap (18–20). In general, a special contribution to quantum coherence arises thanks to the spatial indistinguishability of identical particles, which cannot exist for nonidentical (or distinguishable) particles (18). Recently, it has been found that the spatial indistinguishability of identical particles can be exploited for entanglement generation (21), applicable even for spacelike-separated quanta (22) and against preparation and dynamical noises (23–26). The presence of entanglement is a signature that the bipartite system as a whole carries coherence even when the individual particles do not, the amount of this coherence being dependent on the degree of indistinguishability. We name this specific contribution to quantumness of compound systems “indistinguishability-based coherence,” in contrast to the more familiar “single-particle superposition-based coherence.” Indistinguishability-based coherence qualifies in principle as an exploitable resource for quantum metrology (18). However, it requires sophisticated control techniques to be harnessed, especially in view of its nonlocal nature. Moreover, a crucial property of identical particles is the exchange statistics, while its experimental study requiring operating both bosons and fermions in the same setup is generally challenging.In the present work, we investigate the operational contribution of quantum coherence stemming from the spatial indistinguishability of identical particles. The main aim of our experiment is to prove that elementary states of two independent spatially indistinguishable particles can give rise to exploitable quantum coherence, with a measurable effect due to particle statistics. By utilizing our recently developed photonic architecture capable of tuning the indistinguishability of two uncorrelated photons (27), we observe the direct connection between the degree of indistinguishability and the amount of generated coherence and show that indistinguishability-based coherence can be concurrent with single-particle superposition-based coherence. In particular, we demonstrate its operational implications, namely, providing a quantifiable advantage in a phase discrimination task (28, 29), as depicted in Fig. 1. Furthermore, we design a setup capable of testing the impact of particle statistics in coherence production and phase discrimination for both bosons and fermions; this is accomplished by compensating for the exchange phase during state preparation, simulating fermionic states with photons, which leads to statistics-dependent efficiency of the quantum task.Open in a separate windowFig. 1.Illustration of the indistinguishability-activated phase discrimination task. A resource state ρin that contains coherence in a computational basis is generated from spatial indistinguishability. The state then enters a black box which implements a phase unitary U^k=eiG^ϕk,k∈{1,…,n} on ρin. The goal is to determine the ϕk actually applied through the output state ρout: indistinguishability-based coherence provides an operational advantage in this task.     

Layering transitions and metastable structures of cholesteric liquid crystals in cylindrical confinement

Our study of cholesteric lyotropic chromonic liquid crystals in cylindrical confinement reveals the topological aspects of cholesteric liquid crystals. The double-twist configurations we observe exhibit discontinuous layering transitions, domain formation, metastability, and chiral point defects as the concentration of chiral dopant is varied. We demonstrate that these distinct layer states can be distinguished by chiral topological invariants. We show that changes in the layer structure give rise to a chiral soliton similar to a toron, comprising a metastable pair of chiral point defects. Through the applicability of the invariants we describe to general systems, our work has broad relevance to the study of chiral materials.

Chiral liquid crystals (LCs) are ubiquitous, useful, and rich systems (1–4). From the first discovery of the liquid crystalline phase to the variety of chiral structures formed by biomolecules (5–9), the twisted structure, breaking both mirror and continuous spatial symmetries, is omnipresent. The unique structure also makes the chiral nematic (cholesteric) LC, an essential material for applications utilizing the tunable, responsive, and periodic modulation of anisotropic properties.The cholesteric is also a popular model system to study the geometry and topology of partially ordered matter. The twisted ground state of the cholesteric is often incompatible with confinement and external fields, exhibiting a large variety of frustrated and metastable director configurations accompanying topological defects. Besides the classic example of cholesterics in a Grandjean−Cano wedge (10, 11), examples include cholesteric droplets (12–16), colloids (17–19), shells (20–22), tori (23, 24), cylinders (25–29), microfabricated structures (30, 31), and films between parallel plates with external fields (32–40). These structures are typically understood using a combination of nematic (achiral) topology (41, 42) and energetic arguments, for example, the highly successful Landau−de Gennes approach (43). However, traditional extensions of the nematic topological approach to cholesterics are known to be conceptually incomplete and difficult to apply in regimes where the system size is comparable to the cholesteric pitch (41, 44).An alternative perspective, chiral topology, can give a deeper understanding of these structures (45–47). In this approach, the key role is played by the twist density, given in terms of the director field n by n⋅∇×n. This choice is not arbitrary; the Frank free energy prefers n⋅∇×n≈−q0=2π/p0 with a helical pitch p0, and, from a geometric perspective, n⋅∇×n≠0 defines a contact structure (48). This allows a number of new integer-valued invariants of chiral textures to be defined (45). A configuration with a single sign of twist is chiral, and two configurations which cannot be connected by a path of chiral configurations are chirally distinct, and hence separated by a chiral energy barrier. Within each chiral class of configuration, additional topological invariants may be defined using methods of contact topology (45–48), such as layer numbers. Changing these chiral topological invariants requires passing through a nonchiral configuration. Cholesterics serve as model systems for the exploration of chirality in ordered media, and the phenomena we describe here—metastability in chiral systems controlled by chiral topological invariants—has applicability to chiral order generally. This, in particular, includes chiral ferromagnets, where, for example, our results on chiral topological invariants apply to highly twisted nontopological Skyrmions (49, 50) (“Skyrmionium”).Our experimental model to explore the chiral topological invariants is the cholesteric phase of lyotropic chromonic LCs (LCLCs). The majority of experimental systems hitherto studied are based on thermotropic LCs with typical elastic and surface-anchoring properties. The aqueous LCLCs exhibiting unusual elastic properties, that is, very small twist modulus K2 and large saddle-splay modulus K24 (51–56), often leading to chiral symmetry breaking of confined achiral LCLCs (53, 54, 56–61), may enable us to access uncharted configurations and defects of topological interests. For instance, in the layer configuration by cholesteric LCLCs doped with chiral molecules, their small K2 provides energetic flexibility to the thickness of the cholesteric layer, that is, the repeating structure where the director n twists by π. The large K24 affords curvature-induced surface interactions in combination with a weak anchoring strength of the lyotropic LCs (62–64).We present a systematic investigation of the director configuration of cholesteric LCLCs confined in cylinders with degenerate planar anchoring, depending on the chiral dopant concentration. We show that the structure of cholesteric configurations is controlled by higher-order chiral topological invariants. We focus on two intriguing phenomena observed in cylindrically confined cholesterics. First, the cylindrical symmetry renders multiple local minima to the energy landscape and induces discontinuous increase of twist angles, that is, a layering transition, upon the dopant concentration increase. Additionally, the director configurations of local minima coexist as metastable domains with point-like defects between them. We demonstrate that a chiral layer number invariant distinguishes these configurations, protects the distinct layer configurations (45), and explains the existence of the topological defect where the invariant changes.     

Transient conformational fluctuation of TePixD during a reaction

Knowledge of the dynamical behavior of proteins, and in particular their conformational fluctuations, is essential to understanding the mechanisms underlying their reactions. Here, transient enhancement of the isothermal partial molar compressibility, which is directly related to the conformational fluctuation, during a chemical reaction of a blue light sensor protein from the thermophilic cyanobacterium Thermosynechococcus elongatus BP-1 (TePixD, Tll0078) was investigated in a time-resolved manner. The UV-Vis absorption spectrum of TePixD did not change with the application of high pressure. Conversely, the transient grating signal intensities representing the volume change depended significantly on the pressure. This result implies that the compressibility changes during the reaction. From the pressure dependence of the amplitude, the compressibility change of two short-lived intermediate (I₁ and I₂) states were determined to be +(5.6 ± 0.6) × 10⁻² cm³⋅mol⁻¹⋅MPa⁻¹ for I₁ and +(6.6 ± 0.7)×10⁻² cm³⋅mol⁻¹⋅MPa⁻¹ for I₂. This result showed that the structural fluctuation of intermediates was enhanced during the reaction. To clarify the relationship between the fluctuation and the reaction, the compressibility of multiply excited TePixD was investigated. The isothermal compressibility of I₁ and I₂ intermediates of TePixD showed a monotonic decrease with increasing excitation laser power, and this tendency correlated with the reactivity of the protein. This result indicates that the TePixD decamer cannot react when its structural fluctuation is small. We concluded that the enhanced compressibility is an important factor for triggering the reaction of TePixD. To our knowledge, this is the first report showing enhanced fluctuations of intermediate species during a protein reaction, supporting the importance of fluctuations.Proteins often transfer information through changes in domain–domain (or intermolecular) interactions. Photosensor proteins are an important example. They have light-sensing domains and function by using the light-driven changes in domain–domain interactions (1). The sensor of blue light using FAD (BLUF) domain is a light-sensing module found widely among the bacterial kingdom (2). The BLUF domain initiates its photoreaction by the light excitation of the flavin moiety inside the protein, which changes the domain–domain interaction, causing a quaternary structural change and finally transmitting biological signals (3, 4). It has been an important research topic to elucidate how the initial photochemistry occurring in the vicinity of the chromophore leads to the subsequent large conformation change in other domains, which are generally apart from the chromophore.It may be reasonable to consider that the conformation change in the BLUF domain is the driving force in its subsequent reaction; that is, the change in domain–domain interaction. However, sometimes, clear conformational changes have not been observed for the BLUF domain; its conformation is very similar before and after photo-excitation (5–13). The circular dichroism (CD) spectra of BLUF proteins AppA and PixD from thermophilic cyanobacterium Thermosynechococcus elongatus BP-1 (TePixD) did not change on illumination (5, 13). Similarly, solution NMR studies of AppA and BlrB showed only small chemical shifts on excitation (9, 10). The solution NMR structure of BlrP1 showed a clear change, but this was limited in its C-terminal extension region and not core BLUF (11). Furthermore, the diffusion coefficient (D) of the BLUF domain of YcgF was not changed by photo-excitation (12), although D is sensitive to global conformational changes. These results imply that a minor structural change occurs in the BLUF domain. In such cases, how does the BLUF domain control its interdomain interaction? Recently, a molecular dynamics (MD) simulation on another light-sensing domain, the light-oxygen-voltage (LOV) sensing domain, suggested that fluctuation of the LOV core structure could be a key to understanding the mechanism of information transfer (14–16).Because proteins work at room temperature, they are exposed to thermal fluctuations. The importance of such structural fluctuations for biomolecular reactions has been also pointed out: for example, enzymatic activity (17–20). Experimental detections of such conformation fluctuations using single molecular detection (21) or NMR techniques such as the hydrogen-deuterium (H-D) exchange, relaxation dispersion method, and high-pressure NMR (22–24) have succeeded. However, these techniques could not detect the fluctuation of short-lived transient species. Indeed, single molecule spectroscopy can trace the fluctuation in real time, but it is still rather difficult to detect rapid fluctuations for a short-lived intermediate during a reaction. Therefore, information about the fluctuation of intermediates is thus far limited.A thermodynamic measurement is another way to characterize the fluctuation of proteins. In particular, the partial molar isothermal compressibility [K¯T=−(∂V¯/∂P)T] is essential, because this property is directly linked to the mean-square fluctuations of the protein partial molar volume by 〈(V¯−〈V¯〉)2〉≡〈δV¯2〉=kBTK¯T (25). (Here, <X> means the averaged value of a quantity of X.) Therefore, isothermal compressibility is thought to reflect the structural fluctuation of molecules (26). However, experimental measurement of this parameter of proteins in a dilute solution is quite difficult. Indeed, this quantity has been determined indirectly from the theoretical equation using the adiabatic compressibility of a protein solution, which was determined by the sound velocity in the solution (26–31). Although the relation between volume fluctuations and isothermal compressibility is rigorously correct only with respect to the intrinsic part of the volume compressibility, and not the partial molar volume compressibility (32), we considered that this partial molar volume compressibility is still useful for characterizing the fluctuation of the protein structure including its interacting water molecules. In fact, the relationship between β¯T and the volume fluctuation has been often used to discuss the fluctuation of proteins (17, 26–28), and the strong correlation of β¯T of reactants with the functioning for some enzymes (17, 33, 34) has been reported. These studies show the functional importance of the structural fluctuation represented by β¯T. However, thermodynamic techniques lack time resolution, and it has been impossible to measure the fluctuations of short-lived intermediate species.Recently, we developed a time-resolving method for assessing thermodynamic properties using the pulsed laser induced transient grating (TG) method. Using this method, we thus far succeeded in measuring the enthalpy change (ΔH) (35–38), partial molar volume change (ΔV¯) (12, 35, 37), thermal expansion change (Δα¯th) (12, 37), and heat capacity change (ΔC_p) (36–38) for short-lived species. Therefore, in principle, the partial molar isothermal compressibility change (ΔK¯T) of a short-lived intermediate become observable if we conduct the TG experiment under the high-pressure condition and detect ΔV¯ with varying external pressure.There are several difficulties in applying the traditional high-pressure cell to the TG method to measure thermodynamic parameters quantitatively. The most serious problem is ensuring the quantitative performance of the intensity of TG signals measured under the high-pressure condition. On this point, our group has developed a new high-pressure cell specially designed for TG spectroscopy (39) and overcome this problem. In this paper, by applying this high-pressure TG system to the BLUF protein TePixD, we report the first measurement, to our knowledge, of ΔK¯T of short-lived intermediates to investigate the mechanism underlying signal transmission by BLUF proteins, from the view point of the transient fluctuation.TePixD is a homolog of the BLUF protein PixD, which regulates the phototaxis of cyanobacterium (40) and exists in a thermophilic cyanobacterium Thermocynechococcus elongates BP-1 (Tll0078). TePixD is a relatively small (17 kDa) protein that consists only of the BLUF domain with two extended helices in the C-terminal region. In crystals and solutions, it forms a decamer that consists of two pentameric rings (41). The photochemistry of TePixD is typical among BLUF proteins (42–45); on blue light illumination, the absorption spectrum shifts toward red by about 10 nm within a nanosecond. The absorption spectrum does not change further, and the dark state is recovered with a time constant of ∼5 s at room temperature (40, 43). The spectral red shift was explained by the rearrangement of the hydrogen bond network around the chromophore (6, 46–48). The TG method has revealed the dynamic photoreaction mechanism, which cannot be detected by conventional spectroscopic methods. The TG signal of TePixD (Fig. S1) showed that there are two spectrally silent reaction phases: a partial molar volume expansion with the time constant of ∼40 μs and the diffusion coefficient (D) change with a time constant of ∼4 ms. Furthermore, it was reported that the pentamer and decamer states of TePixD are in equilibrium and that the final photoproduct of the decamer is pentamers generated by its dissociation (13, 49). On the basis of these studies, the reaction scheme has been identified as shown in Fig. 1. Here, I₁ is the intermediate of the spectrally red-shifted species (generated within a nanosecond) and I₂ is the one created on the subsequent volume expansion process of +4 cm³⋅mol⁻¹ (∼40 μs). Furthermore, an experiment of the excitation laser power dependence of its TG signal revealed that the TePixD decamer undergoes the original dissociation reaction when only one monomer in the decamer is excited (50). In this study, we investigated the transient compressibility of the intermediates I₁ and I₂ of the photoreaction of TePixD and found a direct link between their fluctuation and reactivity.Open in a separate windowFig. 1.Schematic illustration of the photoreaction of TePixD. Yellow circles represent the TePixD monomer in the ground state, which constructs the decamer and pentamer states. In the dark state, these two forms are in equilibrium. The excited, spectral red-shifted state of the TePixD monomer is indicated by a red circle. The square represents the I₂ state of the monomer, which is created by the volume expansion process.     

Optimized river diversion scenarios promote sustainability of urbanized deltas

Socioeconomic viability of fluvial-deltaic systems is limited by natural processes of these dynamic landforms. An especially impactful occurrence is avulsion, whereby channels unpredictably shift course. We construct a numerical model to simulate artificial diversions, which are engineered to prevent channel avulsion, and direct sediment-laden water to the coastline, thus mitigating land loss. We provide a framework that identifies the optimal balance between river diversion cost and civil disruption by flooding. Diversions near the river outlet are not sustainable, because they neither reduce avulsion frequency nor effectively deliver sediment to the coast; alternatively, diversions located halfway to the delta apex maximize landscape stability while minimizing costs. We determine that delta urbanization generates a positive feedback: infrastructure development justifies sustainability and enhanced landform preservation vis-à-vis diversions.

Deltaic environments are critical for societal wellbeing because these landscapes provide an abundance of natural resources that promote human welfare (1, 2). However, the sustainability of deltas is uncertain due to sea-level rise (3, 4), sediment supply reduction (4–6), and land subsidence (7, 8). Additionally, river avulsion, the process of sudden channel relocation (9, 10), presents a dichotomy to delta sustainability: the unanticipated civil disruption associated with flooding brought by channel displacement is at odds with society’s desire for landscape stability, yet channel relocation is needed to deliver nutrients and sediment to various locations along the deltaic coastline (11, 12). Indeed, for many of the world’s megadeltas, channel engineering practices have sought to restrict channel mobility and limit floodplain connectivity (13, 14), which in turn prevents sediment dispersal that is necessary to sustain deltas; as a consequence, land loss has ensued (15). Despite providing near-term stability (13–15), engineering of deltaic channels is a long-term detrimental practice (11, 15–17).To maximize societal benefit, measures that promote delta sustainability must balance engineering infrastructure cost and impact on delta morphology with benefits afforded by maintaining and developing deltaic landscapes (1, 2, 11, 12, 16–19). For example, channel diversions, costing millions to billions of dollars (20–22), are now planned worldwide to both prevent unintended avulsions and ensure coastal sustainability through enhanced sediment delivery (e.g., Fig. 1A) (20, 21, 23–26).Open in a separate windowFig. 1.(A) Satellite image of Yellow River delta (Landsat, 1978) showing coastline response to a diversion in 1976 at the open circle, which changed the channel course from the north (Diaokou lobe) to the east (Qingshuigou lobe) and produced flooding over the stripe-hatched area (30). (B and C) Planform view (B) and along-channel cross-section view (C) of conceptual model for numerical simulations and societal benefit formulation. In the diagrams, a diversion at LD≈0.8Lb floods an area (af) defined by Lf and θ, diverting sediment away from the deltaic lobe (with length Ll). Aggradation of the former channel bed (dashed line) is variable; hence, diversion length influences the propensity for subsequent avulsion setup.In this article, we consider the benefits and costs of such engineered river diversions and determine how these practices most effectively sustain deltaic landscapes, by assessing optimal placement and timing for river diversions. Addressing these points requires combining two modeling frameworks: a morphodynamic approach—evolving the landscape over time and space by evaluating the interactions of river fluid flow and sediment transport—and a decision-making framework (21, 22, 27, 28). The former simulates deltaic channel diversions by assessing the nonlinear relationships between channel diversion length (LD) and the frequency (timing) of avulsions (TA), while the latter incorporates a societal benefit model that approximates urbanization by considering the cost of flooding a landscape that would otherwise generate revenue. The aim is to optimize timing and placement of channel diversions, by giving consideration to morphodynamic operations and societal wellbeing. Interestingly, optimal societal benefit indicates that urbanization justifies enhanced sustainability measures, which contradicts existing paradigms that label development and sustainability mutually exclusive (3, 7, 12). Ultimately, the societal benefit model should be an integrated component in decision-making frameworks. This will help locate diversions and promote sustainable and equitable decisions considering historical, ethical, and environmental contexts for river management decisions (29).     

Nature of dynamic gradients,glass formation,and collective effects in ultrathin freestanding films

Molecular, polymeric, colloidal, and other classes of liquids can exhibit very large, spatially heterogeneous alterations of their dynamics and glass transition temperature when confined to nanoscale domains. Considerable progress has been made in understanding the related problem of near-interface relaxation and diffusion in thick films. However, the origin of “nanoconfinement effects” on the glassy dynamics of thin films, where gradients from different interfaces interact and genuine collective finite size effects may emerge, remains a longstanding open question. Here, we combine molecular dynamics simulations, probing 5 decades of relaxation, and the Elastically Cooperative Nonlinear Langevin Equation (ECNLE) theory, addressing 14 decades in timescale, to establish a microscopic and mechanistic understanding of the key features of altered dynamics in freestanding films spanning the full range from ultrathin to thick films. Simulations and theory are in qualitative and near-quantitative agreement without use of any adjustable parameters. For films of intermediate thickness, the dynamical behavior is well predicted to leading order using a simple linear superposition of thick-film exponential barrier gradients, including a remarkable suppression and flattening of various dynamical gradients in thin films. However, in sufficiently thin films the superposition approximation breaks down due to the emergence of genuine finite size confinement effects. ECNLE theory extended to treat thin films captures the phenomenology found in simulation, without invocation of any critical-like phenomena, on the basis of interface-nucleated gradients of local caging constraints, combined with interfacial and finite size-induced alterations of the collective elastic component of the structural relaxation process.

Spatially heterogeneous dynamics in glass-forming liquids confined to nanoscale domains (1–7) play a major role in determining the properties of molecular, polymeric, colloidal, and other glass-forming materials (8), including thin films of polymers (9, 10) and small molecules (11–15), small-molecule liquids in porous media (2, 4, 16, 17), semicrystalline polymers (18, 19), polymer nanocomposites (20–22), ionomers (23–25), self-assembled block and layered (26–33) copolymers, and vapor-deposited ultrastable molecular glasses (34–36). Intense interest in this problem over the last 30 y has also been motivated by the expectation that its understanding could reveal key insights concerning the mechanism of the bulk glass transition.Considerable progress has been made for near-interface altered dynamics in thick films, as recently critically reviewed (1). Large amplitude gradients of the structural relaxation time, τ(z,T), converge to the bulk value, τ_bulk(T), in an intriguing double-exponential manner with distance, z, from a solid or vapor interface (13, 37–42). This implies that the corresponding effective activation barrier, F_total(z,T,H) (where H is film thickness), varies exponentially with z, as does the glass transition temperature, T_g (37). Thus the fractional reduction in activation barrier, ε(z,H), obeys the equation ε(z,H)≡1−Ftotal(z,T,H)/Ftotal,bulk(T)=ε0⁡exp(−z/ξF), where F_total,bulk(T) is the bulk temperature-dependent barrier and ξ_F a length scale of modest magnitude. Although the gradient of reduction in absolute activation barriers becomes stronger with cooling, the amplitude of the fractional reduction of the barrier gradient, quantified by ε₀, and the range ξ_F of this gradient, exhibit a weak or absent temperature dependence at the lowest temperatures accessed by simulations (typically with the strength of temperature dependence of ξ_F decreasing rather than increasing on cooling), which extend to relaxation timescales of order 10⁵ ps. This finding raises questions regarding the relevance of critical-phenomena–like ideas for nanoconfinement effects (1). Partially due to this temperature invariance, coarse-grained and all-atom simulations (1, 37, 42, 43) have found a striking empirical fractional power law decoupling relation between τ(z,T) and τ_bulk(T):τ(T,z)τbulk(T)∝(τbulk(T))−ε(z).[1]Recent theoretical analysis suggests (44) that this behavior is consistent with a number of experimental data sets as well (45, 46). Eq. 1 also corresponds to a remarkable factorization of the temperature and spatial location dependences of the barrier:Ftotal(z,T)=[1−ε(z)]Ftotal,bulk(T).[2]This finding indicates that the activation barrier for near-interface relaxation can be factored into two contributions: a z-dependent, but T-independent, “decoupling exponent,” ε(z), and a temperature-dependent, but position-insensitive, bulk activation barrier, F_total,bulk(T). Eq. 2 further emphasizes that ε(z) is equivalent to an effective fractional barrier reduction factor (for a vapor interface), 1−Ftotal(z,T,H)/Ftotal,bulk(T), that can be extracted from relaxation data.In contrast, the origin of “nanoconfinement effects” in thin films, and how much of the rich thick-film physics survives when dynamic gradients from two interfaces overlap, is not well understood. The distinct theoretical efforts for aspects of the thick-film phenomenology (44, 47–50) mostly assume an additive summation of one-interface effects in thin films, thereby ignoring possibly crucial cooperative and whole film finite size confinement effects. If the latter involve phase-transition–like physics as per recent speculations (14, 51), one can ask the following: do new length scales emerge that might be truncated by finite film size? Alternatively, does ultrathin film phenomenology arise from a combination of two-interface superposition of the thick-film gradient physics and noncritical cooperative effects, perhaps in a property-, temperature-, and/or thickness-dependent manner?Here, we answer these questions and establish a mechanistic understanding of thin-film dynamics for the simplest and most universal case: a symmetric freestanding film with two vapor interfaces. We focus on small molecules (modeled theoretically as spheres) and low to medium molecular weight unentangled polymers, which empirically exhibit quite similar alterations in dynamics under “nanoconfinement.” We do not address anomalous phenomena [e.g., much longer gradient ranges (29), sporadic observation of two distinct glass transition temperatures (52, 53)] that are sometimes reported in experiments with very high molecular weight polymers and which may be associated with poorly understood chain connectivity effects that are distinct from general glass formation physics (54–56).We employ a combination of molecular dynamics simulations with a zero-parameter extension to thin films of the Elastically Cooperative Nonlinear Langevin Equation (ECNLE) theory (57, 58). This theory has previously been shown to predict well both bulk activated relaxation over up to 14 decades (44–46) and the full single-gradient phenomenology in thick films (1). Here, we extend this theory to treat films of finite thickness, accounting for coupled interface and geometric confinement effects. We compare predictions of ECNLE theory to our previously reported (37, 43) and new simulations, which focus on translational dynamics of films comprised of a standard Kremer–Grest-like bead-spring polymer model (see SI Appendix). These simulations cover a wide range of film thicknesses (H, from 4 to over 90 segment diameters σ) and extend to low temperatures where the bulk alpha time is ∼0.1 μs (10⁵ Lennard Jones time units τ_LJ).The generalized ECNLE theory is found to be in agreement with simulation for all levels of nanoconfinement. We emphasize that this theory does not a priori assume any of the empirically established behaviors discovered using simulation (e.g., fractional power law decoupling, double-exponential barrier gradient, gradient flattening) but rather predicts these phenomena based upon interfacial modifications of the two coupled contributions to the underlying activation barrier– local caging constraints and a long-ranged collective elastic field. It is notable that this strong agreement is found despite the fact the dynamical ideas are approximate, and a simple hard sphere fluid model is employed in contrast to the bead-spring polymers employed in simulation. The basic unit of length in simulation (bead size σ) and theory (hard sphere diameter d) are expected to be proportional to within a prefactor of order unity, which we neglect in making comparisons.As an empirical matter, we find from simulation that many features of thin-film behavior can be described to leading order by a linear superposition of the thick-film gradients in activation barrier, that is:ε(z,H)=1−Ftotal(z,T,H)/Ftotal,bulk(T)≈ε0[exp(−z/ξF)+exp(−(H−z)/ξF)],[3]where the intrinsic decay length ξ_F is unaltered from its thick-film value and where ε₀ is a constant that, in the hypothesis of literal gradient additivity, is invariant to temperature and film thickness. We employ this functional form [originally suggested by Binder and coworkers (59)], which is based on a simple superposition of the two single-interface gradients, as a null hypothesis throughout this study: this form is what one expects if no new finite-size physics enters the thin-film problem relative to the thick film.However, we find that the superposition approximation progressively breaks down, and eventually entirely fails, in ultrathin films as a consequence of the emergence of a finite size confinement effect. The ECNLE theory predicts that this failure is not tied to a phase-transition–like mechanism but rather is a consequence of two key coupled physical effects: 1) transfer of surface-induced reduction of local caging constraints into the film, and 2) interfacial truncation and nonadditive modifications of the collective elastic contribution to the activation barrier.     

Empirical social triad statistics can be explained with dyadic homophylic interactions

The remarkable robustness of many social systems has been associated with a peculiar triangular structure in the underlying social networks. Triples of people that have three positive relations (e.g., friendship) between each other are strongly overrepresented. Triples with two negative relations (e.g., enmity) and one positive relation are also overrepresented, and triples with one or three negative relations are drastically suppressed. For almost a century, the mechanism behind these very specific (“balanced”) triad statistics remained elusive. Here, we propose a simple realistic adaptive network model, where agents tend to minimize social tension that arises from dyadic interactions. Both opinions of agents and their signed links (positive or negative relations) are updated in the dynamics. The key aspect of the model resides in the fact that agents only need information about their local neighbors in the network and do not require (often unrealistic) higher-order network information for their relation and opinion updates. We demonstrate the quality of the model on detailed temporal relation data of a society of thousands of players of a massive multiplayer online game where we can observe triangle formation directly. It not only successfully predicts the distribution of triangle types but also explains empirical group size distributions, which are essential for social cohesion. We discuss the details of the phase diagrams behind the model and their parameter dependence, and we comment on to what extent the results might apply universally in societies.

Recognizing the fundamental role of triadic interactions in shaping social structures, Heider (1) introduced the notion of balanced and unbalanced triads. A triad (triangle) of individuals is balanced if it includes zero or two negative links; otherwise, it is unbalanced. Heider (1) hypothesized that social networks have a tendency to reduce the number of unbalanced triangles over time such that balanced triads would dominate in a stationary situation. This theory of “social balance” has been confirmed empirically in many different contexts, such as schools (2), monasteries (3), social media (4), or computer games (5). Social balance theory and its generalizations (6–8) have been studied extensively for more than a half century for their importance in understanding polarization of societies (9), global organization of social networks (10), evolution of the network of international relations (11), opinion formation (12, 13), epidemic spreading (14, 15), government formation (16), and decision-making processes (17).Following Heider’s intuition (18–41), current approaches toward social balance often account for the effect of triangles on social network formation in one way or another. For example, the models in refs. 22 and 23 consider a reduction of the number of unbalanced triads either in the neighborhood of a node or in the whole network. The latter process sometimes leads to imbalance due to the existence of so-called jammed states (42). In order to reach social balance, individuals can also update their links according to their relations to common neighbors (18–21) or adjust link weights via opinion updates (24, 25) or via a minimization of social stress based on triadic interactions (37–44). These works not only ignore the difficulty of individuals to know the social interactions beyond their direct neighbors in reality, so far, they also have not considered the detailed statistical properties of the over- or underrepresentation of the different types of triads, such as those reported in refs. 4 and 5, with the exception of refs. 43 and 44.It is generally believed that the similarity of individuals plays a crucial role in the formation of social ties in social networks, something that has been called homophily (45–48). This means that to form a positive or negative tie with another person, people compare only pairwise overlaps in their individual opinions (dyadic interaction). It has also been argued that social link formation takes into account a tendency in people to balance their local interaction networks in the sense that they introduce friends to each other, that they do give up friendships if two mutual friends have negative attitudes toward each other, and that they tend to avoid situations where everyone feels negatively about the others. This is the essence of social balance theory (1). Obviously, link formation following social balance is cognitively much more challenging than homophily-based link formation since in the former, one has to keep in mind the many mutual relations between all your neighbors in a social network. While social balance–driven link formation certainly occurs in the context of close friendships, it is less realistic to assume that this mechanism is at work in social link formation in general. In Fig. 1, we schematically show the situation in a portion of a social network. It is generally hard for node i to know all the relations between his neighbors j, k, and l.Open in a separate windowFig. 1.Schematic view of opinion and link updates in a society. Every individual has an opinion vector whose components represent (binary) opinions on G=5 different subjects. Red (blue) links denote positive (negative) relationships. The question marks denote unknown relationships between i’s neighbors. As an agent i flips one of its opinions (red circle), si1, from 1 to –1, i can either decrease or increase its individual stress, H(i), depending on the value of the parameter α (Eq. 1). For instance, H(i) would increase if α=1 but would decrease for α=0. For high “rationality” values of individuals w.r.t. social stress, as quantified by β, the latter is more likely to be accepted, resulting in a reduction of the number of unbalanced triads in i’s neighborhood.Here, assuming that it is generally unrealistic for individuals to know their social networks at the triadic level, we aim to understand the emergence and the concrete statistics of balanced triads on the basis of dyadic or one-to-one interactions. Therefore, we use a classic homophily rule (45, 46) to define a “stress level” between any pair of individuals based on the similarity (or overlap) of their individual opinions. Here, the opinions of an individual i are represented by a vector with G components, si, that we show in Fig. 1. Homophily implies that i and j tend to become friends if the overlap (e.g., scalar product of their opinion vectors) is positive, and they become enemies if the overlap is negative. Such a specification of homophily is often referred to as an attraction–repulsion or assimilation–differentiation rule (49, 50). Assuming that, generally, social relations rearrange such as to minimize individual social stress on average, we will show that balanced triads naturally emerge from purely dyadic homophilic interactions without any explicit selection mechanisms for specific triads. We formulate the opinion link dynamics leading to social balance within a transparent physics-inspired framework. In particular, we observe a dynamic transition between two different types of balanced steady states that correspond to different compositions of balanced triads.Explaining the empirical statistics of triangles in social systems is a challenge. Early works considered groups of a few monks in a monastery (3) or a few students in classrooms (51). The studies suffered from limited data and small network sizes. Large-scale studies were first performed in online platforms (4) and in the society of players of the massive multiplayer online game (MMOG) Pardus. Players in Pardus engage in a form of economic life, such as trade and mining, and in social activities, such as communication on a number of channels, forming friendships and enmities (details are in refs. 5, 52, and 53). In the social networks of this game, balanced triads were once more confirmed to be overrepresented compared with what is expected by chance. Similar patterns of triad statistics were also observed in Epinion, Slashdot, and Wikipedia (4). More details on the Pardus society are in Materials and Methods. This dataset gives us the unique possibility to validate the model and compare the predictions with actual triangle statistics and formation of positively connected groups that are foundational to social cohesion.     

Dynamic rupture initiation and propagation in a fluid-injection laboratory setup with diagnostics across multiple temporal scales

Fluids are known to trigger a broad range of slip events, from slow, creeping transients to dynamic earthquake ruptures. Yet, the detailed mechanics underlying these processes and the conditions leading to different rupture behaviors are not well understood. Here, we use a laboratory earthquake setup, capable of injecting pressurized fluids, to compare the rupture behavior for different rates of fluid injection, slow (megapascals per hour) versus fast (megapascals per second). We find that for the fast injection rates, dynamic ruptures are triggered at lower pressure levels and over spatial scales much smaller than the quasistatic theoretical estimates of nucleation sizes, suggesting that such fast injection rates constitute dynamic loading. In contrast, the relatively slow injection rates result in gradual nucleation processes, with the fluid spreading along the interface and causing stress changes consistent with gradually accelerating slow slip. The resulting dynamic ruptures propagating over wetted interfaces exhibit dynamic stress drops almost twice as large as those over the dry interfaces. These results suggest the need to take into account the rate of the pore-pressure increase when considering nucleation processes and motivate further investigation on how friction properties depend on the presence of fluids.

The close connection between fluids and faulting has been revealed by a large number of observations, both in tectonic settings and during human activities, such as wastewater disposal associated with oil and gas extraction, geothermal energy production, and CO₂ sequestration (1–11). On and around tectonic faults, fluids also naturally exist and are added at depths due to rock-dehydration reactions (12–15) Fluid-induced slip behavior can range from earthquakes to slow, creeping motion. It has long been thought that creeping and seismogenic fault zones have little to no spatial overlap. Nonetheless, growing evidence suggests that the same fault areas can exhibit both slow and dynamic slip (16–19). The existence of large-scale slow slip in potentially seismogenic areas has been revealed by the presence of transient slow-slip events in subduction zones (16, 18) and proposed by studies investigating the physics of foreshocks (20–22).Numerical and laboratory modeling has shown that such complex fault behavior can result from the interaction of fluid-related effects with the rate-and-state frictional properties (9, 14, 19, 23, 24); other proposed rheological explanations for complexities in fault stability include combinations of brittle and viscous rheology (25) and friction-to-flow transitions (26). The interaction of frictional sliding and fluids results in a number of coupled and competing mechanisms. The fault shear resistance τres is typically described by a friction model that linearly relates it to the effective normal stress σ^n via a friction coefficient f:τres=f σ^n=f (σn−p),[1]where σn is the normal stress acting across the fault and p is the pore pressure. Clearly, increasing pore pressure p would reduce the fault frictional resistance, promoting the insurgence of slip. However, such slip need not be fast enough to radiate seismic waves, as would be characteristic of an earthquake, but can be slow and aseismic. In fact, the critical spatial scale h* for the slipping zone to reach in order to initiate an unstable, dynamic event is inversely proportional to the effective normal stress (27, 28) and hence increases with increasing pore pressure, promoting stable slip. This stabilizing effect of increasing fluid pressure holds for both linear slip-weakening and rate-and-state friction; it occurs because lower effective normal stress results in lower fault weakening during slip for the same friction properties. For example, the general form for two-dimensional (2D) theoretical estimates of this so-called nucleation size, h*, on rate-and-state faults with steady-state, velocity-weakening friction is given by:h*=(μ*DRS)/[F(a,b)(σn−p)],[2]where μ*=μ/(1−ν) for modes I and II, and μ*=μ for mode III (29); DRS is the characteristic slip distance; and F(a, b) is a function of the rate-and-state friction parameters a and b. The function F(a, b) depends on the specific assumptions made to obtain the estimate: FRR(a,b)=4(b−a)/π (ref. 27, equation 40) for a linearized stability analysis of steady sliding, or FRA(a,b)=[π(b−a)2]/2b, with a/b>1/2 for quasistatic crack-like expansion of the nucleation zone (ref. 30, equation 42).Hence, an increase in pore pressure induces a reduction in the effective normal stress, which both promotes slip due to lower frictional resistance and increases the critical length scale h*, potentially resulting in slow, stable fault slip instead of fast, dynamic rupture. Indeed, recent field and laboratory observations suggest that fluid injection triggers slow slip first (4, 9, 11, 31). Numerical modeling based on these effects, either by themselves or with an additional stabilizing effect of shear-layer dilatancy and the associated drop in fluid pressure, have been successful in capturing a number of properties of slow-slip events observed on natural faults and in field fluid-injection experiments (14, 24, 32–34). However, understanding the dependence of the fault response on the specifics of pore-pressure increase remains elusive. Several studies suggest that the nucleation size can depend on the loading rate (35–38), which would imply that the nucleation size should also depend on the rate of friction strength change and hence on the rate of change of the pore fluid pressure. The dependence of the nucleation size on evolving pore fluid pressure has also been theoretically investigated (39). However, the commonly used estimates of the nucleation size (Eq. 2) have been developed for faults under spatially and temporally uniform effective stress, which is clearly not the case for fluid-injection scenarios. In addition, the friction properties themselves may change in the presence of fluids (40–42). The interaction between shear and fluid effects can be further affected by fault-gauge dilation/compaction (40, 43–45) and thermal pressurization of pore fluids (42, 46–48).Recent laboratory investigations have been quite instrumental in uncovering the fundamentals of the fluid-faulting interactions (31, 45, 49–57). Several studies have indicated that fluid-pressurization rate, rather than injection volume, controls slip, slip rate, and stress drop (31, 49, 57). Rapid fluid injection may produce pressure heterogeneities, influencing the onset of slip. The degree of heterogeneity depends on the balance between the hydraulic diffusion rate and the fluid-injection rate, with higher injection rates promoting the transition from drained to locally undrained conditions (31). Fluid pressurization can also interact with friction properties and produce dynamic slip along rate-strengthening faults (50, 51).In this study, we investigate the relation between the rate of pressure increase on the fault and spontaneous rupture nucleation due to fluid injection by laboratory experiments in a setup that builds on and significantly develops the previous generations of laboratory earthquake setup of Rosakis and coworkers (58, 59). The previous versions of the setup have been used to study key features of dynamic ruptures, including sub-Rayleigh to supershear transition (60); rupture directionality and limiting speeds due to bimaterial effects (61); pulse-like versus crack-like behavior (62); opening of thrust faults (63); and friction evolution (64). A recent innovation in the diagnostics, featuring ultrahigh-speed photography in conjunction with digital image correlation (DIC) (65), has enabled the quantification of the full-field behavior of dynamic ruptures (66–68), as well as the characterization of the local evolution of dynamic friction (64, 69). In these prior studies, earthquake ruptures were triggered by the local pressure release due to an electrical discharge. This nucleation procedure produced only dynamic ruptures, due to the nearly instantaneous normal stress reduction.To study fault slip triggered by fluid injection, we have developed a laboratory setup featuring a hydraulic circuit capable of injecting pressurized fluid onto the fault plane of a specimen and a set of experimental diagnostics that enables us to detect both slow and fast fault slip and stress changes. The range of fluid-pressure time histories produced by this setup results in both quasistatic and dynamic rupture nucleation; the diagnostics allows us to capture the nucleation processes, as well as the resulting dynamic rupture propagation. In particular, here, we explore two injection techniques: procedure 1, a gradual, and procedure 2, a sharp fluid-pressure ramp-up. An array of strain gauges, placed on the specimen’s surface along the fault, can capture the strain (translated into stress) time histories over a wide range of temporal scales, spanning from microseconds to tens of minutes. Once dynamic ruptures nucleate, an ultrahigh-speed camera records images of the propagating ruptures, which are turned into maps of full-field displacements, velocities, and stresses by a tailored DIC) analysis. One advantage of using a specimen made of an analog material, such as poly(methyl meth-acrylate) (PMMA) used in this study, is its transparency, which allows us to look at the interface through the bulk and observe fluid diffusion over the interface. Another important advantage of using PMMA is that its much lower shear modulus results in much smaller nucleation sizes h* than those for rocks, allowing the experiments to produce both slow and fast slip in samples of manageable sizes.We start by describing the laboratory setup and the diagnostics monitoring the pressure evolution and the slip behavior. We then present and discuss the different slip responses measured as a result of slow versus fast fluid injection and interpret our measurements by using the rate-and-state friction framework and a pressure-diffusion model.     

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.RNAs perform critical cellular functions at the level of gene expression and regulation (1–4). RNA functions are determined not only by RNA structure or structure motifs [e.g., tetraloop hairpins (5, 6)] but also by conformational distributions and dynamics and kinetics of conformational changes. For example, riboswitches can adopt different conformations in response to specific conditions of the cellular environment (7, 8). Understanding the kinetics, such as the rate and pathways for the conformational changes, is critical for deciphering the mechanism of RNA function (9–19). Extensive experimental and theoretical studies on RNA folding kinetics have provided significant insights into the kinetic mechanism of RNA functions (19–36). However, due to the complexity of the RNA folding energy landscape (37–46) and the limitations of experimental tools (47–55), many fundamental problems, including single base flipping and base pair formation and fraying, remain unresolved. These unsolved fundamental problems have hampered our ability to resolve other important issues, such as RNA hairpin and larger structure folding kinetics. Several key questions remain unanswered, such as whether the hairpin folding is rate-limited by the conformational search of the native base pairs, whose formation leads to fast downhill folding of the whole structure, or by the breaking of misfolded base pairs before refolding to the native structure (18, 19, 54–73).Motivated by the need to understand the basic steps of nucleic acids folding, Hagan et al. (74) performed forty-three 200-ps unfolding trajectories at 400 K and identified both on- and off-pathway intermediates and two dominant unfolding pathways for a terminal C-G base pair in a DNA duplex. In one of the pathways, base pairing and stacking interactions are broken concomitantly, whereas in the other pathway, base stacking is broken after base pairing is disrupted. Furthermore, the unfolding requires that the Cyt diffuse away from the pairing Gua to a distance such that the C-G hydrogen bond cannot reform easily. More recently, Colizzi and Bussi (75) performed molecular dynamics (MD) pulling simulations for an RNA duplex and construct free energy landscape from the pulling simulation. The simulation showed that the base pair opening reaction starts with the unbinding of the 5′-base, followed by the unbinding of the 3′-base (i.e., the 5′-base is less stable than the 3′-base). These previous unfolding simulations offered significant insights into the pathways and transition states. However, as shown below, several important issues remain.One intriguing problem is the rate model for base pairing. There are currently three main types of models. In the first type of model, the barrier ΔG+‡ for closing a base pair is dominated by the entropic cost ΔS for positioning the nucleotides to the base-paired configuration and the barrier ΔG−‡ for opening a base pair is the enthalpic cost ΔH for disrupting the hydrogen bonds and base stacking interactions (18, 59, 60). In the second type of model, ΔG+‡ is the net free energy change for base pairing ΔG = ΔH ? TΔS and ΔG−‡ is zero (76, 77). In the third type of model, ΔG±‡=±ΔG/2 is used (78). In addition to the above three main types, other models, such as more sophisticated hybrid rate models, have been proposed (29).In this paper, we report a hybrid method (see Fig. 1) to investigate the single base pairing process. In contrast to the previous simulations for temperature- or force-induced unfolding reactions, we directly model the folding process here (i.e., the base pair closing process). Specifically, we use MD simulations to identify the conformational clusters. Based on the network of the conformational clusters as a reduced conformational ensemble, we apply kinetic Monte Carlo (KMC) and master equation (ME) methods to elucidate the detailed roles of base pairing and stacking interactions, as well as the roles of water and ions (79–82). The study reveals previously unidentified kinetics pathways, misfolded states, and rate-limiting steps. A clear understanding of the microscopic details of the elementary kinetic move is a prerequisite for further rigorous study of large-scale RNA kinetic studies. The method described here may provide a feasible way to develop a rate model for the base pair/stack-based kinetic move set. Furthermore, the mechanism of RNA single base folding may provide useful insights into many biologically significant processes, such as nucleotide flipping (83) in helicases and base pair fraying (84) (as the possible first step for nucleic duplex melting in nucleic acid enzymatic processes).Open in a separate windowFig. 1.(A) Folding of a single nucleotide (G1, red) from the unfolded (Left) to the native folded (Right) state. (B) Exhaustive sampling for the (discrete) conformations of the G1 nucleotide (Right) through enumeration of the torsion angles (formed by the blue bonds). (C) Schematic plot shows the trajectories on the energy landscape (depicted with two reaction coordinates for clarity) explored by the MD simulations. The lines, open circles, and hexagons denote the trajectories; the initial states; and the (centroid structures of the) clusters, respectively. (D) Conformational network based on six clusters. (E) The rmsds to the different clusters provide information about the structural changes in a MD trajectory.     

Decoupling between Shockley partials and stacking faults strengthens multiprincipal element alloys

Mechanical properties are fundamental to structural materials, where dislocations play a decisive role in describing their mechanical behavior. Although the high-yield stresses of multiprincipal element alloys (MPEAs) have received extensive attention in the last decade, the relation between their mechanistic origins remains elusive. Our multiscale study of density functional theory, atomistic simulations, and high-resolution microscopy shows that the excellent mechanical properties of MPEAs have diverse origins. The strengthening effects through Shockley partials and stacking faults can be decoupled in MPEAs, breaking the conventional wisdom that low stacking fault energies are coupled with wide partial dislocations. This study clarifies the mechanistic origins for the strengthening effects, laying the foundation for physics-informed predictive models for materials design.

Multiprincipal element alloys (MPEAs) have triggered ever-increasing interest from the physics and materials science community due to their huge unexplored compositional space and superior physical, mechanical, and functional properties (1–12). They also provide an ideal platform to study fundamental physical mechanisms (6, 9, 13, 14). With the rise of MPEAs, understanding their mechanical properties has become a central topic in materials science in the last decade. In face-centered cubic (fcc) MPEAs, the motion of partial dislocations (Shockley partials) and their associated stacking faults (SF) defines their mechanical properties. Alloys with low SF energies (SFEs) have more extended SFs, which are generally believed to have more strength and ductility through twinning-induced plasticity (TWIP) and transformation-induced plasticity (TRIP) mechanisms (15–17).Although extensive endeavors have been made, the commonalities in the origins of high-yield stresses shared by many MPEAs remain elusive. Among the most common intrinsic contributions of yield stresses are the lattice friction (or Peierls stress) and solute solution strengthening (18–22). Since the birth of MPEAs, it has been a controversy about the relative importance of Peierls stress among the other contributions of yield stress, including the solid-solution strengthening effect (18, 21–23). Many researchers assume small Peierls stresses based on the common wisdom of conventional alloys and pure metals (24, 25) and the low SFEs in MPEAs. Low SFEs usually accompany small Peierls stresses. Overall, this controversy originates from the lack of accurate dislocation geometry in MPEAs, which allows for a direct, critical evaluation of the Peierls stress. There are reports on the dislocation geometry in MPEAs, but almost all of them focused on the widths of SFs (26–28). In contrast, the core widths of Shockley partials are rarely reported for MPEAs, partly due to the difficulty in measurements and partly due to unawareness of its importance. To address this issue, we need very accurate determination of the core width of the Shockley partials. It is an important input parameter for mechanical simulations and various theories and models (21, 29–31). Here, we adopt three of the most extensively studied MPEAs, NiCoCr, VCoNi, and CoCrFeNiMn, and their only common fcc element, Ni, to address the above issues.The commonalities in the origins of high-yield stresses shared by the MPEAs can be indicated by the minimum energy profile along the dislocation motion path, i.e., the increased energies introduced by generalized SFEs (GSFEs; Fig. 1A). The local minima of the curves are SFEs, and the maxima are the theoretical energy barriers for pure shearing, which is a good indicator of the changes of Peierls stresses. Assisted by the accurate density functional theory (DFT), we compute GSFE curves for several representative MPEAs and their common fcc component Ni. This identifies a surprising fact: One of the representative MPEAs, NiCoCr, has a decoupled strengthening effect, i.e., it has a narrower dislocation core of Shockley partial than pure Ni, although its SF is much wider than Ni. Usually, in fcc alloys, when SFE is lower, its unstable SFE (USFE) (maximal GSFE) is also lower, which is coupled. Examples include the two other MPEAs, VCoNi and CoCrFeNiMn, and many Mg alloys (basal plane dislocations) (25) and Al alloys (32). However, NiCoCr does not follow this convention. The understanding from multiscale simulations, atomistic simulations, and the high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) images rationalizes the narrow core of Shockley partials. These results clearly reveal the diverse and decoupled mechanistic origins for the strengthening effects in the MPEAs with excellent mechanical properties.Open in a separate windowFig. 1.GSFEs of three representative MPEAs and pure Ni. (A) The schematic for the generation of GSFs along the slip direction. The displacement 0.75 is equivalent to –0.25 due to the adopted periodic boundary condition. (B) The atom models at two representative displacements for GSFs. (C) The dashed lines are the fitting of the data points to equation γ=γ0sin 2(πx)+(γu−γ0/2)sin 2(2πx) (64, 65). (D) The GSFEs in C are along the path indicated by the white arrows on the gamma surface, i.e., the minimum energy projected along the path denoted by the orange arrow. The GSFE curves reveal the origin for the wide SF and smaller half-width of Shockley partial of NiCoCr than Ni. We need to decrease SFE, while increasing γ_u, in order to optimize the mechanical properties.     

SOD1 aggregation in ALS mice shows simplistic test tube behavior

A longstanding challenge in studies of neurodegenerative disease has been that the pathologic protein aggregates in live tissue are not amenable to structural and kinetic analysis by conventional methods. The situation is put in focus by the current progress in demarcating protein aggregation in vitro, exposing new mechanistic details that are now calling for quantitative in vivo comparison. In this study, we bridge this gap by presenting a direct comparison of the aggregation kinetics of the ALS-associated protein superoxide dismutase 1 (SOD1) in vitro and in transgenic mice. The results based on tissue sampling by quantitative antibody assays show that the SOD1 fibrillation kinetics in vitro mirror with remarkable accuracy the spinal cord aggregate buildup and disease progression in transgenic mice. This similarity between in vitro and in vivo data suggests that, despite the complexity of live tissue, SOD1 aggregation follows robust and simplistic rules, providing new mechanistic insights into the ALS pathology and organism-level manifestation of protein aggregation phenomena in general.So far, the difficulty to experimentally measure protein aggregation in live tissue has focused many researchers to infer mechanistic details of neurodegenerative disease from molecular studies in vitro. An important outcome of this in vitro development is the establishment of rational protocols for quantitative assessment of protein aggregation data (1–4), which now start to consolidate our view of what is happening (5). Protein aggregation follows general and simplistic rules dictated by the amino acid sequence. However, the sheer number of competing aggregation sites within a typical protein chain (6) makes the process intrinsically malleable and dependent on experimental conditions (7). The nagging concern is then to what extent these already complex in vitro data are transferable to the even more complex situation in vivo? Here, we shed light on this question by comparing directly in vitro aggregation kinetics with corresponding data from transgenic mice using a recently developed in vivo quantification strategy based on antibodies (8). Our model system is the aggregation of superoxide dismutase 1 (SOD1) associated with the motor neuron disease ALS (8) (Fig. 1). A key feature of this system is that the immature apoSOD1 monomer, which is also implicated as a precursor in human pathology (9–12), needs to be globally unfolded to fibrillate in vitro (7) (Fig. 1). This simplistic behavior presents the experimental advantage that the fibrillation kinetics of apoSOD1 show simple dependence on structural stability (13, 14):ΔGD−N=−RTlnKD−N=−RTln[N][D],[1]where N is the soluble native structure, and D is the aggregation-competent unfolded state. Accordingly, it has been shown that the in vitro fibrillation of apoSOD1 displays the characteristic fingerprint of fragmentation-assisted growth (15) with a square root dependence on [D] (7), consistent with the requirement of sample agitation to expedite the reaction (1–4, 10). Analogous fibrillation behavior is found for β2-microglobulin (2), yeast prions Sup35 (16) and Ure2p (17), insulin (18), WW domain (19), TI 127 (20), and α-synuclein (21). The main difference between these proteins seems to be that some are intrinsically disordered and constantly aggregation-competent by lacking the ability to hide sticky sequence material by folding. In this study, we see that this simplistic in vitro behavior also translates to the more complex conditions in live tissue: the survival times of ALS mice expressing SOD1 variants of different stabilities are directly correlated with the in vivo levels of globally unfolded protein. Also, spinal cords of mice expressing the human SOD1 mutation G93A show exponential buildup of SOD1 aggregates with a square root dependence on log[D] indistinguishable from the fibrillation kinetics observed in agitated test tubes. The data raise fundamental questions about not only the striking resemblance between mouse and test tube aggregation but also, the apparent differences with human ALS pathology, which seems to have less ordered progression. Clues to the latter, however, are hinted in data from homozygous D90A mice showing two strains of structurally distinct SOD1 aggregates.Open in a separate windowFig. 1.SOD1 aggregation in vitro and in ALS mice. (A) Aggregation of SOD in test tubes yields fibrillar structures similar to those of other proteins (7). (B) Immunohistochemistry of the ventral horn in the terminal hSOD1^G93A mouse showing characteristics of aggresomes (44). (C) Competition between SOD1 folding and fibrillation in vitro, where elongation occurs by unfolded monomers through an encounter complex (7). The question that we ask is how do the in vitro and in vivo aggregations compare mechanistically. (D) Agitation-induced fibrillation in vitro with representative data from an SOD1 mutant in 0 (blue) and 5 M (red) urea with the associated statistics of τ_1/2 for repeated measures. To account for this statistical variation, we use the distribution average (Table S1). (E) Log plot of ν_max vs. τ_1/2 for all individual measures in this study showing uniform behavior of the various SOD1 mutants and a slope of one characteristic for exponential growth (1–4). ALS-associated SOD1 mutations examined in ALS mice (red) (Table S1), other ALS-associated mutations (blue) (Table S1), and SOD1 control mutations (black) (Table S1).     

Structures and topological defects in pressure-driven lyotropic chromonic liquid crystals

Lyotropic chromonic liquid crystals are water-based materials composed of self-assembled cylindrical aggregates. Their behavior under flow is poorly understood, and quantitatively resolving the optical retardance of the flowing liquid crystal has so far been limited by the imaging speed of current polarization-resolved imaging techniques. Here, we employ a single-shot quantitative polarization imaging method, termed polarized shearing interference microscopy, to quantify the spatial distribution and the dynamics of the structures emerging in nematic disodium cromoglycate solutions in a microfluidic channel. We show that pure-twist disclination loops nucleate in the bulk flow over a range of shear rates. These loops are elongated in the flow direction and exhibit a constant aspect ratio that is governed by the nonnegligible splay-bend anisotropy at the loop boundary. The size of the loops is set by the balance between nucleation forces and annihilation forces acting on the disclination. The fluctuations of the pure-twist disclination loops reflect the tumbling character of nematic disodium cromoglycate. Our study, including experiment, simulation, and scaling analysis, provides a comprehensive understanding of the structure and dynamics of pressure-driven lyotropic chromonic liquid crystals and might open new routes for using these materials to control assembly and flow of biological systems or particles in microfluidic devices.

Lyotropic chromonic liquid crystals (LCLCs) are aqueous dispersions of organic disk-like molecules that self-assemble into cylindrical aggregates, which form nematic or columnar liquid crystal phases under appropriate conditions of concentration and temperature (1–6). These materials have gained increasing attention in both fundamental and applied research over the past decade, due to their distinct structural properties and biocompatibility (4, 7–14). Used as a replacement for isotropic fluids in microfluidic devices, nematic LCLCs have been employed to control the behavior of bacteria and colloids (13, 15–20).Nematic liquid crystals form topological defects under flow, which gives rise to complex dynamical structures that have been extensively studied in thermotropic liquid crystals (TLCs) and liquid crystal polymers (LCPs) (21–29). In contrast to lyotropic liquid crystals that are dispersed in a solvent and whose phase can be tuned by either concentration or temperature, TLCs do not need a solvent to possess a liquid-crystalline state and their phase depends only on temperature (30). Most TLCs are shear-aligned nematics, in which the director evolves toward an equilibrium out-of-plane polar angle. Defects nucleate beyond a critical Ericksen number due to the irreconcilable alignment of the directors from surface anchoring and shear alignment in the bulk flow (24, 31–33). With an increase in shear rate, the defect type can transition from π-walls (domain walls that separate regions whose director orientation differs by an angle of π) to ordered disclinations and to a disordered chaotic regime (34). Recent efforts have aimed to tune and control the defect structures by understanding the relation between the selection of topological defect types and the flow field in flowing TLCs. Strategies to do so include tuning the geometry of microfluidic channels, inducing defect nucleation through the introduction of isotropic phases or designing inhomogeneities in the surface anchoring (35–39). LCPs are typically tumbling nematics for which α2α3 < 0, where α2 and α3 are the Leslie viscosities. This leads to a nonzero viscous torque for any orientation of the director, which allows the director to rotate in the shear plane (22, 29, 30, 40). The tumbling character of LCPs facilitates the nucleation of singular topological defects (22, 40). Moreover, the molecular rotational relaxation times of LCPs are longer than those of TLCs, and they can exceed the timescales imposed by the shear rate. As a result, the rheological behavior of LCPs is governed not only by spatial gradients of the director field from the Frank elasticity, but also by changes in the molecular order parameter (25, 41–43). With increasing shear rate, topological defects in LCPs have been shown to transition from disclinations to rolling cells and to worm-like patterns (25, 26, 43).Topological defects occurring in the flow of nematic LCLCs have so far received much more limited attention (44, 45). At rest, LCLCs exhibit unique properties distinct from those of TLCs and LCPs (1, 2, 4–6, 44). In particular, LCLCs have significant elastic anisotropy compared to TLCs; the twist Frank elastic constant, K2, is much smaller than the splay and bend Frank elastic constants, K1 and K3. The resulting relative ease with which twist deformations can occur can lead to a spontaneous symmetry breaking and the emergence of chiral structures in static LCLCs under spatial confinement, despite the achiral nature of the molecules (4, 46–51). When driven out of equilibrium by an imposed flow, the average director field of LCLCs has been reported to align predominantly along the shear direction under strong shear but to reorient to an alignment perpendicular to the shear direction below a critical shear rate (52–54). A recent study has revealed a variety of complex textures that emerge in simple shear flow in the nematic LCLC disodium cromoglycate (DSCG) (44). The tumbling nature of this liquid crystal leads to enhanced sensitivity to shear rate. At shear rates γ˙<1 s−1, the director realigns perpendicular to the flow direction adapting a so-called log-rolling state characteristic of tumbling nematics. For 1 s−1<γ˙<10 s−1, polydomain textures form due to the nucleation of pure-twist disclination loops, for which the rotation vector is parallel to the loop normal, and mixed wedge-twist disclination loops, for which the rotation vector is perpendicular to the loop normal (44, 55). Above γ˙>10 s−1, the disclination loops gradually transform into periodic stripes in which the director aligns predominantly along the flow direction (44).Here, we report on the structure and dynamics of topological defects occurring in the pressure-driven flow of nematic DSCG. A quantitative evaluation of such dynamics has so far remained challenging, in particular for fast flow velocities, due to the slow image acquisition rate of current quantitative polarization-resolved imaging techniques. Quantitative polarization imaging traditionally relies on three commonly used techniques: fluorescence confocal polarization microscopy, polarizing optical microscopy, and LC-Polscope imaging. Fluorescence confocal polarization microscopy can provide accurate maps of birefringence and orientation angle, but the fluorescent labeling may perturb the flow properties (56). Polarizing optical microscopy requires a mechanical rotation of the polarizers and multiple measurements, which severely limits the imaging speed. LC-Polscope, an extension of conventional polarization optical microscopy, utilizes liquid crystal universal compensators to replace the compensator used in conventional polarization microscopes (57). This leads to an enhanced imaging speed and better compensation for polarization artifacts of the optical system. The need for multiple measurements to quantify retardance, however, still limits the acquisition rate of LC-Polscopes.We overcome these challenges by using a single-shot quantitative polarization microscopy technique, termed polarized shearing interference microscopy (PSIM). PSIM combines circular polarization light excitation with off-axis shearing interferometry detection. Using a custom polarization retrieval algorithm, we achieve single-shot mapping of the retardance, which allows us to reach imaging speeds that are limited only by the camera frame rate while preserving a large field-of-view and micrometer spatial resolution. We provide a brief discussion of the optical design of PSIM in Materials and Methods; further details of the measurement accuracy and imaging performance of PSIM are reported in ref. 58.Using a combination of experiments, numerical simulations and scaling analysis, we show that in the pressure-driven flow of nematic DSCG solutions in a microfluidic channel, pure-twist disclination loops emerge for a certain range of shear rates. These loops are elongated in the flow with a fixed aspect ratio. We demonstrate that the disclination loops nucleate at the boundary between regions where the director aligns predominantly along the flow direction close to the channel walls and regions where the director aligns predominantly perpendicular to the flow direction in the center of the channel. The large elastic stresses of the director gradient at the boundary are then released by the formation of disclination loops. We show that both the characteristic size and the fluctuations of the pure-twist disclination loops can be tuned by controlling the flow rate.     

Unusual magnetotransport in twisted bilayer graphene

We present transport measurements of bilayer graphene with a 1.38^∘ interlayer twist. As with other devices with twist angles substantially larger than the magic angle of 1.1^∘, we do not observe correlated insulating states or band reorganization. However, we do observe several highly unusual behaviors in magnetotransport. For a large range of densities around half filling of the moiré bands, magnetoresistance is large and quadratic. Over these same densities, the magnetoresistance minima corresponding to gaps between Landau levels split and bend as a function of density and field. We reproduce the same splitting and bending behavior in a simple tight-binding model of Hofstadter’s butterfly on a triangular lattice with anisotropic hopping terms. These features appear to be a generic class of experimental manifestations of Hofstadter’s butterfly and may provide insight into the emergent states of twisted bilayer graphene.

The mesmerizing Hofstadter butterfly spectrum arises when electrons in a two-dimensional periodic potential are immersed in an out-of-plane magnetic field. When the magnetic flux Φ through a unit cell is a rational multiple p / q of the magnetic flux quantum Φ0=h/e, each Bloch band splits into q subbands (1). The carrier densities corresponding to gaps between these subbands follow straight lines when plotted as a function of normalized density n/ns and magnetic field (2). Here, n_s is the density of carriers required to fill the (possibly degenerate) Bloch band. These lines can be described by the Diophantine equation (n/ns)=t(Φ/Φ0)+s for integers s and t. In experiments, they appear as minima or zeros in longitudinal resistivity coinciding with Hall conductivity quantized at σxy=te2/h (3, 4). Hofstadter originally studied magnetosubbands emerging from a single Bloch band on a square lattice. In the following decades, other authors considered different lattices (5–7), the effect of anisotropy (6, 8–10), next-nearest-neighbor hopping (11–15), interactions (16, 17), density wave states (9), and graphene moirés (18, 19).It took considerable ingenuity to realize clean systems with unit cells large enough to allow conventional superconducting magnets to reach Φ/Φ0∼1. The first successful observation of the butterfly in electrical transport measurements was in GaAs/AlGaAs heterostructures with lithographically defined periodic potentials (20–22). These experiments demonstrated the expected quantized Hall conductance in a few of the largest magnetosubband gaps. In 2013, three groups mapped out the full butterfly spectrum in both density and field in heterostructures based on monolayer (23, 24) and bilayer (25) graphene. In all three cases, the authors made use of the 2% lattice mismatch between their graphene and its encapsulating hexagonal boron nitride (hBN) dielectric. With these layers rotationally aligned, the resulting moiré pattern was large enough in area that gated structures studied in available high-field magnets could simultaneously approach normalized carrier densities and magnetic flux ratios of 1. Later work on hBN-aligned bilayer graphene showed that, likely because of electron–electron interactions, the gaps could also follow lines described by fractional s and t (26).In twisted bilayer graphene (TBG), a slight interlayer rotation creates a similar-scale moiré pattern. Unlike with graphene–hBN moirés, in TBG there is a gap between lowest and neighboring moiré subbands (27). As the twist angle approaches the magic angle of 1.1^∘ the isolated moiré bands become flat (28, 29), and strong correlations lead to fascinating insulating (30–37), superconducting (31–33, 35–37), and magnetic (34, 35, 38) states. The strong correlations tend to cause moiré subbands within a fourfold degenerate manifold to move relative to each other as one tunes the density, leading to Landau levels that project only toward higher magnitude of density from charge neutrality and integer filling factors (37, 39). This correlated behavior obscures the single-particle Hofstadter physics that would otherwise be present.In this work, we present measurements from a TBG device twisted to 1.38^∘. When we apply a perpendicular magnetic field, a complicated and beautiful fan diagram emerges. In a broad range of densities on either side of charge neutrality, the device displays large, quadratic magnetoresistance. Within the magnetoresistance regions, each Landau level associated with ν=±8,±12,±16,… appears to split into a pair, and these pairs follow complicated paths in field and density, very different from those predicted by the usual Diophantine equation. Phenomenology similar in all qualitative respects appears in measurements on several regions of this same device with similar twist angles and in two separate devices, one at 1.59^∘ and the other at 1.70^∘ (see SI Appendix for details).We reproduce the unusual features of the Landau levels (LLs) in a simple tight-binding model on a triangular lattice with anisotropy and a small energetic splitting between two species of fermions. At first glance, this is surprising, because that model does not represent the symmetries of the experimental moiré structure. We speculate that the unusual LL features we experimentally observe can generically emerge from spectra of Hofstadter models that include the same ingredients we added to the triangular lattice model. With further theoretical work it may be possible to use our measurements to gain insight into the underlying Hamiltonian of TBG near the magic angle.     

Sharp thresholds limit the benefit of defector avoidance in cooperation on networks

Consider a cooperation game on a spatial network of habitat patches, where players can relocate between patches if they judge the local conditions to be unfavorable. In time, the relocation events may lead to a homogeneous state where all patches harbor the same relative densities of cooperators and defectors, or they may lead to self-organized patterns, where some patches become safe havens that maintain an elevated cooperator density. Here we analyze the transition between these states mathematically. We show that safe havens form once a certain threshold in connectivity is crossed. This threshold can be analytically linked to the structure of the patch network and specifically to certain network motifs. Surprisingly, a forgiving defector avoidance strategy may be most favorable for cooperators. Our results demonstrate that the analysis of cooperation games in ecological metacommunity models is mathematically tractable and has the potential to link topics such as macroecological patterns, behavioral evolution, and network topology.

Cooperation, behavior that leads to benefits for others at a cost to oneself, is widespread across biological systems, ranging from cells cooperating to form organisms, to cooperation among individuals in populations and among microbiotic and macrobiotic taxa in ecosystems. In many cases the costs of cooperation are high. Hence, how cooperative behavior persists in a population represents a fundamental question in biology (1–8). In general, cooperation is most likely to evolve and persist if there are mechanisms that directly or indirectly benefit cooperators’ reproductive success. Examples include kin selection, punishment of defectors who forgo the cooperative investment, or a direct self-benefit such as in cases of investment into a common good (4).Among the most general mechanisms that can favor cooperation is the notion of network or spatial reciprocity (1, 9–11). In classical examples of reciprocity, cooperation creates favorable conditions for other proximal cooperators (4). A result is the emergence of cooperative havens, where the rewards generated by mutual cooperation have enriched some physical or topological neighborhoods. The formation of cooperative neighborhoods in structured populations, where individuals interact with only a limited subset of the population, has traditionally been studied on networks, where each node represents an individual agent and an edge means that the two connected individuals play against each other (1, 10, 12–19). By assuming weak selection and treating space implicitly, the resulting systems can often be analyzed mathematically. Although this framework has become a powerful tool for conceptual understanding, it represents a strong abstraction from real-world ecology where interactions, and hence cooperative behaviors, occur often randomly within a location that is itself embedded in a larger spatial context (20–23). By focusing on spatially explicit models of cooperation, we gain the opportunity to understand feedbacks between the rules of the game, movement strategies, and long-term persistence of cooperation at larger scales (10, 20, 23–27).Here we study a model of cooperation in spatially structured populations inspired by ecological metacommunities (21–23, 28), where network nodes—instead of individuals—represent habitat patches containing many interacting individuals, and edges mean that two patches are connected by dispersal of those individuals (Fig. 1A). Each patch is a location where games are played, harboring cooperator and defector subpopulations which grow and shrink in time due to internal interactions and movement among locations. Metacommunity models allow one to represent the effects of physical spatial structure directly and explicitly. Moreover, they can be analyzed using master stability functions, which can be used to untangle the impacts of local dynamics and network structure (28–30). We use this ability to explore how different movement strategies impact the outcomes of a cooperation game as a function of network structure.Open in a separate windowFig. 1.Emergence of a heterogeneous stationary state on a two-patch network. (A) Schematic of the spatial game, showing local payoff (Π) relationships among cooperators and defectors occupying the same patch (gray circles) and the dispersal route between them. (B) Difference in equilibrium densities of both types in patches 1 and 2 as link strength is varied. Arrows refer to the example time series shown in C and D. Initial conditions were uniformly drawn from [10−4,10−3], and the patch with the largest initial cooperator density is patch 1. (C) The homogeneous steady state, with the same equilibrium densities of C and D across locations. (Inset) Network showing the proportions of each type in each patch. (D) The same game but with faster diffusion (larger δ), showing emergence of a heterogeneous steady state with higher cooperator densities in patch 1. Parameters are R = 3, S=2, T = 5, P=0.2, μ = 1, and α=3.