Although our existing one-dimensional (1D) model provides a successful quantitative description of rupture events, a 1D description is somewhat limited. We therefore derive a two-dimensional (2D) model which allows us to investigate characteristics of earthquakes (EQs) and slow slip events (SSEs) that are only apparent in a second dimension. We find that the leading edge of an EQ rupture in the direction of the global shear stress (x-direction) is wider in the plane of the crustal fault (y-direction) than the trailing edge. The direction of the slip velocity is primarily in the x-direction. EQ ruptures expand in both the x- and y-directions. In SSEs, the rupture also expands in both directions for a short period of time, then, after pulses are formed, there is no further expandsion, i.e., the pulse shape remains practically unchanged. The 2D simulations show the seismic moment (M) versus time (T) scaling law may be expressed by the relation , where . For EQs, is approximately equal to 3, while for SSEs, tends to 1 as the pulses develop and become fully established. The 2D model quantitatively describes the parameters of a seismic pulse in episodic tremor and slip (ETS) phenomena and also explains important features of non-volcanic tremor associated with this pulse, such as reverse tremor migration. The 2D model confirms the basic findings obtained previously by the 1D model: (i) The type of a seismic event, EQ or SSE, and the type of a rupture, crack-like or pulse-like, are determined solely by the fault strength, the ratio of the shear to normal stress, and the gradient in this ratio; (ii) The shape of a rupture and position of the maximum slip is determined by the spatial distribution of the initial stress.

In vivo Nuclear Magnetic Resonance (NMR) spectroscopy has great potential to interpret the biochemical response of organisms to their environment, thus making it an essential tool in understanding toxic mechanisms. However, magnetic susceptibility distortions lead to 1D NMR spectra of living organisms with lines that are too broad to identify and quantify metabolites, necessitating the use of 2D 1H–13C Heteronuclear Single Quantum Coherence (HSQC) as a primary tool. While quantitative 2D HSQC is well established, to our knowledge it has yet to be applied in vivo. This study represents a simple pilot study that compares two of the most popular quantitative 2D HSQC approaches to determine if quantitative results can be directly obtained in vivo in isotopically enriched Daphnia magna (water flea). The results show the perfect-HSQC experiment performs very well in vivo, but the decoupling scheme used is critical for accurate quantitation. An improved decoupling approach derived using optimal control theory is presented here that improves the accuracy of metabolite concentrations that can be extracted in vivo down to micromolar concentrations. When combined with 2D Electronic Reference To access In vivo Concentrations (ERETIC) protocols, the protocol allows for the direct extraction of in vivo metabolite concentrations without the use of internal standards that can be detrimental to living organisms. Extracting absolute metabolic concentrations in vivo is an important first step and should, for example, be important for the parameterization as well as the validation of metabolic flux models in the future.

We have developed a model that describes the major characteristics of a rupture, ranging from regular earthquakes (EQs) to slow slip events (SSEs), including episodic tremor and slip (ETS). Previous model predictions, while accurate, are based on a highly idealized initial stress distribution and a simple velocity-dependent expression for friction. The full scope of the model has, therefore, not been fully demonstrated. Further developments, presented here, include more physically realistic treatments of both the initial conditions and friction. Model predictions are: (1) The type of a seismic event, i.e. regular EQ or SSE, is determined by the fault strength, the shear to normal stress ratio, and the gradient in the ratio. Quantitative values for these crucial parameters are also obtained here; (2) Rupture velocities for regular EQs range from a fraction of the shear wave velocity up to the supershear velocity. The maximum slip velocity for regular EQs is typically on the order of 1 m/s. For SSEs, the slip velocity ranges widely from a few cm/year up to 0.1 m/s. The magnitude of the stress drop may vary in the range from 1% to 10% of the initial shear stress for regular EQs and from 0.1% to 10% for the SSEs; (3) The rupture can expend as a crack-like mode or as a self-healing pulse mode. The type of rupture mode is determined by the stress ratio and its gradient. If stress heterogeneity has a step-like shape, the rupture is always crack-like. If the heterogeneity is localized, the rupture is crack-like for sufficiently large values of the stress and its gradient and is pulse-like for smaller values; (4) After an instability develops, rupture dynamics do not depend on the relative values of the rate-state a and b parameters, i.e., it does not matter if frictional sliding is velocity-weakening or velocity-strengthening. Models based on the rate-state are not consistent with this result.

The Bloch equation and its variants constitute the fundamental dynamical model for arbitrary two-level systems. Many important processes, including those in more complicated systems, can be modeled and understood through the two-level approximation. It is therefore of widespread relevance, especially as it relates to understanding dissipative processes in current cutting-edge applications of quantum mechanics. Although the Bloch equation has been the subject of considerable analysis in the 70 years since its inception, there is still, perhaps surprisingly, significant work that can be done. This paper extends the scope of previous analyses. It provides a framework for more fully understanding the dynamics of dissipative two-level systems. A solution is derived that is compact, tractable, and completely general, in contrast to previous results. Any solution of the Bloch equation depends on three roots of a cubic polynomial that are crucial to the time dependence of the system. The roots are typically only sketched out qualitatively, with no indication of their dependence on the physical parameters of the problem. Degenerate roots, which modify the solutions, have been ignored altogether. Here the roots are obtained explicitly in terms of a single real-valued root that is expressed as a simple function of the system parameters. For the conventional Bloch equation, a simple graphical representation of this root is presented that makes evident the explicit time dependence of the system for each point in the parameter space. Several intuitive, visual models of system dynamics are developed. A Euclidean coordinate system is identified in which any generalized Bloch equation is separable, i.e., the sum of commuting rotation and relaxation operators. The time evolution in this frame is simply a rotation followed by relaxation at modified rates that play a role similar to the standard longitudinal and transverse rates. These rates are functions of the applied field, which provides information towards control of the dissipative process. The Bloch equation also describes a system of three coupled harmonic oscillators, providing additional perspective on dissipative systems.

In this work, optimal control theory was used to design efficient excitation schemes in highly conductive materials, where both the radio frequency field strength and phase vary as a function of penetration depth. A pulse was designed to achieve phase alignment between signals at different depths within the conductor and thus to obtain higher signals from that region. In addition, an efficient suppression pulse was designed by insuring mutual suppression between the signals from various depths in the sample. The performance of the new approach was demonstrated experimentally for a bulk lithium sample for the excitation problem and for a biphasic metal/liquid sample for the selective suppression pulse.

A binary mixture of superparamagnetic colloidal particles is confined between glass plates such that the large particles become fixed and provide a two-dimensional disorderedmatrix for the still mobile small particles, which form a fluid. By varying fluid and matrix area fractions and tuning the interactions between the superparamagnetic particles via an external magnetic field, different regions of the state diagram are explored. The mobile particles exhibit delocalized dynamics at smallmatrix area fractions and localized motion at highmatrix area fractions, and the localization transition is rounded by the soft interactions [T. O. E. Skinner et al., Phys. Rev. Lett. 111, 128301 ( 2013)]. Expanding on previous work, we find the dynamics of the tracers to be strongly heterogeneous and show that molecular dynamics simulations of an ideal gas confined in a fixed matrix exhibit similar behavior. The simulations show how these soft interactions make the dynamics more heterogeneous compared to the disordered Lorentz gas and lead to strong non-Gaussian fluctuations.

The Bloch equation is the fundamental dynamical model applicable to arbitrary two-level systems. Analytical solutions to date are incomplete for a number of reasons that motivate further investigation. The solution obtained here for the propagator, which generates the time evolution of the system and embodies all the system dynamics, is compact and completely general. The parameter space that results in division by zero in previous treatments is explicitly defined and accommodated in the solution. Polynomial roots required for the solution are expressed in terms of a single real root obtained using simple functional forms. A simple graphical rendition of this root is developed that clarifies and characterizes its dependence on the physical parameters of the problem. As a result, the explicit time dependence of the system as a function of its physical parameters is immediately evident. Several intuitive models of system dynamics are also developed. In particular, the Bloch equation is separable in the proper coordinate system, written as the sum of a relaxation operator and either a null operator or a commuting rotation. The propagator thus drives either pure relaxation or relaxation followed by a rotation. The paper provides a basis for increased physical insight into the Bloch equation and its widespread applications.

The Frenkel-Kontorova (FR) model and its continuum approximation, the sine-Gordon (SG) equation, are widely used to model a variety of important nonlinear physical systems. Many practical applications require the wave-train solution, which includes many solitons. In such cases, an important and relevant extension of these models applies Whitham's averaging procedure to the SG equation. The resulting SG modulation equations describe the behavior of important measurable system parameters that are the average of the small-scale solutions given by the SG equation.A fundamental problem of modern physics that is the topic of this paper is the description of the transitional process from a static to a dynamic frictional regime. We have shown that the SG modulation equations are a suitable apparatus for describing this transition. The model provides relations between kinematic (rupture and slip velocities) and dynamic (shear and normal stresses) parameters of the transition process. A particular advantage of the model is its ability to describe frictional processes over a wide range of rupture and slip velocities covering seismic events ranging from regular earthquakes, with rupture velocities on the order of a few km/s, to slow slip events, with rupture velocities on the order of a few km/day. (C) 2016 Elsevier B.V. All rights reserved.

The dynamics of quantum N-level systems is equivalent to rotation of the state vector in a real, Euclidean vector space. The relevance of this result to general magnetic resonance systems is considered here, with a focus on heteronuclear decoupling. An exact general theory for spin-1/2 IS systems is presented. This article provides a vehicle for more comprehensive details on the derivation of the theoretical model, insights into its application, new visual representations, and computational efficiencies for simulations. The theory provides a unifying mechanism for decoupling and has clarified possibilities for developing improved decoupling sequences, unconstrained by former limitations of periodic and cyclic sequences.

Description of the transitional process from a static to a dynamic frictional regime is a fundamental problem of modern physics. Previously, we developed a model based on the well-known Frenkel-Kontorova model to describe dry macroscopic friction. Here, this model has been modified to include the effect of dissipation in derived relations between the kinematic and dynamic parameters of a transition process. The main (somewhat counterintuitive) result is a demonstration that the rupture (i.e., detachment front) velocity of the slip pulse which arises during the transition does not depend on friction. The only parameter (besides the elastic and plastic properties of the medium) controlling the rupture velocity is the shear to normal stress ratio. In contrast to the rupture velocity, the slip velocity does depend on friction. The model we have developed describes these processes over a wide range of rupture and slip velocities (up to 7 orders of magnitude) allowing, in particular, the consideration of seismic events ranging from regular earthquakes, with rupture velocities on the order of a few km/s, to slow slip events, with rupture velocities of a few km/day.

Accurate prediction of the saturation dependence of different modes of transport in porous media, such as those due to conductivity, air permeability, and diffusion, is of broad interest in engineering and natural resources management. Most current predictions use a "bundle of capillary tubes" concept, which, despite its widespread use, is a severely distorted idealization of natural porous media. In contrast, percolation theory provides a reliable and powerful means to model interconnectivity of disordered networks and porous materials. In this study, we invoke scaling concepts from percolation theory and effective medium theory to predict the saturation dependence of modes of transport - hydraulic and electrical conductivity, air permeability, and gas diffusion - in two disturbed soils. Universal scaling from percolation theory predicts the saturation dependence of air permeability and gas diffusion accurately, even when the percolation threshold for airflow is estimated from the porosity. We also find that the non-universal scaling obtained from the critical path analysis (CPA) of percolation theory can make excellent predictions of hydraulic and electrical conductivity under partially saturated conditions.

Transport in porous media is quite complex, and still yields occasional surprises. In geological porous media, the rate at which chemical reactions (e.g., weathering and dissolution) occur is found to diminish by orders of magnitude with increasing time or distance. The temporal rates of laboratory experiments and field observations differ, and extrapolating from laboratory experiments (in months) to field rates (in millions of years) can lead to order-of-magnitude errors. The reactions are transport-limited, but characterizing them using standard solute transport expressions can yield results in agreement with experiment only if spurious assumptions and parameters are introduced. We previously developed a theory of non-reactive solute transport based on applying critical path analysis to the cluster statistics of percolation. The fractal structure of the clusters can be used to generate solute distributions in both time and space. Solute velocities calculated from the temporal evolution of that distribution have the same time dependence as reaction-rate scaling in a wide range of field studies and laboratory experiments, covering some 10 decades in time. The present theory thus both explains a wide range of experiments, and also predicts changes in the scaling behavior in individual systems with increasing time and/or length scales. No other theory captures these variations in scaling by invoking a single physical mechanism. Because the successfully predicted chemical reactions include known results for silicate weathering rates, our theory provides a framework for understanding changes in the global carbon cycle, including its effects on extinctions, climate change, soil production, and denudation rates. It further provides a basis for understanding the fundamental time scales of hydrology and shallow geochemistry, as well as the basis of industrial agriculture.

An exact general theory of heteronuclear decoupling is presented for spin-1/2 IS systems. RF irradiation applied to the I spins both modifies and generates additional couplings between states of the system. The recently derived equivalence between the dynamics of any N-level quantum system and a system of classical coupled harmonic oscillators makes explicit the exact physical couplings between states. Decoupling is thus more properly viewed as a complex intercoupling. The sign of antiphase magnetization plays a fundamental role in decoupling. A one-to-one correspondence is demonstrated between ±2SyIz and the sense of the S-spin coupling evolution. Magnetization Sx is refocused to obtain the desired decoupled state when ∫2SyIzdt=0. The exact instantaneous coupling at any time during the decoupling sequence is readily obtained in terms of the system states, showing that the creation of two-spin coherence is crucial for reducing the effective scalar coupling, as required for refocusing to occur. Representative examples from new aperiodic sequences as well as standard cyclic, periodic composite-pulse and adiabatic decoupling sequences illustrate the decoupling mechanism. The more general aperiodic sequences, obtained using optimal control, realize the potential inherent in the theory for significantly improved decoupling.

Ultra-high-field NMR spectroscopy requires an increased bandwidth for heteronuclear decoupling, especially in biomolecular NMR applications. Composite pulse decoupling cannot provide sufficient bandwidth at practical power levels, and adiabatic pulse decoupling with sufficient bandwidth is compromised by sideband artifacts. A novel low-power, broadband heteronuclear decoupling pulse is presented that generates minimal, ultra-low sidebands. The pulse was derived using optimal control theory and represents a new generation of decoupling pulses free from the constraints of periodic and cyclic sequences. In comparison to currently available state-of-the-art methods this novel pulse provides greatly improved decoupling performance that satisfies the demands of high-field biomolecular NMR spectroscopy.

The saturated hydraulic conductivity, Ks, is a fundamental characteristic of subsurface flow and the hydrologic cycle. However, direct measurement of Ks is time consuming. Recently, air permeability measurements at 50 and 100 cm of H2O tension, a relatively dry condition in coarse and medium‐textured soils, have been used to estimate Ks. In this short communication, we present a theoretical framework for relating Ks to ka(ϕ), the air permeability as a function of porosity, ϕ, under completely dry conditions. We used power‐law scaling from continuum percolation theory to develop a theoretical relationship between Ks and ka(ϕ). Our result, similar in form to a published logarithmic equation with empirical coefficients, gives a physical interpretation to these empirical coefficients.

The porosity dependence of the formation factor for geologic media is examined from the perspective of universal scaling laws from percolation and effective medium theories. Over much of the range of observed porosity, the expected percolation scaling is observed, but the values of the numerical prefactor do not conform to the simple predictions from percolation theory. Combining effective medium and percolation theories produces a numerical prefactor whose value depends on both the threshold porosity and the porosity above which the formation factor crosses from percolation to effective medium scaling. This change allows extraction of a numerical value of the prefactor, which is reasonably close to experimental values. Subsequent evaluation of the porosity dependence of the formation factor shows that difficulties in prior comparisons of theory and experiment are largely removed when percolation scaling is allowed to transition to effective medium scaling far above the percolation threshold.

Distance and spins: A new broadband, frequency-selective inversion pulse has been devised for pulsed electron–electron double resonance (PELDOR) spectroscopy, a unique method used to determine distance distributions between two or more paramagnetic centers 2–10 nm apart. The PELDOR modulation depth can be increased by substituting the rectangular pump pulse by a shaped pulse. As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re-organized for online delivery, but are not copy-edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors. Please note: The publisher is not responsible for the content or functionality of any supporting information supplied by the authors. Any queries (other than missing content) should be directed to the corresponding author for the article.

The dynamics of states representing arbitrary N-level quantum systems, including dissipative systems, can be modelled exactly by the dynamics of classical coupled oscillators. There is a direct one-to-one correspondence between the quantum states and the positions of the oscillators. Quantum coherence, expectation values, and measurement probabilities for system observables can therefore be realized from the corresponding classical states. The time evolution of an N-level system is represented as the rotation of a real state vector in hyperspace, as previously known for density matrix states but generalized here to Schrodinger states. A single rotor in n dimensions is then mapped directly to n oscillators in one physical dimension. The number of oscillators needed to represent N-level systems scales linearly with N for Schrodinger states, in contrast to N^2 for the density matrix formalism. Although the well-known equivalence (SU(2), SO(3) homomorphism) of 2-level quantum dynamics to a rotation in real, physical space cannot be generalized to arbitrary N-level systems, representing quantum dynamics by a system of coupled harmonic oscillators in one physical dimension is general for any N. Values for the classical coupling constants are readily obtained from the system Hamiltonian, allowing construction of classical mechanical systems that can provide visual insight into the dynamics of abstract quantum systems as well as a metric for characterizing the interface between quantum and classical mechanics.

Tortuosity is a property of porous media that is invoked and used in the literature on hydrology, soil science, physics, and engineering. It has been defined in a variety of ways, one of which is a purely geometrical concept. In this study, we focused on the geometrical tortuosity and developed a model based on percolation theory and the finite-size scaling approach. Our result, developed for porous media of any saturation, expresses the tortuosity as a power-law function of the water content, the critical water content, and the system size. The model parameters include the fractal dimension of either the backbone (the flow-carrying part of porous media) or the optimal path (for crossing the system between two opposite faces), which have clear physical meaning. The results may be combined with power laws, which percolation theory provides for the hydraulic or electrical conductivity, to develop the appropriate form that expresses the connectivity as a function of the water content of a pore space. Comparison with numerical simulations and experiments reported in the literature indicated that the model estimates the tortuosity accurately. Our model and results leads us to conclude that pore connectivity and tortuosity should be treated as two distinct properties and that significant uncertainty in results for the saturation dependence of the hydraulic conductivity is traceable to the confusion that mixes the two as a single concept.