skip to main content

Title: Efficient spherical surface integration of Gauss functions in three-dimensional spherical coordinates and the solution for the modified Bessel function of the first kind
Award ID(s):
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of Mathematical Chemistry
Page Range / eLocation ID:
505 to 515
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We describe the discovery of an archaeal virus, one that infects archaea, tentatively named Thermoproteus spherical piliferous virus 1 (TSPV1), which was purified from a Thermoproteales host isolated from a hot spring in Yellowstone National Park (USA). TSPV1 packages an 18.65-kb linear double-stranded DNA (dsDNA) genome with 31 open reading frames (ORFs), whose predicted gene products show little homology to proteins with known functions. A comparison of virus particle morphologies and gene content demonstrates that TSPV1 is a new member of the Globuloviridae family of archaeal viruses. However, unlike other Globuloviridae members, TSPV1 has numerous highly unusual filaments decorating its surface, which can extend hundreds of micrometers from the virion. To our knowledge, similar filaments have not been observed in any other archaeal virus. The filaments are remarkably stable, remaining intact across a broad range of temperature and pH values, and they are resistant to chemical denaturation and proteolysis. A major component of the filaments is a glycosylated 35-kDa TSPV1 protein (TSPV1 GP24). The filament protein lacks detectable homology to structurally or functionally characterized proteins. We propose, given the low host cell densities of hot spring environments, that the TSPV1 filaments serve to increase the probability of virus attachment and entry into host cells. IMPORTANCE High-temperature environments have proven to be an important source for the discovery of new archaeal viruses with unusual particle morphologies and gene content. Our isolation of Thermoproteus spherical piliferous virus 1 (TSPV1), with numerous filaments extending from the virion surface, expands our understanding of viral diversity and provides new insight into viral replication in high-temperature environments. 
    more » « less
  2. Abstract

    We present an exactγ= 5/3 spherical accretion solution that modifies the Bondi boundary condition ofρconstantasr→ ∞ toρ→ 0 asr→ ∞. This change allows for simple power-law solutions on the density and infall velocity fields, ranging from a cold empty freefall condition where pressure tends to zero, to a hot hydrostatic equilibrium limit with no infall velocity. As in the case of the Bondi solution, a maximum accretion rate appears. As in theγ= 5/3 case of the Bondi solution, no sonic radius appears, this time however, because the flow is always characterized by a constant Mach number. This number equals 1 for the case of the maximum accretion rate, diverges toward the cold empty state, and becomes subsonic toward the hydrostatic equilibrium limit. It can be shown that in the limitr→ 0, the Bondi solution tends to the new solution presented, extending the validity of the Bondi accretion value to cases where the accretion density profile does not remain at a fixed constant value out to infinity. We then explore small deviations from sphericity and the presence of angular momentum through an analytic perturbative analysis. Such perturbed solutions yield a rich phenomenology through density and velocity fields in terms of Legendre polynomials, which we begin to explore for simple angular velocity boundary conditions having zeros on the plane and pole. The new solution presented provides complementary physical insight into accretion problems in general.

    more » « less
  3. Abstract The SPECFEM3D_Cartesian code package is widely used in simulating seismic wave propagation on local and regional scales due to its computational efficiency compared with the one-chunk version of the SPECFEM3D_Globe code. In SPECFEM3D_Cartesian, the built-in meshing tool maps a spherically curved cube to a rectangular cube using the Universal Transverse Mercator projection (UTM). Meanwhile, the geodetic east, north, and up directions are assigned as the local x–y–z directions. This causes coordinate orientation issues in simulating waveform propagation in regions larger than 6° × 6° or near the Earth’s polar regions. In this study, we introduce a new code package, named Cartesian Meshing Spherical Earth (CMSE), that can accurately mesh the 3D geometry of the Earth’s surface under the Cartesian coordinate frame, while retaining the geodetic directions. To benchmark our new package, we calculate the residual amplitude of the CMSE synthetics with respect to the reference synthetics calculated by SPECFEM3D_Globe. In the regional scale simulations with an area of 1300 km × 1300 km, we find a maximum of 5% amplitude residual for the SPECFEM3D_Cartesian synthetics using the mesh generated by the CMSE, much smaller than the maximum amplitude residual of 100% for the synthetics based on its built-in meshing tool. Therefore, our new meshing tool CMSE overcomes the limitations of the internal mesher used by SPECFEM3D_Cartesian and can be used for more accurate waveform simulations in larger regions beyond one UTM zone. Furthermore, CMSE can deal with regions at the south and north poles that cannot be handled by the UTM projection. Although other external code packages can be used to mesh the curvature of the Earth, the advantage of the CMSE code is that it is open-source, easy to use, and fully integrated with SPECFEM3D_Cartesian. 
    more » « less