Showing 1–28 of 28 results for author: Tromp, J

Cosserat elasticity as the weak-field limit of Einstein--Cartan relativity

Authors: Matthew Maitra, Jeroen Tromp

Abstract: The weak-field limit of Einstein--Cartan (EC) relativity is studied. The equations of EC theory are rewritten such that they formally resemble those of Einstein General Relativity (EGR); this allows ideas from post-Newtonian theory to be imported without essential change. The equations of motion are then written both at first post-Newtonian (1PN) order and at 1.5PN order. EC theory's 1PN equations… ▽ More The weak-field limit of Einstein--Cartan (EC) relativity is studied. The equations of EC theory are rewritten such that they formally resemble those of Einstein General Relativity (EGR); this allows ideas from post-Newtonian theory to be imported without essential change. The equations of motion are then written both at first post-Newtonian (1PN) order and at 1.5PN order. EC theory's 1PN equations of motion are found to be those of a micropolar/Cosserat elastic medium, along with a decoupled evolution equation for non-classical, spin-related fields. It seems that a necessary condition for these results to hold is that one chooses the non-classical fields to scale with the speed of light in a certain empirically reasonable way. Finally, the 1.5PN equations give greater insight into the coupling between energy-momentum and spin within slowly moving, weakly gravitating matter. Specifically, the weakly relativistic modifications to Cosserat theory involve a gravitational torque and an augmentation of the gravitational force due to a `dynamic mass moment density' with an accompanying `dynamic mass moment density flux', and new forms of linear momentum density captured by a `dynamic mass density flux' and a `dynamic momentum density'. △ Less

Submitted 20 May, 2024; originally announced May 2024.

Minuscule corrections to near-surface solar internal rotation using mode-coupling

Authors: Srijan Bharati Das, Samarth G. Kashyap, Deniz Oktay, Shravan M. Hanasoge, Jeroen Tromp

Abstract: The observed solar oscillation spectrum is influenced by internal perturbations such as flows and structural asphericities. These features induce splitting of characteristic frequencies and distort the resonant-mode eigenfunctions. Global axisymmertric flow -- differential rotation -- is a very prominent perturbation. Tightly constrained rotation profiles as a function of latitude and radius are p… ▽ More The observed solar oscillation spectrum is influenced by internal perturbations such as flows and structural asphericities. These features induce splitting of characteristic frequencies and distort the resonant-mode eigenfunctions. Global axisymmertric flow -- differential rotation -- is a very prominent perturbation. Tightly constrained rotation profiles as a function of latitude and radius are products of established helioseismic pipelines that use observed Dopplergrams to generate frequency-splitting measurements at high precision. However, the inference of rotation using frequency-splittings do not consider the effect of mode-coupling. This approximation worsens for high-angular-degree modes, as they become increasingly proximal in frequency. Since modes with high angular degrees probe the near-surface layers of the Sun, inversions considering coupled modes could potentially lead to more accurate estimates of rotation very close to the surface. In order to investigate if this is indeed the case, we perform inversions for solar differential rotation, considering coupling of modes for angular degrees $160 \leq \ell \leq 300$ in the surface gravity $f$-branch and first-overtone $p$ modes. In keeping with the character of mode coupling, we carry out a non-linear inversion using an eigenvalue solver. Differences in inverted profiles for frequency splitting measurements from MDI and HMI are compared and discussed. We find that corrections to the near-surface differential rotation profile, when accounting for mode-coupling effects, are smaller than 0.003 nHz and hence are insignificant. These minuscule corrections are found to be correlated with the solar cycle. We also present corrections to even-order splitting coefficients, which could consequently impact inversions for structure and magnetic fields. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Comments: 21 pages; 15 figures; Accepted for publication in the Astrophysical Journal Supplement Series

Full-waveform tomography reveals iron spin crossover in Earth lower mantle

Authors: Laura Cobden, Jingyi Zhuang, Wenjie Lei, Renata Wentzcovitch, Jeannot Trampert, Jeroen Tromp

Abstract: Joint interpretation of bulk and shear wave speeds constrains the chemistry of the deep mantle. At all depths, the diversity of wave speeds cannot be explained by an isochemical mantle. Between 1000 and 2500 km depth, hypothetical mantle models containing an electronic spin crossover in (Mg,Fe)O provide a significantly better fit to the wave-speed distributions, as well as more realistic temperatu… ▽ More Joint interpretation of bulk and shear wave speeds constrains the chemistry of the deep mantle. At all depths, the diversity of wave speeds cannot be explained by an isochemical mantle. Between 1000 and 2500 km depth, hypothetical mantle models containing an electronic spin crossover in (Mg,Fe)O provide a significantly better fit to the wave-speed distributions, as well as more realistic temperatures and silica contents, than models without a spin crossover. Below 2500 km, wave speed distributions are explained by enrichment in silica towards the core-mantle-boundary. This silica enrichment may represent the fractionated remains of an ancient basal magma ocean. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Ab initio calculations of third-order elastic coefficients

Authors: Chenxing Luo, Jeroen Tromp, Renata M. Wentzcovitch

Abstract: Third-order elasticity (TOE) theory is predictive of strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the… ▽ More Third-order elasticity (TOE) theory is predictive of strain-induced changes in second-order elastic coefficients (SOECs) and can model elastic wave propagation in stressed media. Although third-order elastic tensors have been determined based on first principles in previous studies, their current definition is based on an expansion of thermodynamic energy in terms of the Lagrangian strain near the natural, or zero pressure, reference state. This definition is inconvenient for predictions of SOECs under significant initial stresses. Therefore, when TOE theory is necessary to study the strain dependence of elasticity, the seismological community has resorted to an empirical version of the theory. This study reviews the thermodynamic definition of the third-order elastic tensor and proposes using an "effective" third-order elastic tensor. An explicit expression for the effective third-order elastic tensor is given and verified. We extend the ab initio approach to calculate third-order elastic tensors under finite pressure and apply it to two cubic systems, namely, NaCl and MgO. As applications and validations, we evaluate (a) strain-induced changes in SOECs and (b) pressure derivatives of SOECs based on ab initio calculations. Good agreement between third-order elasticity-based predictions and numerically calculated values confirms the validity of our theory. △ Less

Submitted 13 November, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Journal ref: Phys.Rev.B 106 (2022) 214104

Inferring solar differential rotation through normal-mode coupling using Bayesian statistics

Authors: Samarth G. Kashyap, Srijan Bharati Das, Shravan M. Hanasoge, Martin F. Woodard, Jeroen Tromp

Abstract: Normal-mode helioseismic data analysis uses observed solar oscillation spectra to infer perturbations in the solar interior due to global and local-scale flows and structural asphericity. Differential rotation, the dominant global-scale axisymmetric perturbation, has been tightly constrained primarily using measurements of frequency splittings via "$a$-coefficients". However, the frequency-splitti… ▽ More Normal-mode helioseismic data analysis uses observed solar oscillation spectra to infer perturbations in the solar interior due to global and local-scale flows and structural asphericity. Differential rotation, the dominant global-scale axisymmetric perturbation, has been tightly constrained primarily using measurements of frequency splittings via "$a$-coefficients". However, the frequency-splitting formalism invokes the approximation that multiplets are isolated. This assumption is inaccurate for modes at high angular degrees. Analysing eigenfunction corrections, which respect cross coupling of modes across multiplets, is a more accurate approach. However, applying standard inversion techniques using these cross-spectral measurements yields $a$-coefficients with a significantly wider spread than the well-constrained results from frequency splittings. In this study, we apply Bayesian statistics to infer $a$-coefficients due to differential rotation from cross spectra for both $f$-modes and $p$-modes. We demonstrate that this technique works reasonably well for modes with angular degrees $\ell=50-291$. The inferred $a_3-$coefficients are found to be within $1$ nHz of the frequency splitting values for $\ell > 200$. We also show that the technique fails at $\ell < 50$ owing to the insensitivity of the measurement to the perturbation. These results serve to further establish mode coupling as an important helioseismic technique with which to infer internal structure and dynamics, both axisymmetric (e.g., meridional circulation) and non-axisymmetric perturbations. △ Less

Submitted 21 January, 2021; originally announced January 2021.

Comments: 20 pages, 11 figures

Journal ref: ApJS 253 (2021) 47

Preconditioned BFGS-based Uncertainty Quantification in elastic Full Waveform Inversion

Authors: Qiancheng Liu, Stephen Beller, Wenjie Lei, Daniel Peter, Jeroen Tromp

Abstract: Full Waveform Inversion (FWI) plays a vital role in reconstructing geophysical structures. The Uncertainty Quantification regarding the inversion results is equally important but has been missing out in most of the current geophysical inversions. Mathematically, uncertainty quantification is involved with the inverse Hessian (or the posterior covariance matrix), which is prohibitive in computation… ▽ More Full Waveform Inversion (FWI) plays a vital role in reconstructing geophysical structures. The Uncertainty Quantification regarding the inversion results is equally important but has been missing out in most of the current geophysical inversions. Mathematically, uncertainty quantification is involved with the inverse Hessian (or the posterior covariance matrix), which is prohibitive in computation and storage for practical geophysical FWI problems. L-BFGS populates as the most efficient Gauss-Newton method; however, in this study, we empower it with the new possibility of accessing the inverse Hessian for uncertainty quantification in FWI. To facilitate the inverse-Hessian retrieval, we put together BFGS (essentially, full-history L-BFGS) with randomized singular value decomposition towards a low-rank approximation of the Hessian inverse. That the rank number equals the number of iterations makes this solution efficient and memory-affordable even for large-scale inversions. Also, based on the adjoint method, we formulate different diagonal Hessian initials as preconditioners and compare their performances in elastic FWI. We highlight our methods with the elastic Marmousi benchmark, demonstrating the applicability of preconditioned BFGS in large-scale FWI and uncertainty quantification. △ Less

Submitted 15 August, 2021; v1 submitted 26 September, 2020; originally announced September 2020.

Sensitivity kernels for inferring Lorentz stresses from normal-mode frequency splittings in the Sun

Authors: Srijan Bharati Das, Tuneer Chakraborty, Shravan M. Hanasoge, Jeroen Tromp

Abstract: Departures from standard spherically symmetric solar models, in the form of perturbations such as global and local-scale flows and structural asphericities, result in the splitting of eigenfrequencies in the observed spectrum of solar oscillations. Drawing from prevalent ideas in normal-mode coupling theory in geophysical literature, we devise a procedure that enables the computation of sensitivit… ▽ More Departures from standard spherically symmetric solar models, in the form of perturbations such as global and local-scale flows and structural asphericities, result in the splitting of eigenfrequencies in the observed spectrum of solar oscillations. Drawing from prevalent ideas in normal-mode coupling theory in geophysical literature, we devise a procedure that enables the computation of sensitivity kernels for general Lorentz stress fields in the Sun. Mode coupling due to any perturbation requires careful consideration of self- and cross-coupling of multiplets. Invoking the isolated-multiplet approximation allows for limiting the treatment to purely self-coupling, requiring significantly less computational resources. We identify the presence of such isolated multiplets under the effect of Lorentz stresses in the Sun. Currently, solar missions allow precise measurements of self-coupling of multiplets via "$a$-coefficients" and the cross-spectral correlation signal which enables the estimation of the "structure coefficients". We demonstrate the forward problem for both self-coupling ($a$-coefficients) and cross-coupling (structure coefficients). In doing so, we plot the self-coupling kernels and estimate $a$-coefficients arising from a combination of deep-toroidal and surface-dipolar axisymmetric fields. We also compute the structure coefficients for an arbitrary general magnetic field (real and solenoidal) and plot the corresponding "splitting function", a convenient way to visualize the splitting of multiplets under 3D internal perturbations. The results discussed in this paper pave the way to formally pose an inverse problem, and infer solar internal magnetic fields. △ Less

Submitted 29 April, 2020; originally announced April 2020.

Comments: 32 pages, 8 figures, Astrophysical Journal

Harnessing the Power of Many: Extensible Toolkit for Scalable Ensemble Applications

Authors: Vivek Balasubramanian, Matteo Turilli, Weiming Hu, Matthieu Lefebvre, Wenjie Lei, Guido Cervone, Jeroen Tromp, Shantenu Jha

Abstract: Many scientific problems require multiple distinct computational tasks to be executed in order to achieve a desired solution. We introduce the Ensemble Toolkit (EnTK) to address the challenges of scale, diversity and reliability they pose. We describe the design and implementation of EnTK, characterize its performance and integrate it with two distinct exemplar use cases: seismic inversion and ada… ▽ More Many scientific problems require multiple distinct computational tasks to be executed in order to achieve a desired solution. We introduce the Ensemble Toolkit (EnTK) to address the challenges of scale, diversity and reliability they pose. We describe the design and implementation of EnTK, characterize its performance and integrate it with two distinct exemplar use cases: seismic inversion and adaptive analog ensembles. We perform nine experiments, characterizing EnTK overheads, strong and weak scalability, and the performance of two use case implementations, at scale and on production infrastructures. We show how EnTK meets the following general requirements: (i) implementing dedicated abstractions to support the description and execution of ensemble applications; (ii) support for execution on heterogeneous computing infrastructures; (iii) efficient scalability up to O(10^4) tasks; and (iv) fault tolerance. We discuss novel computational capabilities that EnTK enables and the scientific advantages arising thereof. We propose EnTK as an important addition to the suite of tools in support of production scientific computing. △ Less

Submitted 16 May, 2018; v1 submitted 23 October, 2017; originally announced October 2017.

A spectral-infinite-element solution of Poisson's equation: an application to self gravity

Authors: Hom Nath Gharti, Jeroen Tromp

Abstract: We solve Poisson's equation by combining a spectral-element method with a mapped infinite-element method. We focus on problems in geostatics and geodynamics, where Earth's gravitational field is determined by Poisson's equation inside the Earth and Laplace's equation in the rest of space. Spectral elements are used to capture the internal field, and infinite elements are used to represent the exte… ▽ More We solve Poisson's equation by combining a spectral-element method with a mapped infinite-element method. We focus on problems in geostatics and geodynamics, where Earth's gravitational field is determined by Poisson's equation inside the Earth and Laplace's equation in the rest of space. Spectral elements are used to capture the internal field, and infinite elements are used to represent the external field. To solve the weak form of Poisson/Laplace equation, we use Gauss-Legendre-Lobatto quadrature in spectral elements inside the domain of interest. Outside the domain, we use Gauss-Radau quadrature in the infinite direction, and Gauss-Legendre-Lobatto quadrature in the other directions. We illustrate the efficiency and accuracy of the method by comparing the gravitational fields of a homogeneous sphere and the Preliminary Reference Earth Model (PREM) with (semi-)analytical solutions. △ Less

Submitted 2 June, 2017; originally announced June 2017.

Comments: 15 pages, 14 figures

Double-difference adjoint seismic tomography

Authors: Yanhua O. Yuan, Frederik J. Simons, Jeroen Tromp

Abstract: We introduce a `double-difference' method for the inversion for seismic wavespeed structure based on adjoint tomography. Differences between seismic observations and model predictions at individual stations may arise from factors other than structural heterogeneity, such as errors in the assumed source-time function, inaccurate timings, and systematic uncertainties. To alleviate the corresponding… ▽ More We introduce a `double-difference' method for the inversion for seismic wavespeed structure based on adjoint tomography. Differences between seismic observations and model predictions at individual stations may arise from factors other than structural heterogeneity, such as errors in the assumed source-time function, inaccurate timings, and systematic uncertainties. To alleviate the corresponding nonuniqueness in the inverse problem, we construct differential measurements between stations, thereby reducing the influence of the source signature and systematic errors. We minimize the discrepancy between observations and simulations in terms of the differential measurements made on station pairs. We show how to implement the double-difference concept in adjoint tomography, both theoretically and in practice. We compare the sensitivities of absolute and differential measurements. The former provide absolute information on structure along the ray paths between stations and sources, whereas the latter explain relative (and thus higher-resolution) structural variations in areas close to the stations. Whereas in conventional tomography a measurement made on a single earthquake-station pair provides very limited structural information, in double-difference tomography one earthquake can actually resolve significant details of the structure. The double-difference methodology can be incorporated into the usual adjoint tomography workflow by simply pairing up all conventional measurements; the computational cost of the necessary adjoint simulations is largely unaffected. Rather than adding to the computational burden, the inversion of double-difference measurements merely modifies the construction of the adjoint sources for data assimilation. △ Less

Submitted 6 July, 2016; originally announced July 2016.

Comments: 21 pages, 17 figures, accepted for publication by the Geophysical Journal International

Anelastic sensitivity kernels with parsimonious storage for adjoint tomography and full waveform inversion

Authors: Dimitri Komatitsch, Zhinan Xie, Ebru Bozdag, Elliott Sales de Andrade, Daniel Peter, Qinya Liu, Jeroen Tromp

Abstract: We introduce a technique to compute exact anelastic sensitivity kernels in the time domain using parsimonious disk storage. The method is based on a reordering of the time loop of time-domain forward/adjoint wave propagation solvers combined with the use of a memory buffer. It avoids instabilities that occur when time-reversing dissipative wave propagation simulations. The total number of required… ▽ More We introduce a technique to compute exact anelastic sensitivity kernels in the time domain using parsimonious disk storage. The method is based on a reordering of the time loop of time-domain forward/adjoint wave propagation solvers combined with the use of a memory buffer. It avoids instabilities that occur when time-reversing dissipative wave propagation simulations. The total number of required time steps is unchanged compared to usual acoustic or elastic approaches. The cost is reduced by a factor of 4/3 compared to the case in which anelasticity is partially accounted for by accommodating the effects of physical dispersion. We validate our technique by performing a test in which we compare the $K_αあるふぁ$ sensitivity kernel to the exact kernel obtained by saving the entire forward calculation. This benchmark confirms that our approach is also exact. We illustrate the importance of including full attenuation in the calculation of sensitivity kernels by showing significant differences with physical-dispersion-only kernels. △ Less

Submitted 30 May, 2016; v1 submitted 19 April, 2016; originally announced April 2016.

Full Waveform Inversion for Time-Distance Helioseismology

Authors: Shravan Hanasoge, Jeroen Tromp

Abstract: Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the so… ▽ More Inferring interior properties of the Sun from photospheric measurements of the seismic wavefield constitutes the helioseismic inverse problem. Deviations in seismic measurements (such as wave travel times) from their fiducial values estimated for a given model of the solar interior imply that the model is inaccurate. Contemporary inversions in local helioseismology assume that properties of the solar interior are linearly related to measured travel-time deviations. It is widely known, however, that this assumption is invalid for sunspots and active regions, and likely for supergranular flows as well. Here, we introduce nonlinear optimization, executed iteratively, as a means of inverting for the sub-surface structure of large-amplitude perturbations. Defining the penalty functional as the $L_2$ norm of wave travel-time deviations, we compute the the total misfit gradient of this functional with respect to the relevant model parameters %(only sound speed in this case) at each iteration around the corresponding model. The model is successively improved using either steepest descent, conjugate gradient, or quasi-Newton limited-memory BFGS. Performing nonlinear iterations requires privileging pixels (such as those in the near-field of the scatterer), a practice not compliant with the standard assumption of translational invariance. Measurements for these inversions, although similar in principle to those used in time-distance helioseismology, require some retooling. For the sake of simplicity in illustrating the method, we consider a 2-D inverse problem with only a sound-speed perturbation. △ Less

Submitted 29 January, 2014; originally announced January 2014.

Comments: 24 pages, 10 figures, to appear in ApJ

Seismic probes of solar interior magnetic structure

Authors: Shravan Hanasoge, Aaron Birch, Laurent Gizon, Jeroen Tromp

Abstract: Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of c… ▽ More Sunspots are prominent manifestations of solar magnetoconvection and imaging their subsurface structure is an outstanding problem of wide physical importance. Travel times of seismic waves that propagate through these structures are typically used as inputs to inversions. Despite the presence of strongly anisotropic magnetic waveguides, these measurements have always been interpreted in terms of changes to isotropic wavespeeds and flow-advection related Doppler shifts. Here, we employ PDE-constrained optimization to determine the appropriate parameterization of the structural properties of the magnetic interior. Seven different wavespeeds fully characterize helioseismic wave propagation: the isotropic sound speed, a Doppler-shifting flow-advection velocity and an anisotropic magnetic velocity. The structure of magnetic media is sensed by magnetoacoustic slow and fast modes and Alfvén waves, each of which propagates at a different wavespeed. We show that even in the case of weak magnetic fields, significant errors may be incurred if these anisotropies are not accounted for in inversions. Translation invariance is demonstrably lost. These developments render plausible the accurate seismic imaging of magnetoconvection in the Sun. △ Less

Submitted 18 July, 2012; originally announced July 2012.

Comments: 4 pages, 4 figures, accepted Physical Review Letters

Detectable seismic consequences of the interaction of a primordial black hole with Earth

Authors: Yang Luo, Shravan Hanasoge, Jeroen Tromp, Frans Pretorius

Abstract: Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here… ▽ More Galaxies observed today are likely to have evolved from density perturbations in the early universe. Perturbations that exceeded some critical threshold are conjectured to have undergone gravitational collapse to form primordial black holes (PBHs) at a range of masses. Such PBHs serve as candidates for cold dark matter and their detection would shed light on conditions in the early universe. Here we propose a mechanism to search for transits of PBHs through/nearby Earth by studying the associated seismic waves. Using a spectral-element method, we simulate and visualize this seismic wave field in Earth's interior. We predict the emergence of two unique signatures, namely, a wave that would arrive almost simultaneously everywhere on Earth's free surface and the excitation of unusual spheroidal modes with a characteristic frequency-spacing in free oscillation spectra. These qualitative characteristics are unaffected by the speed or proximity of the PBH trajectory. The seismic energy deposited by a proximal ${M^{PBH} = 10^{15}}$ g PBH is comparable to a magnitude $M_w=4$ earthquake. The non-seismic collateral damage due to the actual impact of such small PBHs with Earth would be negligible. Unfortunately, the expected collision rate is very low even if PBHs constituted all of dark matter, at ${\sim 10^{-7} {yr}^{-1}}$, and since the rate scales as ${1/M^{PBH}}$, fortunately encounters with larger, Earth-threatening PBHs are exceedingly unlikely. However, the rate at which non-colliding close encounters of PBHs could be detected by seismic activity alone is roughly two orders of magnitude larger --- that is once every hundred thousand years --- than the direct collision rate. △ Less

Submitted 16 March, 2012; originally announced March 2012.

Comments: Accepted, ApJ; 22 pages, 12 figures

The Adjoint Method Applied to Time-Distance Helioseismology

Authors: Shravan Hanasoge, Aaron Birch, Laurent Gizon, Jeroen Tromp

Abstract: For a given {\it misfit function}, a specified optimality measure of a model, its gradient describes the manner in which one may alter properties of the system to march towards a stationary point. The adjoint method, arising from partial-differential-equation-constrained optimization, describes a means of extracting derivatives of a misfit function with respect to model parameters through finite c… ▽ More For a given {\it misfit function}, a specified optimality measure of a model, its gradient describes the manner in which one may alter properties of the system to march towards a stationary point. The adjoint method, arising from partial-differential-equation-constrained optimization, describes a means of extracting derivatives of a misfit function with respect to model parameters through finite computation. It relies on the accurate calculation of wavefields that are driven by two types of sources, namely the average wave-excitation spectrum, resulting in the {\it forward wavefield}, and differences between predictions and observations, resulting in an {\it adjoint wavefield}. All sensitivity kernels relevant to a given measurement emerge directly from the evaluation of an interaction integral involving these wavefields. The technique facilitates computation of sensitivity kernels (Fréchet derivatives) relative to three-dimensional heterogeneous background models, thereby paving the way for non-linear iterative inversions. An algorithm to perform such inversions using as many observations as desired is discussed. △ Less

Submitted 21 May, 2011; originally announced May 2011.

Comments: 23 pages, 9 Figures, accepted by ApJ

The Taurus Spitzer Survey: New Candidate Taurus Members Selected Using Sensitive Mid-Infrared Photometry

Authors: L. M. Rebull, D. L. Padgett, C. -E. McCabe, L. A. Hillenbrand, K. R. Stapelfeldt, A. Noriega-Crespo, S. J. Carey, T. Brooke, T. Huard, S. Terebey, M. Audard, J. -L. Monin, M. Fukagawa, M. Guedel, G. R. Knapp, F. Menard, L. E. Allen, J. R. Angione, C. Baldovin-Saavedra, J. Bouvier, K. Briggs, C. Dougados, N. J. Evans, N. Flagey, S. Guieu , et al. (9 additional authors not shown)

Abstract: We report on the properties of pre-main-sequence objects in the Taurus molecular clouds as observed in 7 mid- and far-infrared bands with the Spitzer Space Telescope. There are 215 previously-identified members of the Taurus star-forming region in our ~44 square degree map; these members exhibit a range of Spitzer colors that we take to define young stars still surrounded by circumstellar dust (… ▽ More We report on the properties of pre-main-sequence objects in the Taurus molecular clouds as observed in 7 mid- and far-infrared bands with the Spitzer Space Telescope. There are 215 previously-identified members of the Taurus star-forming region in our ~44 square degree map; these members exhibit a range of Spitzer colors that we take to define young stars still surrounded by circumstellar dust (noting that ~20% of the bonafide Taurus members exhibit no detectable dust excesses). We looked for new objects in the survey field with similar Spitzer properties, aided by extensive optical, X-ray, and ultraviolet imaging, and found 148 candidate new members of Taurus. We have obtained follow-up spectroscopy for about half the candidate sample, thus far confirming 34 new members, 3 probable new members, and 10 possible new members, an increase of 15-20% in Taurus members. Of the objects for which we have spectroscopy, 7 are now confirmed extragalactic objects, and one is a background Be star. The remaining 93 candidate objects await additional analysis and/or data to be confirmed or rejected as Taurus members. Most of the new members are Class II M stars and are located along the same cloud filaments as the previously-identified Taurus members. Among non-members with Spitzer colors similar to young, dusty stars are evolved Be stars, planetary nebulae, carbon stars, galaxies, and AGN. △ Less

Submitted 16 November, 2009; originally announced November 2009.

Comments: Accepted to ApJS. Two large online-only figures available with the preprint here: http://web.ipac.caltech.edu/staff/rebull/research.html

Journal ref: Astrophys.J.Suppl.186:259-307,2010

Prefix reversals on binary and ternary strings

Authors: Cor Hurkens, Leo van Iersel, Judith Keijsper, Steven Kelk, Leen Stougie, John Tromp

Abstract: Given a permutation pi, the application of prefix reversal f^(i) to pi reverses the order of the first i elements of pi. The problem of Sorting By Prefix Reversals (also known as pancake flipping), made famous by Gates and Papadimitriou (Bounds for sorting by prefix reversal, Discrete Mathematics 27, pp. 47-57), asks for the minimum number of prefix reversals required to sort the elements of a g… ▽ More Given a permutation pi, the application of prefix reversal f^(i) to pi reverses the order of the first i elements of pi. The problem of Sorting By Prefix Reversals (also known as pancake flipping), made famous by Gates and Papadimitriou (Bounds for sorting by prefix reversal, Discrete Mathematics 27, pp. 47-57), asks for the minimum number of prefix reversals required to sort the elements of a given permutation. In this paper we study a variant of this problem where the prefix reversals act not on permutations but on strings over a fixed size alphabet. We determine the minimum number of prefix reversals required to sort binary and ternary strings, with polynomial-time algorithms for these sorting problems as a result; demonstrate that computing the minimum prefix reversal distance between two binary strings is NP-hard; give an exact expression for the prefix reversal diameter of binary strings, and give bounds on the prefix reversal diameter of ternary strings. We also consider a weaker form of sorting called grouping (of identical symbols) and give polynomial-time algorithms for optimally grouping binary and ternary strings. A number of intriguing open problems are also discussed. △ Less

Submitted 4 July, 2006; v1 submitted 21 February, 2006; originally announced February 2006.

Comments: 20 pages. Just submitted for journal publication. Previous version should be considered no longer valid

On the Complexity of the Single Individual SNP Haplotyping Problem

Authors: Rudi Cilibrasi, Leo van Iersel, Steven Kelk, John Tromp

Abstract: We present several new results pertaining to haplotyping. These results concern the combinatorial problem of reconstructing haplotypes from incomplete and/or imperfectly sequenced haplotype fragments. We consider the complexity of the problems Minimum Error Correction (MEC) and Longest Haplotype Reconstruction (LHR) for different restrictions on the input data. Specifically, we look at the gaple… ▽ More We present several new results pertaining to haplotyping. These results concern the combinatorial problem of reconstructing haplotypes from incomplete and/or imperfectly sequenced haplotype fragments. We consider the complexity of the problems Minimum Error Correction (MEC) and Longest Haplotype Reconstruction (LHR) for different restrictions on the input data. Specifically, we look at the gapless case, where every row of the input corresponds to a gapless haplotype-fragment, and the 1-gap case, where at most one gap per fragment is allowed. We prove that MEC is APX-hard in the 1-gap case and still NP-hard in the gapless case. In addition, we question earlier claims that MEC is NP-hard even when the input matrix is restricted to being completely binary. Concerning LHR, we show that this problem is NP-hard and APX-hard in the 1-gap case (and thus also in the general case), but is polynomial time solvable in the gapless case. △ Less

Submitted 6 September, 2005; v1 submitted 11 August, 2005; originally announced August 2005.

Comments: 26 pages. Related to the WABI2005 paper, "On the Complexity of Several Haplotyping Problems", but with more/different results. This papers has just been submitted to the IEEE/ACM Transactions on Computational Biology and Bioinformatics and we are awaiting a decision on acceptance. It differs from the mid-August version of this paper because here we prove that 1-gap LHR is APX-hard. (In the earlier version of the paper we could prove only that it was NP-hard.)

On the Complexity of Several Haplotyping Problems

Authors: Rudi Cilibrasi, Leo van Iersel, Steven Kelk, John Tromp

Abstract: In this paper we present a collection of results pertaining to haplotyping. The first set of results concerns the combinatorial problem of reconstructing haplotypes from incomplete and/or imperfectly sequenced haplotype data. More specifically, we show that an interesting, restricted case of Minimum Error Correction (MEC) is NP-hard, point out problems in earlier claims about a related problem,… ▽ More In this paper we present a collection of results pertaining to haplotyping. The first set of results concerns the combinatorial problem of reconstructing haplotypes from incomplete and/or imperfectly sequenced haplotype data. More specifically, we show that an interesting, restricted case of Minimum Error Correction (MEC) is NP-hard, point out problems in earlier claims about a related problem, and present a polynomial-time algorithm for the ungapped case of Longest Haplotype Reconstruction (LHR). Secondly, we present a polynomial time algorithm for the problem of resolving genotype data using as few haplotypes as possible (the Pure Parsimony Haplotyping Problem, PPH) where each genotype has at most two ambiguous positions, thus solving an open problem posed by Lancia et al in "Haplotyping Populations by Pure Parsimony: Complexity of Exact and Approximation Algorithms." △ Less

Submitted 13 May, 2005; originally announced May 2005.

Comments: 12 pages. Submitted to WABI2005, are awaiting a decision on whether it has been accepted

Limits of Rush Hour Logic Complexity

Authors: John Tromp, Rudi Cilibrasi

Abstract: Rush Hour Logic was introduced in [Flake&Baum99] as a model of computation inspired by the ``Rush Hour'' toy puzzle, in which cars can move horizontally or vertically within a parking lot. The authors show how the model supports polynomial space computation, using certain car configurations as building blocks to construct boolean circuits for a cpu and memory. They consider the use of cars of le… ▽ More Rush Hour Logic was introduced in [Flake&Baum99] as a model of computation inspired by the ``Rush Hour'' toy puzzle, in which cars can move horizontally or vertically within a parking lot. The authors show how the model supports polynomial space computation, using certain car configurations as building blocks to construct boolean circuits for a cpu and memory. They consider the use of cars of length 3 crucial to their construction, and conjecture that cars of size 2 only, which we'll call `Size 2 Rush Hour', do not support polynomial space computation. We settle this conjecture by showing that the required building blocks are constructible in Size 2 Rush Hour. Furthermore, we consider Unit Rush Hour, which was hitherto believed to be trivial, show its relation to maze puzzles, and provide empirical support for its hardness. △ Less

Submitted 15 February, 2005; originally announced February 2005.

ACM Class: F.1.3; F.2

Sharpening Occam's Razor

Authors: Ming Li, John Tromp, Paul Vitanyi

Abstract: We provide a new representation-independent formulation of Occam's razor theorem, based on Kolmogorov complexity. This new formulation allows us to: (i) Obtain better sample complexity than both length-based and VC-based versions of Occam's razor theorem, in many applications. (ii) Achieve a sharper reverse of Occam's razor theorem than previous work. Specifically, we weaken the assumption… ▽ More We provide a new representation-independent formulation of Occam's razor theorem, based on Kolmogorov complexity. This new formulation allows us to: (i) Obtain better sample complexity than both length-based and VC-based versions of Occam's razor theorem, in many applications. (ii) Achieve a sharper reverse of Occam's razor theorem than previous work. Specifically, we weaken the assumptions made in an earlier publication, and extend the reverse to superpolynomial running times. △ Less

Submitted 10 October, 2002; v1 submitted 8 January, 2002; originally announced January 2002.

Comments: LaTeX 13 pages; Proc 8th COCOON, LNCS 2387, Springer-Verlag, Berlin, 2002, 411--419

Report number: CWI Manuscript 1994 ACM Class: F.2, E.4, I.2

Kolmogorov Random Graphs and the Incompressibility Method

Authors: Harry Buhrman, Ming Li, John Tromp, Paul Vitanyi

Abstract: We investigate topological, combinatorial, statistical, and enumeration properties of finite graphs with high Kolmogorov complexity (almost all graphs) using the novel incompressibility method. Example results are: (i) the mean and variance of the number of (possibly overlapping) ordered labeled subgraphs of a labeled graph as a function of its randomness deficiency (how far it falls short of th… ▽ More We investigate topological, combinatorial, statistical, and enumeration properties of finite graphs with high Kolmogorov complexity (almost all graphs) using the novel incompressibility method. Example results are: (i) the mean and variance of the number of (possibly overlapping) ordered labeled subgraphs of a labeled graph as a function of its randomness deficiency (how far it falls short of the maximum possible Kolmogorov complexity) and (ii) a new elementary proof for the number of unlabeled graphs. △ Less

Submitted 18 October, 2001; originally announced October 2001.

Comments: LaTeX 9 pages

MSC Class: 05C78; 94A17; 05C80; 05C70; 05C30; 05C35; 68R99

Journal ref: H. Buhrman, M. Li, J. Tromp and P.M.B. Vitanyi, Kolmogorov random graphs and the incompressibility method, SIAM J. Comput., 29:2(2000), 590--599

Randomized Two-Process Wait-Free Test-and-Set

Authors: John Tromp, Paul Vitanyi

Abstract: We present the first explicit, and currently simplest, randomized algorithm for 2-process wait-free test-and-set. It is implemented with two 4-valued single writer single reader atomic variables. A test-and-set takes at most 11 expected elementary steps, while a reset takes exactly 1 elementary step. Based on a finite-state analysis, the proofs of correctness and expected length are compressed i… ▽ More We present the first explicit, and currently simplest, randomized algorithm for 2-process wait-free test-and-set. It is implemented with two 4-valued single writer single reader atomic variables. A test-and-set takes at most 11 expected elementary steps, while a reset takes exactly 1 elementary step. Based on a finite-state analysis, the proofs of correctness and expected length are compressed into one table. △ Less

Submitted 15 March, 2002; v1 submitted 28 June, 2001; originally announced June 2001.

Comments: 9 pages, 4 figures, LaTeX source; Submitted

ACM Class: C.2; F.2.2

Time and Space Bounds for Reversible Simulation

Authors: Harry Buhrman, J. Tromp, Paul Vitanyi

Abstract: We prove a general upper bound on the tradeoff between time and space that suffices for the reversible simulation of irreversible computation. Previously, only simulations using exponential time or quadratic space were known. The tradeoff shows for the first time that we can simultaneously achieve subexponential time and subquadratic space. The boundary values are the exponential time with h… ▽ More We prove a general upper bound on the tradeoff between time and space that suffices for the reversible simulation of irreversible computation. Previously, only simulations using exponential time or quadratic space were known. The tradeoff shows for the first time that we can simultaneously achieve subexponential time and subquadratic space. The boundary values are the exponential time with hardly any extra space required by the Lange-McKenzie-Tapp method and the ($\log 3$)th power time with square space required by the Bennett method. We also give the first general lower bound on the extra storage space required by general reversible simulation. This lower bound is optimal in that it is achieved by some reversible simulations. △ Less

Submitted 19 April, 2001; v1 submitted 29 January, 2001; originally announced January 2001.

Comments: 11 pages LaTeX, Proc ICALP 2001, Lecture Notes in Computer Science, Vol xxx Springer-Verlag, Berlin, 2001

Journal ref: Journal of Physics A: Mathematical and General, 34(2001), 6821--6830.

CaB_6: a new semiconducting material for spin electronics

Authors: H. J. Tromp, P. van Gelderen, P. J. Kelly, G. Brocks, P. A. Bobbert

Abstract: Ferromagnetism was recently observed at unexpectedly high temperatures in La-doped CaB_6. The starting point of all theoretical proposals to explain this observation is a semimetallic electronic structure calculated for CaB_6 within the local density approximation. Here we report the results of parameter-free quasiparticle calculations of the single-particle excitation spectrum which show that C… ▽ More Ferromagnetism was recently observed at unexpectedly high temperatures in La-doped CaB_6. The starting point of all theoretical proposals to explain this observation is a semimetallic electronic structure calculated for CaB_6 within the local density approximation. Here we report the results of parameter-free quasiparticle calculations of the single-particle excitation spectrum which show that CaB_6 is not a semimetal but a semiconductor with a band gap of 0.8 eV. Magnetism in La_xCa_{1-x}B_6 occurs just on the metallic side of a Mott transition in the La-induced impurity band. △ Less

Submitted 7 November, 2000; originally announced November 2000.

Comments: 4 pages, 1 postscript figure

Journal ref: Phys. Rev. Lett. 87, 016401 (2001)

Algorithmic Statistics

Authors: Peter Gacs, John Tromp, Paul Vitanyi

Abstract: While Kolmogorov complexity is the accepted absolute measure of information content of an individual finite object, a similarly absolute notion is needed for the relation between an individual data sample and an individual model summarizing the information in the data, for example, a finite set (or probability distribution) where the data sample typically came from. The statistical theory based… ▽ More While Kolmogorov complexity is the accepted absolute measure of information content of an individual finite object, a similarly absolute notion is needed for the relation between an individual data sample and an individual model summarizing the information in the data, for example, a finite set (or probability distribution) where the data sample typically came from. The statistical theory based on such relations between individual objects can be called algorithmic statistics, in contrast to classical statistical theory that deals with relations between probabilistic ensembles. We develop the algorithmic theory of statistic, sufficient statistic, and minimal sufficient statistic. This theory is based on two-part codes consisting of the code for the statistic (the model summarizing the regularity, the meaningful information, in the data) and the model-to-data code. In contrast to the situation in probabilistic statistical theory, the algorithmic relation of (minimal) sufficiency is an absolute relation between the individual model and the individual data sample. We distinguish implicit and explicit descriptions of the models. We give characterizations of algorithmic (Kolmogorov) minimal sufficient statistic for all data samples for both description modes--in the explicit mode under some constraints. We also strengthen and elaborate earlier results on the ``Kolmogorov structure function'' and ``absolutely non-stochastic objects''--those rare objects for which the simplest models that summarize their relevant information (minimal sufficient statistics) are at least as complex as the objects themselves. We demonstrate a close relation between the probabilistic notions and the algorithmic ones. △ Less

Submitted 9 October, 2001; v1 submitted 30 June, 2000; originally announced June 2000.

Comments: LaTeX, 22 pages, 1 figure, with correction to the published journal version

MSC Class: 62B05; 62B10; 68Q32; 68Q30; 60AXX; 68T04

Journal ref: IEEE Transactions on Information Theory, Vol. 47, No. 6, September 2001, pp 2443-2463

Mutual Search

Authors: Harry Buhrman, Matthew Franklin, Juan A. Garay, Jaap-Henk Hoepman, John Tromp, Paul Vitanyi

Abstract: We introduce a search problem called ``mutual search'' where $k$ \agents, arbitrarily distributed over $n$ sites, are required to locate one another by posing queries of the form ``Anybody at site $i$?''. We ask for the least number of queries that is necessary and sufficient. For the case of two \agents using deterministic protocols we obtain the following worst-case results: In an oblivious se… ▽ More We introduce a search problem called ``mutual search'' where $k$ \agents, arbitrarily distributed over $n$ sites, are required to locate one another by posing queries of the form ``Anybody at site $i$?''. We ask for the least number of queries that is necessary and sufficient. For the case of two \agents using deterministic protocols we obtain the following worst-case results: In an oblivious setting (where all pre-planned queries are executed) there is no savings: $n-1$ queries are required and are sufficient. In a nonoblivious setting we can exploit the paradigm of ``no news is also news'' to obtain significant savings: in the synchronous case $0.586n$ queries suffice and $0.536n$ queries are required; in the asynchronous case $0.896n$ queries suffice and a fortiori 0.536 queries are required; for $o(\sqrt{n})$ \agents using a deterministic protocol less than $n$ queries suffice; there is a simple randomized protocol for two \agents with worst-case expected $0.5n$ queries and all randomized protocols require at least $0.125n$ worst-case expected queries. The graph-theoretic framework we formulate for expressing and analyzing algorithms for this problem may be of independent interest. △ Less

Submitted 2 February, 1999; originally announced February 1999.

Comments: 18 pages, Latex, 5 figures, J. Assoc. Comp. Mach., To appear

ACM Class: F.2,C.2,E,1,D.4.4

Reversible Simulation of Irreversible Computation by Pebble Games

Authors: Ming Li, John Tromp, Paul Vitanyi

Abstract: Reversible simulation of irreversible algorithms is analyzed in the stylized form of a `reversible' pebble game. While such simulations incur little overhead in additional computation time, they use a large amount of additional memory space during the computation. The reacheable reversible simulation instantaneous descriptions (pebble configurations) are characterized completely. As a corollary… ▽ More Reversible simulation of irreversible algorithms is analyzed in the stylized form of a `reversible' pebble game. While such simulations incur little overhead in additional computation time, they use a large amount of additional memory space during the computation. The reacheable reversible simulation instantaneous descriptions (pebble configurations) are characterized completely. As a corollary we obtain the reversible simulation by Bennett and that among all simulations that can be modelled by the pebble game, Bennett's simulation is optimal in that it uses the least auxiliary space for the greatest number of simulated steps. One can reduce the auxiliary storage overhead incurred by the reversible simulation at the cost of allowing limited erasing leading to an irreversibility-space tradeoff. We show that in this resource-bounded setting the limited erasing needs to be performed at precise instants during the simulation. We show that the reversible simulation can be modified so that it is applicable also when the simulated computation time is unknown. △ Less

Submitted 11 April, 1997; v1 submitted 6 March, 1997; originally announced March 1997.

Comments: 11 pages, Latex, Submitted to Physica D

Report number: CWI Tech Report 1996

Journal ref: Physica D120 (1998) 168-176