Critical Collapse in Einstein-Gauss-Bonnet Gravity in Five and Six Dimensions

Einstein-Gauss-Bonnet gravity (EGB) provides a natural higher dimensional and higher order curvature generalization of Einstein gravity. It contains a new, presumably microscopic, length scale that should affect short distance properties of the dynamics, such as Choptuik scaling. We present the results of a numerical analysis in generalized flat slice co-ordinates of self-gravitating massless scalar spherical collapse in five and six dimensional EGB gravity near the threshold of black hole formation. Remarkably, the behaviour is universal (i.e. independent of initial data) but qualitatively different in five and six dimensions. In five dimensions there is a minimum horizon radius, suggestive of a first order transition between black hole and dispersive initial data. In six dimensions no radius gap is evident. Instead, below the GB scale there is a change in the critical exponent and echoing period.Comment: 21 pages, 39 figures, a couple of references and two new figures adde

Huge First-Order Metamagnetic Transition in the Paramagnetic Heavy-Fermion System CeTiGe

We report on the observation of large, step-like anomalies in the magnetization ($\Delta M = 0.74$\,$\mu_{\rm B}$/Ce), in the magnetostriction ($\Delta l/l_{0} = 2.0 \cdot 10^{-3}$), and in the magnetoresistance in polycrystals of the paramagnetic heavy-fermion system CeTiGe at a critical magnetic field $\mu_0 H_c \approx$ 12.5\,T at low temperatures. The size of these anomalies is much larger than those reported for the prototypical heavy-fermion metamagnet CeRu$_2$Si$_2$. Furthermore, hysteresis between increasing and decreasing field data indicate a real thermodynamic, first-order type of phase transition, in contrast to the crossover reported for CeRu$_2$Si$_2$. Analysis of the resistivity data shows a pronounced decrease of the electronic quasiparticle mass across $H_c$. These results establish CeTiGe as a new metamagnetic Kondo-lattice system, with an exceptionally large, metamagnetic transition of first-order type at a moderate field.Comment: 5 pages, 4 figure

Neuromorphic Twins for Networked Control and Decision-Making

We consider the problem of remotely tracking the state of and unstable linear time-invariant plant by means of data transmitted through a noisy communication channel from an algorithmic point of view. Assuming the dynamics of the plant are known, does there exist an algorithm that accepts a description of the channel's characteristics as input, and returns 'Yes' if the transmission capabilities permit the remote tracking of the plant's state, 'No' otherwise? Does there exist an algorithm that, in case of a positive answer, computes a suitable encoder/decoder-pair for the channel? Questions of this kind are becoming increasingly important with regards to future communication technologies that aim to solve control engineering tasks in a distributed manner. In particular, they play an essential role in digital twinning, an emerging information processing approach originally considered in the context of Industry 4.0. Yet, the abovementioned questions have been answered in the negative with respect to algorithms that can be implemented on idealized digital hardware, i.e., Turing machines. In this article, we investigate the remote state estimation problem in view of the Blum-Shub-Smale computability framework. In the broadest sense, the latter can be interpreted as a model for idealized analog computation. Especially in the context of neuromorphic computing, analog hardware has experienced a revival in the past view years. Hence, the contribution of this work may serve as a motivation for a theory of neuromorphic twins as a counterpart to digital twins for analog hardware

Improving Gravitational Wave Extraction Using Weyl Characteristic Fields

We present a detailed methodology for extracting the full set of Newman-Penrose Weyl scalars from numerically generated spacetimes without requiring a tetrad that is completely orthonormal or perfectly aligned to the principal null directions. We also describe an extrapolation technique for computing the Weyl scalars' contribution at asymptotic null infinity in post-processing. These methods have continued to be used to produce $\Psi_4$ and $h$ waveforms for the Simulating eXtreme Spacetimes Waveform Catalog and now have been expanded to produce the entire set of Weyl scalars. These new waveform quantities are critical for the future of gravitational wave astronomy in order to understand the finite-amplitude gauge differences that can occur in numerical waveforms. We also present a new analysis of the accuracy of waveforms produced by the Spectral Einstein Code. While ultimately we expect Cauchy Characteristic Extraction to yield more accurate waveforms, the extraction techniques described here are far easier to implement and have already proven to be a viable way to produce production-level waveforms that can meet the demands of current gravitational-wave detectors

Quantum Criticality in doped CePd_1-xRh_x Ferromagnet

CePd_1-xRh_x alloys exhibit a continuous evolution from ferromagnetism (T_C= 6.5 K) at x = 0 to a mixed valence (MV) state at x = 1. We have performed a detailed investigation on the suppression of the ferromagnetic (F) phase in this alloy using dc-(\chi_dc) and ac-susceptibility (\chi_ac), specific heat (C_m), resistivity (\rho) and thermal expansion (\beta) techniques. Our results show a continuous decrease of T_C (x) with negative curvature down to T_C = 3K at x*= 0.65, where a positive curvature takes over. Beyond x*, a cusp in cac is traced down to T_C* = 25 mK at x = 0.87, locating the critical concentration between x = 0.87 and 0.90. The quantum criticality of this region is recognized by the -log(T/T_0) dependence of C_m/T, which transforms into a T^-q (~0.5) one at x = 0.87. At high temperature, this system shows the onset of valence instability revealed by a deviation from Vegard's law (at x_V~0.75) and increasing hybridization effects on high temperature \chi_dc and \rho. Coincidentally, a Fermi liquid contribution to the specific heat arises from the MV component, which becomes dominant at the CeRh limit. In contrast to antiferromagnetic systems, no C_m/T flattening is observed for x > x_cr rather the mentioned power law divergence, which coincides with a change of sign of \beta. The coexistence of F and MV components and the sudden changes in the T dependencies are discussed in the context of randomly distributed magnetic and Kondo couplings.Comment: 11 pages, 11 figure

Low temperature magnetic phase diagram of the cubic non-Fermi liquid system CeIn_(3-x)Sn_x

In this paper we report a comprehensive study of the magnetic susceptibility (\chi), resistivity (\rho), and specific heat (C_P), down to 0.5 K of the cubic CeIn_(3-x)Sn_x alloy. The ground state of this system evolves from antiferromagnetic (AF) in CeIn_3(T_N=10.2 K) to intermediate-valent in CeSn_3, and represents the first example of a Ce-lattice cubic non-Fermi liquid (NFL) system where T_N(x) can be traced down to T=0 over more than a decade of temperature. Our results indicate that the disappearance of the AF state occurs near x_c ~ 0.7, although already at x ~ 0.4 significant modifications of the magnetic ground state are observed. Between these concentrations, clear NFL signatures are observed, such as \rho(T)\approx \rho_0 + A T^n (with n<1.5) and C_P(T)\propto -T ln(T) dependencies. Within the ordered phase a first order phase transition occurs for 0.25 < x < 0.5. With larger Sn doping, different weak \rho(T) dependencies are observed at low temperatures between x=1 and x=3 while C_P/T shows only a weak temperature dependence.Comment: 7 pages, 7 figures. Accepted in Eur. J. Phys.

Fixing the BMS Frame of Numerical Relativity Waveforms with BMS Charges

The Bondi-van der Burg-Metzner-Sachs (BMS) group, which uniquely describes the symmetries of asymptotic infinity and therefore of the gravitational waves that propagate there, has become increasingly important for accurate modeling of waveforms. In particular, waveform models, such as post-Newtonian (PN) expressions, numerical relativity (NR), and black hole perturbation theory, produce results that are in different BMS frames. Consequently, to build a model for the waveforms produced during the merging of compact objects, which ideally would be a hybridization of PN, NR, and black hole perturbation theory, one needs a fast and robust method for fixing the BMS freedoms. In this work, we present the first means of fixing the entire BMS freedom of NR waveforms to match the frame of either PN waveforms or black hole perturbation theory. We achieve this by finding the BMS transformations that change certain charges in a prescribed way -- e.g., finding the center-of-mass transformation that maps the center-of-mass charge to a mean of zero. We find that this new method is 20 times faster, and more correct when mapping to the superrest frame, than previous methods that relied on optimization algorithms. Furthermore, in the course of developing this charge-based frame fixing method, we compute the PN expression for the Moreschi supermomentum to 3PN order without spins and 2PN order with spins. This Moreschi supermomentum is effectively equivalent to the energy flux or the null memory contribution at future null infinity $\mathscr{I}^{+}$. From this PN calculation, we also compute oscillatory ($m\not=0$ modes) and spin-dependent memory terms that have not been identified previously or have been missing from strain expressions in the post-Newtonian literature. <br

Fixing the BMS Frame of Numerical Relativity Waveforms

Understanding the Bondi-Metzner-Sachs (BMS) frame of the gravitational waves produced by numerical relativity is crucial for ensuring that analyses on such waveforms are performed properly. It is also important that models are built from waveforms in the same BMS frame. Up until now, however, the BMS frame of numerical waveforms has not been thoroughly examined, largely because the necessary tools have not existed. In this paper, we show how to analyze and map to a suitable BMS frame for numerical waveforms calculated with the Spectral Einstein Code (SpEC). However, the methods and tools that we present are general and can be applied to any numerical waveforms. We present an extensive study of 13 binary black hole systems that broadly span parameter space. From these simulations, we extract the strain and also the Weyl scalars using both SpECTRE's Cauchy-characteristic extraction module and also the standard extrapolation procedure with a displacement memory correction applied during post-processing. First, we show that the current center-of-mass correction used to map these waveforms to the center-of-mass frame is not as effective as previously thought. Consequently, we also develop an improved correction that utilizes asymptotic Poincar\'e charges instead of a Newtonian center-of-mass trajectory. Next, we map our waveforms to the post-Newtonian (PN) BMS frame using a PN strain waveform. This helps us find the unique BMS transformation that minimizes the $L^{2}$ norm of the difference between the numerical and PN strain waveforms during the early inspiral phase. We find that once the waveforms are mapped to the PN BMS frame, they can be hybridized with a PN strain waveform much more effectively than if one used any of the previous alignment schemes, which only utilize the Poincar\'e transformations

High-accuracy numerical models of Brownian thermal noise in thin mirror coatings

Brownian coating thermal noise in detector test masses is limiting the sensitivity of current gravitational-wave detectors on Earth. Therefore, accurate numerical models can inform the ongoing effort to minimize Brownian coating thermal noise in current and future gravitational-wave detectors. Such numerical models typically require significant computational resources and time, and often involve closed-source commercial codes. In contrast, open-source codes give complete visibility and control of the simulated physics and enable direct assessment of the numerical accuracy. In this article, we use the open-source SpECTRE numerical-relativity code and adopt a novel discontinuous Galerkin numerical method to model Brownian coating thermal noise. We demonstrate that SpECTRE achieves significantly higher accuracy than a previous approach at a fraction of the computational cost. Furthermore, we numerically model Brownian coating thermal noise in multiple sub-wavelength crystalline coating layers for the first time. Our new numerical method has the potential to enable fast exploration of realistic mirror configurations, and hence to guide the search for optimal mirror geometries, beam shapes and coating materials for gravitational-wave detectors