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

Lee Lindblom

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-phphysics.comp-phmath-phmath.DGmath.MPmath.NANumerical Analysis

Catalog footprint

What is connected

21works
9topics
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

21 published item(s)

preprint2024arXiv

Building Three-Dimensional Differentiable Manifolds Numerically II: Limitations

Methods were developed in Ref. [1] for constructing reference metrics (and from them differentiable structures) on three-dimensional manifolds with topologies specified by suitable triangulations. This note generalizes those methods by expanding the class of suitable triangulations, significantly increasing the number of manifolds to which these methods apply. These new results show that this expanded class of triangulations is still a small subset of all possible triangulations. This demonstrates that fundamental changes to these methods are needed to further expand the collection of manifolds on which differentiable structures can be constructed numerically.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Building Three-Dimensional Differentiable Manifolds Numerically

A method is developed here for building differentiable three-dimensional manifolds on multicube structures. This method constructs a sequence of reference metrics that determine differentiable structures on the cubic regions that serve as non-overlapping coordinate charts on these manifolds. It uses solutions to the two- and three-dimensional biharmonic equations in a sequence of steps that increase the differentiability of the reference metrics across the interfaces between cubic regions. This method is algorithmic and has been implemented in a computer code that automatically generates these reference metrics. Examples of three-manifolds constructed in this way are presented here, including representatives from five of the eight Thurston geometrization classes, plus the well-known Hantzsche-Wendt, the Poincare dodecahedral space, and the Seifert-Weber space.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

First Constraining Upper Limits on Gravitational Wave Emission from NS 1987A in SNR 1987A

We report on a search for continuous gravitational waves (GWs) from NS 1987A, the neutron star born in SN 1987A. The search covered a frequency band of 75-275 Hz, included a wide range of spin-down parameters for the first time, and coherently integrated 12.8 days of data below 125 Hz and 8.7 days of data above 125 Hz from the second Advanced LIGO observing run. We found no astrophysical signal. We set upper limits on GW emission as tight as an intrinsic strain of $2\times10^{-25}$ at 90\% confidence. The large spin-down parameter space makes this search the first astrophysically consistent one for continuous GWs from NS 1987A. Our upper limits are the first consistent ones to beat an analog of the spin-down limit based on the age of the neutron star, and hence are the first GW observations to put new constraints on NS 1987A.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Directed searches for continuous gravitational waves from twelve supernova remnants in data from Advanced LIGO's second observing run

We describe directed searches for continuous gravitational waves from twelve well localized non-pulsing candidate neutron stars in young supernova remnants using data from Advanced LIGO's second observing run. We assumed that each neutron star is isolated and searched a band of frequencies from 15 to 150\,Hz, consistent with frequencies expected from known young pulsars. After coherently integrating spans of data ranging from 12.0 to 55.9 days using the F-statistic and applying data-based vetoes, we found no evidence of astrophysical signals. We set upper limits on intrinsic gravitational wave amplitude in some cases stronger than 10^{-25}, generally about a factor of two better than upper limits on the same objects from Advanced LIGO's first observing run.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constructing Reference Metrics on Multicube Representations of Arbitrary Manifolds

Reference metrics are used to define the differential structure on multicube representations of manifolds, i.e., they provide a simple and practical way to define what it means globally for tensor fields and their derivatives to be continuous. This paper introduces a general procedure for constructing reference metrics automatically on multicube representations of manifolds with arbitrary topologies. The method is tested here by constructing reference metrics for compact, orientable two-dimensional manifolds with genera between zero and five. These metrics are shown to satisfy the Gauss-Bonnet identity numerically to the level of truncation error (which converges toward zero as the numerical resolution is increased). These reference metrics can be made smoother and more uniform by evolving them with Ricci flow. This smoothing procedure is tested on the two-dimensional reference metrics constructed here. These smoothing evolutions (using volume-normalized Ricci flow with DeTurck gauge fixing) are all shown to produce reference metrics with constant scalar curvatures (at the level of numerical truncation error).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Model Waveform Accuracy Requirements for the Allen $χ^2$ Discriminator

This paper derives accuracy standards for model gravitational waveforms required to ensure proper use of the Allen $χ^2$ discriminator in gravitational wave (GW) data analysis. These standards are different from previously established requirements for detection and waveform parameter measurement based on signal-to-noise optimization. We present convenient formulae for evaluating and interpreting the contribution of model errors to measured values of this $χ^2$ statistic. The new accuracy standards derived here are needed to ensure the reliability of measured values of the Allen $χ^2$ statistic, both in their traditional role as vetoes and in their current role as elements in evaluating the significance of candidate detections.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spectral Approach to the Relativistic Inverse Stellar Structure Problem II

The inverse stellar structure problem determines the equation of state of the matter in stars from a knowledge of their macroscopic observables (e.g. their masses and radii). This problem was solved in a previous paper by constructing a spectral representation of the equation of state whose stellar models match a prescribed set of macroscopic observables. This paper improves and extends that work in two significant ways: i) The method is made more robust by accounting for an unexpected feature of the enthalpy based representations of the equations of state used in this work. After making the appropriate modifications, accurate initial guesses for the spectral parameters are no longer needed so Monte-Carlo techniques can now be used to ensure the best fit to the observables. ii) The method is extended here to use masses and tidal deformabilities (which will be measured by gravitational wave observations of neutron-star mergers) as the macroscopic observables instead of masses and radii. The accuracy and reliability of this extended and more robust spectral method is evaluated in this paper using mock data for observables from stars based on 34 different theoretical models of the high density neutron-star equation of state. In qualitative agreement with earlier work, these tests suggest the high density part of the neutron-star equation of state could be determined at the few-percent accuracy level using high quality measurements of the masses and radii (or masses and tidal deformabilities) of just two or three neutron stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Solving Einstein's Equation Numerically on Manifolds With Arbitrary Spatial Topologies

This paper develops a method for solving Einstein's equation numerically on multi-cube representations of manifolds with arbitrary spatial topologies. This method is designed to provide a set of flexible, easy to use computational procedures that make it possible to explore the never before studied properties of solutions to Einstein's equation on manifolds with arbitrary toplogical structures. A new covariant, first-order symmetric-hyperbolic representation of Einstein's equation is developed for this purpose, along with the needed boundary conditions at the interfaces between adjoining cubic regions. Numerical tests are presented that demonstrate the long-term numerical stability of this method for evolutions of a complicated, time-dependent solution of Einstein's equation coupled to a complex scalar field on a manifold with spatial topology S^3. The accuracy of these numerical test solutions is evaluated by performing convergence studies and by comparing the full non-linear numerical results to the analytical perturbation solutions, which are also derived here.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Relativistic Inverse Stellar Structure Problem

The observable macroscopic properties of relativistic stars (whose equations of state are known) can be predicted by solving the stellar structure equations that follow from Einstein's equation. For neutron stars, however, our knowledge of the equation of state is poor, so the direct stellar structure problem can not be solved without modeling the highest density part of the equation of state in some way. This talk will describe recent work on developing a model independent approach to determining the high-density neutron-star equation of state by solving an inverse stellar structure problem. This method uses the fact that Einstein's equation provides a deterministic relationship between the equation of state and the macroscopic observables of the stars which are composed of that material. This talk illustrates how this method will be able to determine the high-density part of the neutron-star equation of state with few percent accuracy when high quality measurements of the masses and radii of just two or three neutron stars become available. This talk will also show that this method can be used with measurements of other macroscopic observables, like the masses and tidal deformabilities, which can (in principle) be measured by gravitational wave observations of binary neutron-star mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Spectral Approach to the Relativistic Inverse Stellar Structure Problem

A new method for solving the relativistic inverse stellar structure problem is presented. This method determines a spectral representation of the unknown high density portion of the stellar equation of state from a knowledge of the total masses M and radii R of the stars. Spectral representations of the equation of state are very efficient, generally requiring only a few spectral parameters to achieve good accuracy. This new method is able, therefore, to determine the high density equation of state quite accurately from only a few accurately measured [M,R] data points. This method is tested here by determining the equations of state from mock [M,R] data computed from tabulated "realistic" neutron-star equations of state. The spectral equations of state obtained from these mock data are shown to agree on average with the originals to within a few percent (over the entire high density range of the neutron-star interior) using only two [M,R] data points. Higher accuracies are achieved when more data are used. The accuracies of the equations of state determined in these examples are shown to be nearly optimal, in the sense that their errors are comparable to the errors of the best-fit spectral representations of these realistic equations of state.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solving Partial Differential Equations Numerically on Manifolds with Arbitrary Spatial Topologies

A multi-cube method is developed for solving systems of elliptic and hyperbolic partial differential equations numerically on manifolds with arbitrary spatial topologies. It is shown that any three-dimensional manifold can be represented as a set of non-overlapping cubic regions, plus a set of maps to identify the faces of adjoining regions. The differential structure on these manifolds is fixed by specifying a smooth reference metric tensor. Matching conditions that ensure the appropriate levels of continuity and differentiability across region boundaries are developed for arbitrary tensor fields. Standard numerical methods are then used to solve the equations with the appropriate boundary conditions, which are determined from these inter-region matching conditions. Numerical examples are presented which use pseudo-spectral methods to solve simple elliptic equations on multi-cube representations of manifolds with the topologies T^3, S^2 x S^1 and S^3. Examples are also presented of numerical solutions of simple hyperbolic equations on multi-cube manifolds with the topologies R x T^3, R x S^2 x S^1 and R x S^3.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Improved Time-Domain Accuracy Standards for Model Gravitational Waveforms

Model gravitational waveforms must be accurate enough to be useful for detection of signals and measurement of their parameters, so appropriate accuracy standards are needed. Yet these standards should not be unnecessarily restrictive, making them impractical for the numerical and analytical modelers to meet. The work of Lindblom, Owen, and Brown [Phys. Rev. D 78, 124020 (2008)] is extended by deriving new waveform accuracy standards which are significantly less restrictive while still ensuring the quality needed for gravitational-wave data analysis. These new standards are formulated as bounds on certain norms of the time-domain waveform errors, which makes it possible to enforce them in situations where frequency-domain errors may be difficult or impossible to estimate reliably. These standards are less restrictive by about a factor of 20 than the previously published time-domain standards for detection, and up to a factor of 60 for measurement. These new standards should therefore be much easier to use effectively.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Spectral Representations of Neutron-Star Equations of State

Methods are developed for constructing spectral representations of cold (barotropic) neutron-star equations of state. These representations are faithful in the sense that every physical equation of state has a representation of this type, and conversely every such representation satisfies the minimal thermodynamic stability criteria required of any physical equation of state. These spectral representations are also efficient, in the sense that only a few spectral coefficients are generally required to represent neutron-star equations of state quiet accurately. This accuracy and efficiency is illustrated by constructing spectral fits to a large collection of "realistic" neutron-star equations of state.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Optimal Calibration Accuracy for Gravitational Wave Detectors

Calibration errors in the response function of a gravitational wave detector degrade its ability to detect and then to measure the properties of any detected signals. This paper derives the needed levels of calibration accuracy for each of these data-analysis tasks. The levels derived here are optimal in the sense that lower accuracy would result in missed detections and/or a loss of measurement precision, while higher accuracy would be made irrelevant by the intrinsic noise level of the detector. Calibration errors affect the data-analysis process in much the same way as errors in theoretical waveform templates. The optimal level of calibration accuracy is expressed therefore as a joint limit on modeling and calibration errors: increased accuracy in one reduces the accuracy requirement in the other.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Simulations of Binary Black Hole Mergers Using Spectral Methods

Several improvements in numerical methods and gauge choice are presented that make it possible now to perform simulations of the merger and ringdown phases of "generic" binary black-hole evolutions using the pseudo-spectral evolution code SpEC. These improvements include the use of a new damped-wave gauge condition, a new grid structure with appropriate filtering that improves stability, and better adaptivity in conforming the grid structures to the shapes and sizes of the black holes. Simulations illustrating the success of these new methods are presented for a variety of binary black-hole systems. These include fairly ``generic'' systems with unequal masses (up to 2:1 mass ratios), and spins (with magnitudes up to 0.4 M^2) pointing in various directions.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A New Generalized Harmonic Evolution System

A new representation of the Einstein evolution equations is presented that is first order, linearly degenerate, and symmetric hyperbolic. This new system uses the generalized harmonic method to specify the coordinates, and exponentially suppresses all small short-wavelength constraint violations. Physical and constraint-preserving boundary conditions are derived for this system, and numerical tests that demonstrate the effectiveness of the constraint suppression properties and the constraint-preserving boundary conditions are presented.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Second-order rotational effects on the r-modes of neutron stars

Techniques are developed here for evaluating the r-modes of rotating neutron stars through second order in the angular velocity of the star. Second-order corrections to the frequencies and eigenfunctions for these modes are evaluated for neutron star models. The second-order eigenfunctions for these modes are determined here by solving an unusual inhomogeneous hyperbolic boundary-value problem. The numerical techniques developed to solve this unusual problem are somewhat non-standard and may well be of interest beyond the particular application here. The bulk-viscosity coupling to the r-modes, which appears first at second order, is evaluated. The bulk-viscosity timescales are found here to be longer than previous estimates for normal neutron stars, but shorter than previous estimates for strange stars. These new timescales do not substantially affect the current picture of the gravitational radiation driven instability of the r-modes either for neutron stars or for strange stars.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Stability of the r-modes in white dwarf stars

Stability of the r-modes in rapidly rotating white dwarf stars is investigated. Improved estimates of the growth times of the gravitational-radiation driven instability in the r-modes of the observed DQ Her objects are found to be longer (probably considerably longer) than 6x10^9y. This rules out the possibility that the r-modes in these objects are emitting gravitational radiation at levels that could be detectable by LISA. More generally it is shown that the r-mode instability can only be excited in a very small subset of very hot (T>10^6K), rather massive (M>0.9M_sun) and very rapidly rotating (P_min<P<1.2P_min) white dwarf stars. Further, the growth times of this instability are so long that these conditions must persist for a very long time (t>10^9y) to allow the amplitude to grow to a dynamically significant level. This makes it extremely unlikely that the r-mode instability plays a significant role in any real white dwarf stars.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Gravitational Radiation Instability in Hot Young Neutron Stars

We show that gravitational radiation drives an instability in hot young rapidly rotating neutron stars. This instability occurs primarily in the l=2 r-mode and will carry away most of the angular momentum of a rapidly rotating star by gravitational radiation. On the timescale needed to cool a young neutron star to about T=10^9 K (about one year) this instability can reduce the rotation rate of a rapidly rotating star to about 0.076Ω_K, where Ω_K is the Keplerian angular velocity where mass shedding occurs. In older colder neutron stars this instability is suppressed by viscous effects, allowing older stars to be spun up by accretion to larger angular velocities.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Gravitational Helioseismology?

The magnitudes of the external gravitational perturbations associated with the normal modes of the Sun are evaluated to determine whether these solar oscillations could be observed with the proposed Laser Interferometer Space Antenna (LISA), a network of satellites designed to detect gravitational radiation. The modes of relevance to LISA---the $l=2$, low-order $p$, $f$ and $g$-modes---have not been conclusively observed to date. We find that the energy in these modes must be greater than about $10^{30} \rm{ergs}$ in order to be observable above the LISA detector noise. These mode energies are larger than generally expected, but are much smaller than the current observational upper limits. LISA may be confusion-limited at the relevant frequencies due to the galactic background from short-period white dwarf binaries. Present estimates of the number of these binaries would require the solar modes to have energies above about $10^{33} \rm{ergs}$ to be observable by LISA.

0 citations0 reviewsNo code linkedOpen access