BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Paul A. Johnson

Paul A. Johnson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.geo-phphysics.chem-phcond-mat.str-elcond-mat.dis-nncond-mat.mtrl-scicond-mat.softcond-mat.stat-mechmath-phmath.MPnlin.PSphysics.comp-phphysics.data-an

Catalog footprint

What is connected

15works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2022arXiv

Bivariational Principle for an Antisymmetrized Product of Nonorthogonal Geminals Appropriate for Strong Electron Correlation

We develop a bivariational principle for an antisymmetric product of nonorthogonal geminals. Special cases reduce to the antisymmetric product of strongly-orthogonal geminals (APSG), the generalized valence bond-perfect pairing (GVB-PP), and the antisymmetrized geminal power (AGP) wavefunctions. The presented method employs wavefunctions of the same type as Richardson-Gaudin (RG) states, but which are not eigenvectors of a model Hamiltonian which would allow for more freedom in the mean-field. The general idea is to work with the same state in a primal picture in terms of pairs, and in a dual picture in terms of pair-holes. This leads to an asymmetric energy expression which may be optimized bivariationally, and is strictly variational when the two representations are consistent. The general approach may be useful in other contexts, such as for computationally feasible variational coupled-cluster methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Density Matrices of Seniority-Zero Geminal Wavefunctions

Scalar products and density matrix elements of closed-shell pair geminal wavefunctions are evaluated directly in terms of the pair amplitudes, resulting in an analogue of Wick's theorem for fermions or bosons. This expression is in general intractable, but it is shown how it becomes feasible in three distinct ways for Richardson-Gaudin (RG) states, the antisymmetrized geminal power, and the antisymmetrized product of strongly-orthogonal geminals. Dissociation curves for hydrogen chains are computed with off-shell RG states and the antisymmetrized product of interacting geminals. Both are near exact suggesting that the incorrect results observed with ground state RG states are fixable using a different RG state.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Near-exact treatment of seniority-zero ground and excited states with a Richardson-Gaudin mean-field

Eigenvectors of the reduced Bardeen-Cooper-Schrieffer Hamiltonian, Richardson-Gaudin (RG) states, are used as a variational wavefunction Ansatz for strongly-correlated electronic systems. These states are geminal products whose coefficients are solutions of non-linear equations. Previous results showed un-physical behaviour but in this contribution it is shown that with only the variational solution for the ground state, all the seniority-zero states are quite well approximated. The difficulty is in choosing the correct RG state. The systems studied showed a clear choice and we expect it should always be possible to reason physically which state to choose.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Shear-wave-induced softening and simultaneous compaction in dense granular media through acoustic lubrication at flow heterogeneities

We report the simultaneous softening and compaction of a confined dense granular pack in acoustic resonance experiments. Elastic softening is manifested by a reduction of the shear-wave speed, as the wave amplitude increases beyond some threshold. No macroscopic rearrangement of grains or dilatancy is observed; instead, elastic softening is accompanied by a tiny amount of compaction on the scale of grain asperities. We explain these apparent contradictory observations using a theoretical model, based on shear transformation zones (STZs), of soft spots and slipping contacts. It predicts a linear shear stress-strain response with negligible macro-plastic deformation due to the small-amplitude acoustic oscillation. However, these waves reduce the interparticle friction and contact stiffness through the acoustic lubrication of grain contacts, resulting in an increase in the structural disorder or compactivity and softening of dynamic modulus. The compaction associated with this microscopic friction decrease is consistent with the prediction by an Ising-like correlation between STZs in the subyield regime.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The temporal limits of predicting fault failure

Machine learning models using seismic emissions can predict instantaneous fault characteristics such as displacement in laboratory experiments and slow slip in Earth. Here, we address whether the acoustic emission (AE) from laboratory experiments contains information about near-future frictional behavior. The approach uses a convolutional encoder-decoder containing a transformer layer. We use as input progressively larger AE input time windows and progressively larger output friction time windows. The attention map from the transformer is used to interpret which regions of the AE contain hidden information corresponding to future frictional behavior. We find that very near-term predictive information is indeed contained in the AE signal, but farther into the future the predictions are progressively worse. Notably, information for predicting near future frictional failure and recovery are found to be contained in the AE signal. This first effort predicting future fault frictional behavior with machine learning will guide efforts for applications in Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Predicting Fault Slip via Transfer Learning

Data-driven machine-learning for predicting instantaneous and future fault-slip in laboratory experiments has recently progressed markedly due to large training data sets. In Earth however, earthquake interevent times range from 10's-100's of years and geophysical data typically exist for only a portion of an earthquake cycle. Sparse data presents a serious challenge to training machine learning models. Here we describe a transfer learning approach using numerical simulations to train a convolutional encoder-decoder that predicts fault-slip behavior in laboratory experiments. The model learns a mapping between acoustic emission histories and fault-slip from numerical simulations, and generalizes to produce accurate results using laboratory data. Notably slip-predictions markedly improve using the simulation-data trained-model and training the latent space using a portion of a single laboratory earthquake-cycle. The transfer learning results elucidate the potential of using models trained on numerical simulations and fine-tuned with small geophysical data sets for potential applications to faults in Earth.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Autonomous Extraction of Millimeter-scale Deformation in InSAR Time Series Using Deep Learning

Systematic characterization of slip behaviours on active faults is key to unraveling the physics of tectonic faulting and the interplay between slow and fast earthquakes. Interferometric Synthetic Aperture Radar (InSAR), by enabling measurement of ground deformation at a global scale every few days, may hold the key to those interactions. However, atmospheric propagation delays often exceed ground deformation of interest despite state-of-the art processing, and thus InSAR analysis requires expert interpretation and a priori knowledge of fault systems, precluding global investigations of deformation dynamics. Here we show that a deep auto-encoder architecture tailored to untangle ground deformation from noise in InSAR time series autonomously extracts deformation signals, without prior knowledge of a fault's location or slip behaviour. Applied to InSAR data over the North Anatolian Fault, our method reaches 2 mm detection, revealing a slow earthquake twice as extensive as previously recognized. We further explore the generalization of our approach to inflation/deflation-induced deformation, applying the same methodology to the geothermal field of Coso, California.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reduced density matrices of Richardson-Gaudin states in the Gaudin algebra basis

Eigenvectors of the reduced Bardeen-Cooper-Schrieffer Hamiltonian have recently been employed as a variational wavefunction ansatz in quantum chemistry. This wavefunction is a mean-field of pairs of electrons (geminals). In this contribution we report optimal expressions for their reduced density matrices in both the original physical basis and the basis of the Richardson-Gaudin pairs. Physical basis expressions were originally reported by Gorohovsky and Bettelheim. In each case, the expressions scale like $\mathcal{O}(N^4)$, with the most expensive step the solution of linear equations. Analytic gradients are also reported in the physical basis. These expressions are an important step towards practical mean-field methods to treat strongly-correlated electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Earthquake Arrival Association with Backprojection and Graph Theory

The association of seismic wave arrivals with causative earthquakes becomes progressively more challenging as arrival detection methods become more sensitive, and particularly when earthquake rates are high. For instance, seismic waves arriving across a monitoring network from several sources may overlap in time, false arrivals may be detected, and some arrivals may be of unknown phase (e.g., P- or S-waves). We propose an automated method to associate arrivals with earthquake sources and obtain source locations applicable to such situations. To do so we use a pattern detection metric based on the principle of backprojection to reveal candidate sources, followed by graph-theory-based clustering and an integer linear optimization routine to associate arrivals with the minimum number of sources necessary to explain the data. This method solves for all sources and phase assignments simultaneously, rather than in a sequential greedy procedure as is common in other association routines. We demonstrate our method on both synthetic and real data from the Integrated Plate Boundary Observatory Chile (IPOC) seismic network of northern Chile. For the synthetic tests we report results for cases with varying complexity, including rates of 500 earthquakes/day and 500 false arrivals/station/day, for which we measure true positive detection accuracy of > 95%. For the real data we develop a new catalog between January 1, 2010 - December 31, 2017 containing 817,548 earthquakes, with detection rates on average 279 earthquakes/day, and a magnitude-of-completion of ~M1.8. A subset of detections are identified as sources related to quarry and industrial site activity, and we also detect thousands of foreshocks and aftershocks of the April 1, 2014 Mw 8.2 Iquique earthquake. During the highest rates of aftershock activity, > 600 earthquakes/day are detected in the vicinity of the Iquique earthquake rupture zone.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly to simulate a sheared granular fault system. We focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements measured at locations on the upper and lower plates just adjacent to the gouge in the direction parallel to the shear direction (x direction). The simulations show that at slips the plate velocities are proportional to the normal load and may be inversely proportional to the square root of the plate's Young's modulus; whereas the driving shear velocity does not show distinct influence on the plate velocities. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than the half of the driving shear velocity, and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory scale fault gouge system and may help reveal the complexities of earthquake frictional dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Probing slow earthquakes with deep learning

Slow earthquakes may trigger failure on neighboring locked faults that are stressed enough to break, and slow slip patterns may evolve before a nearby great earthquake. However, even in the clearest cases such as Cascadia, slow earthquakes and associated tremor have only been observed in intermittent and discrete bursts. By training a convolutional neural network to detect known tremor on a single seismic station in Cascadia, we isolate and identify tremor and slip preceding and following known larger slow events. The deep neural network can be used for the detection of quasi-continuous tremor, providing a proxy that quantifies the slow slip rate. Furthermore, the model trained in Cascadia recognizes tremor in other subduction zones and also along the San Andreas Fault at Parkfield, suggesting a universality of waveform characteristics and source processes, as posited from experiments and theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Acoustically-induced slip in sheared granular layers: application to dynamic earthquake triggering

A fundamental mystery in earthquake physics is ``how can an earthquake be triggered by distant seismic sources?'' Here, we use discrete element method simulations of a granular layer, during stick-slip, that is subject to transient vibrational excitation to gain further insight into the physics of dynamic earthquake triggering. Using Coulomb friction law for grains interaction, we observe delayed triggering of slip in the granular gouge. We find that at a critical vibrational amplitude (strain) there is an abrupt transition from negligible time-advanced slip (clock advance) to full clock advance, {\it i.e.}, transient vibration and triggered slip are simultaneous. The critical strain is order of $10^{-6}$, similar to observations in the laboratory and in Earth. The transition is related to frictional weakening of the granular layer due to a dramatic decrease in coordination number and the weakening of the contact force network. Associated with this frictional weakening is a pronounced decrease in the elastic modulus of the layer. The study has important implications for mechanisms of triggered earthquakes and induced seismic events and points out the underlying processes in response of the fault gouge to dynamic transient stresses.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Projected seniority-two orbital optimization of the Antisymmetric Product of one-reference orbital Geminal

We present a new, non-variational orbital-optimization scheme for the Antisymmetric Product of one-reference orbital Geminal wave function. Our approach is motivated by the observation that an orbital-optimized seniority-zero configuration interaction (CI) expansion yields similar results to an orbital-optimized seniority-zero-plus-two CI expansion [J. Chem. Phys., 135, 044119 (2011)]. A numerical analysis is performed for the C$_2$, LiF and CH$_2$ molecules as well as for the symmetric stretching of hypothetical (linear) hydrogen chains. For these test cases, the proposed orbital-optimization protocol yields similar results to its variational orbital optimization counterpart, but prevents symmetry-breaking of molecular orbitals in most cases.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Are megaquakes clustered?

We study statistical properties of the number of large earthquakes over the past century. We analyze the cumulative distribution of the number of earthquakes with magnitude larger than threshold M in time interval T, and quantify the statistical significance of these results by simulating a large number of synthetic random catalogs. We find that in general, the earthquake record cannot be distinguished from a process that is random in time. This conclusion holds whether aftershocks are removed or not, except at magnitudes below M = 7.3. At long time intervals (T = 2-5 years), we find that statistically significant clustering is present in the catalog for lower magnitude thresholds (M = 7-7.2). However, this clustering is due to a large number of earthquakes on record in the early part of the 20th century, when magnitudes are less certain.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Nonequilibrium and Nonlinear Dynamics in Geomaterials I : The Low Strain Regime

Members of a wide class of geomaterials are known to display complex and fascinating nonlinear and nonequilibrium dynamical behaviors over a wide range of bulk strains, down to surprisingly low values, e.g., 10^{-7}. In this paper we investigate two sandstones, Berea and Fontainebleau, and characterize their behavior under the influence of very small external forces via carefully controlled resonant bar experiments. By reducing environmental effects due to temperature and humidity variations, we are able to systematically and reproducibly study dynamical behavior at strains as low as 10^{-9}. Our study establishes the existence of two strain thresholds, the first, epsilon_L, below which the material is essentially linear, and the second, epsilon_M, below which the material is nonlinear but where quasiequilibrium thermodynamics still applies as evidenced by the success of Landau theory and a simple macroscopic description based on the Duffing oscillator. At strains above epsilon_M the behavior becomes truly nonequilibrium -- as demonstrated by the existence of material conditioning -- and Landau theory no longer applies. The main focus of this paper is the study of the region below the second threshold, but we also comment on how our work clarifies and resolves previous experimental conflicts, as well as suggest new directions of research.

0 citations0 reviewsNo code linkedOpen access