Bayesian analysis of multimessenger M-R data with interpolated hybrid EoS

We introduce a family of equations of state (EoS) for hybrid neutron star (NS) matter that is obtained by a two-zone parabolic interpolation between a soft hadronic EoS at low densities and a set of stiff quark matter EoS at high densities within a finite region of chemical potentials $$\mu _H< \mu < \mu _Q$$ μ H < μ < μ Q . Fixing the hadronic EoS as the APR one and choosing the color-superconducting, nonlocal NJL model with two free parameters for the quark phase, we perform Bayesian analyses with this two-parameter family of hybrid EoS. Using three different sets of observational constraints that include the mass of PSR J0740+6620, the tidal deformability for GW170817, and the mass-radius relation for PSR J0030+0451 from NICER as obligatory (set 1), while set 2 uses the possible upper limit on the maximum mass from GW170817 as an additional constraint and set 3 instead of the possibility that the lighter object in the asymmetric binary merger GW190814 is a neutron star. We confirm that in any case, the quark matter phase has to be color superconducting with the dimensionless diquark coupling approximately fulfilling the Fierz relation $$\eta _D=0.75$$ η D = 0.75 and the most probable solutions exhibiting a proportionality between $$\eta _D$$ η D and $$\eta _V$$ η V , the coupling of the repulsive vector interaction that is required for a sufficiently large maximum mass. We used the Bayesian analysis to investigate with the method of fictitious measurements the consequences of anticipating different radii for the massive $$2~M_\odot $$ 2 M ⊙ PSR J0740+6220 for the most likely equation of state. With the actual outcome of the NICER radius measurement on PSR J0740+6220 we could conclude that for the most likely hybrid star EoS would not support a maximum mass as large as $$2.5~M_\odot $$ 2.5 M ⊙ so that the event GW190814 was a binary black hole merger.

A Two-Way Nesting Unstructured Quadrilateral Grid, Finite-Differencing, Estuarine and Coastal Ocean Model with High-Order Interpolation Schemes

The balance between the need of improving horizontal resolution in simulating local small-scale ocean processes and computational costs makes it desirable to refine model mesh locally. A three-dimensional, two-way nesting unstructured quadrilateral grid, primitive equations, finite-differencing, estuarine and coastal ocean model is developed for multi-scale modeling. Because the model grid is capable of multi-area nesting and multi-level refinement at each subdomain, the model is highly compatible with simulations involved in complex topography and studies of local small-scale ocean processes. The two-way information exchange is achieved by a virtual grid method, and its basic idea is to implement numerical integrations of variables at nesting interfaces with the support of virtual grid variables, which are interpolated or updated from actual grid variables. The model is novel for two interpolation schemes: the high-order spatial interpolation at the middle temporal level (HSIMT) parabolic interpolation scheme and HSIMT advection-equivalent interpolation scheme, and they have high-order accuracy and good consistency with the advection scheme applied to solving the tracer equations. The conservation of both volume and tracer contents is ensured via a flux correction algorithm. The two original interpolation schemes are examined in an ideal salinity advection experiment in the peak preservation skill, stability, and conservation properties. A realistic application to the Deep Waterway Project area in the Changjiang Estuary showed that the nested grid model can reproduce the hydrodynamic processes at the observed sites successfully while it failed to maintain the performance with the structured grid model in simulating the variance of salinity, for which the enforced conservation had primary responsibility. The HSIMT parabolic interpolation scheme was distinguished from other schemes for its outstanding performances in simulating salinity.

An Efficient Method for Forming Parabolic Curves and Surfaces

A new method for the formation of parabolic curves and surfaces is proposed. It does not impose restrictions on the relative positions in space of the sequence of reference points relative to each other, meaning it compares favorably with other prototypes. The disadvantages of the Overhauser and Brever–Anderson methods are noted. The method allows one to effectively form and edit curves and surfaces when changing the coordinates of any given point. This positive effect is achieved due to the appropriate choice of local coordinate systems for the mathematical description of each parabola, which together define a composite interpolation curve or surface. The paper provides a detailed mathematical description of the method of parabolic interpolation of curves and surfaces based on the use of matrix calculations. Analytical descriptions of a composite parabolic curve and its first and second derivatives are given, and continuity analysis of these factors is carried out. For the matrix of points of the defining polyhedron, expressions are presented that describe the corresponding surfaces, as well as the unit normal at any point. The comparative table of the required number of pseudo-codes for calculating the coordinates of one point for constructing a parabolic curve for the three methods is given.