Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

An advanced expiratory circuit for the recovery of perfluorocarbon liquid from non-saturated perfluorocarbon vapour during partial liquid ventilation: an experimental model

BACKGROUND: The loss of perfluorocarbon (PFC) vapour in the expired gases during partial liquid ventilation should be minimized both to prevent perfluorocarbon vapour entering the atmosphere and to re-use the recovered PFC liquid. Using a substantially modified design of our previously described condenser, we aimed to determine how much perfluorocarbon liquid could be recovered from gases containing PFC and water vapour, at concentrations found during partial liquid ventilation, and to determine if the amount recovered differed with background flow rate (at flow rates suitable for use in neonates). METHODS: The expiratory line of a standard ventilator circuit set-up was mimicked, with the addition of two condensers. Perfluorocarbon (30 mL of FC-77) and water vapour, at concentrations found during partial liquid ventilation, were passed through the circuit at a number of flow rates and the percentage recovery of the liquids measured. RESULTS: From 14.2 mL (47%) to 27.3 mL (91%) of the infused 30 mL of FC-77 was recovered at the flow rates studied. Significantly higher FC-77 recovery was obtained at lower flow rates (ANOVA with Bonferroni's multiple comparison test, p < 0.0001). As a percentage of the theoretical maximum recovery, 64 to 95% of the FC-77 was recovered. Statistically significantly less FC-77 was recovered at 5 Lmin(-1 )(ANOVA with Bonferroni's multiple comparison test, p < 0.0001). Amounts of perfluorocarbon vapour recovered were 47%, 50%, 81% and 91% at flow rates of 10, 5, 2 and 1 Lmin(-1), respectively. CONCLUSION: Using two condensers in series 47% to 91% of perfluorocarbon liquid can be recovered, from gases containing perfluorocarbon and water vapour, at concentrations found during partial liquid ventilation

On the (not so) constant proportional trade-off in TTO

Abstract. Purpose: The linear and power QALY models require that people in Time Trade-off (TTO) exercises sacrifice the same proportion of lifetime to obtain a health improvement, irrespective of the absolute amount. However, evidence on these constant proportional trade-offs (CPTOs) is mixed, indicating that these versions of the QALY model do not represent preferences. Still, it may be the case that a more general version of the QALY model represents preferences. This version has the property that people want to sacrifice the same proportion of utilities of lifetime for a health improvement, irrespective of the amount of this lifetime. Methods: We use a new method to correct TTO scores for utility of life duration and test whether decision makers trade off utility of duration and quality at the same rate irrespective of duration. Results: We find a robust violation of CPTO for both uncorrected and corrected TTO scores. Remarkably, we find higher values for longer durations, contrary to most previous studies. This represents the only study correcting for utility of life duration to find such a violation. Conclusions: It seems that the trade-off of life years is indeed not so constantly proportional and, therefore, that health state valuations depend on durations

Structural shape optimization using Cartesian grids and automatic h-adaptive mesh projection

[EN] We present a novel approach to 3D structural shape optimization that leans on an Immersed Boundary Method. A boundary tracking strategy based on evaluating the intersections between a fixed Cartesian grid and the evolving geometry sorts elements as internal, external and intersected. The integration procedure used by the NURBS-Enhanced Finite Element Method accurately accounts for the nonconformity between the fixed embedding discretization and the evolving structural shape, avoiding the creation of a boundary-fitted mesh for each design iteration, yielding in very efficient mesh generation process. A Cartesian hierarchical data structure improves the efficiency of the analyzes, allowing for trivial data sharing between similar entities or for an optimal reordering of thematrices for the solution of the system of equations, among other benefits. Shape optimization requires the sufficiently accurate structural analysis of a large number of different designs, presenting the computational cost for each design as a critical issue. The information required to create 3D Cartesian h- adapted mesh for new geometries is projected from previously analyzed geometries using shape sensitivity results. Then, the refinement criterion permits one to directly build h-adapted mesh on the new designs with a specified and controlled error level. Several examples are presented to show how the techniques here proposed considerably improve the computational efficiency of the optimization process.The authors wish to thank the Spanish Ministerio de Economia y Competitividad for the financial support received through the project DPI2013-46317-R and the FPI program (BES-2011-044080), and the Generalitat Valenciana through the project PROMETEO/2016/007.Marco, O.; Ródenas, J.; Albelda Vitoria, J.; Nadal, E.; Tur Valiente, M. (2017). Structural shape optimization using Cartesian grids and automatic h-adaptive mesh projection. Structural and Multidisciplinary Optimization. 1-21. https://doi.org/10.1007/s00158-017-1875-1S121MATLAB version 8.3.0.532 (R2014a) (2014) Documentation. The Mathworks, Inc., Natick, MassachusettsAbel JF, Shephard MS (1979) An algorithm for multipoint constraints in finite element analysis. Int J Numer Methods Eng 14(3):464–467Amestoy P, Davis T, Duff I (1996) An approximate minimum degree ordering algorithm. SIAM J Matrix Anal Appl 17(4):886–905Barth W, Stürzlinger W (1993) Efficient ray tracing for Bezier and B-spline surfaces. Comput Graph 17 (4):423–430Bennett J A, Botkin M E (1985) Structural shape optimization with geometric problem description and adaptive mesh refinement. AIAA J 23(3):459–464Braibant V, Fleury C (1984) Shape optimal design using b-splines. Comput Methods Appl Mech Eng 44 (3):247–267Bugeda G, Oliver J (1993) A general methodology for structural shape optimization problems using automatic adaptive remeshing. Int J Numer Methods Eng 36(18):3161–3185Bugeda G, Ródenas J J, Oñate E (2008) An integration of a low cost adaptive remeshing strategy in the solution of structural shape optimization problems using evolutionary methods. Comput Struct 86(13–14):1563–1578Chang K, Choi K K (1992) A geometry-based parameterization method for shape design of elastic solids. Mech Struct Mach 20(2):215–252Cho S, Ha S H (2009) Isogeometric shape design optimization: exact geometry and enhanced sensitivity. Struct Multidiscip Optim 38(1):53–70Belegundu D, Zhang YMS, Salagame R (1991) The natural approach for shape optimization with mesh distortion control. Tech. rep., Penn State UniversityDavis T A, Gilbert J R, Larimore S, Ng E (2004) An approximate column minimum degree ordering algorithm. ACM Trans Math Softw 30(3):353–376Doctor L J, Torborg J G (1981) Display techniques for octree-encoded objects. IEEE Comput Graph Appl 1(3):29–38Dunning P D, Kim H A, Mullineux G (2011) Investigation and improvement of sensitivity computation using the area-fraction weighted fixed grid FEM and structural optimization. Finite Elem Anal Des 47(8):933–941Düster A, Parvizian J, Yang Z, Rank E (2008) The finite cell method for three-dimensional problems of solid mechanics. Comput Methods Appl Mech Eng 197(45-48):3768–3782Escobar J M, Montenegro R, Rodríguez E, Cascón J M (2014) The meccano method for isogeometric solid modeling and applications. Eng Comput 30(3):331–343Farhat C, Lacour C, Rixen D (1998) Incorporation of linear multipoint constraints in substructure based iterative solvers. Part 1: a numerically scalable algorithm. Int J Numer Methods Eng 43(6):997–1016Fries T P, Omerović S (2016) Higher-order accurate integration of implicit geometries. Int J Numer Methods Eng 106(5):323–371Fuenmayor F J, Oliver J L (1996) Criteria to achieve nearly optimal meshes in the h-adaptive finite element mehod. Int J Numer Methods Eng 39(23):4039–4061Fuenmayor F J, Oliver J L, Ródenas J J (1997) Extension of the Zienkiewicz-Zhu error estimator to shape sensitivity analysis. Int J Numer Methods Eng 40(8):1413–1433García-Ruíz M J, Steven G P (1999) Fixed grid finite elements in elasticity problems. Eng Comput 16 (2):145–164Gill P, Murray W, Saunders M, Wright M (1984) Procedures for optimization problems with a mixture of bounds and general linear constraints. ACM Trans Math Software 10:282–298González-Estrada O A, Nadal E, Ródenas J J, Kerfriden P, Bordas S P A, Fuenmayor F J (2014) Mesh adaptivity driven by goal-oriented locally equilibrated superconvergent patch recovery. Comput Mech 53(5):957–976Ha S H, Choi K K, Cho S (2010) Numerical method for shape optimization using T-spline based isogeometric method. Struct Multidiscip Optim 42(3):417–428Haftka R T, Grandhi R V (1986) Structural shape optimization: A survey. Comput Methods Appl Mech Eng 57(1):91–106Haslinger J, Jedelsky D (1996) Genetic algorithms and fictitious domain based approaches in shape optimization. Struc Optim 12:257–264Hughes T J R, Cottrell J A, Bazilevs Y (2005) Isogeometric Analysis: CAD, Finite Elements, NURBS, Exact Geometry, and Mesh Refinement. Comput Methods Appl Mech Eng 194:4135–4195Jackins C L, Tanimoto S L (1980) Oct-tree and their use in representing three-dimensional objects. Comput Graphics Image Process 14(3):249–270Kajiya J T (1982) Ray Tracing Parametric Patches. SIGGRAPH Comput Graph 16(3):245–254van Keulen F, Haftka R T, Kim N (2005) Review of options for structural design sensitivity analysis. Part I: linear systems. Comput Methods Appl Mech Eng 194(30-33):3213–3243Kibsgaard S (1992) Sensitivity analysis-the basis for optimization. Int J Numer Methods Eng 34(3):901–932Kikuchi N, Chung K Y, Torigaki T, Taylor J E (1986) Adaptive finite element methods for shape optimization of linearly elastic structures. Comput Methods Appl Mech Eng 57(1):67–89Kim N H, Chang Y (2005) Eulerian shape design sensitivity analysis and optimization with a fixed grid. Comput Methods Appl Mech Eng 194(30–33):3291–3314Kudela L, Zander N, Kollmannsberger S, Rank E (2016) Smart octrees: Accurately integrating discontinuous functions in 3d. Comput Methods Appl Mech Eng 306(1):406–426Kunisch K, Peichl G (1996) Numerical gradients for shape optimization based on embedding domain techniques. Comput Optim 18:95–114Li K, Qian X (2011) Isogeometric analysis and shape optimization via boundary integral. Computer-Aided Design 43(11):1427–1437Lian H, Kerfriden P, Bordas S P A (2016) Implementation of regularized isogeometric boundary element methods for gradient-based shape optimization in two-dimensional linear elasticity. Int J Numer Methods Eng 106 (12):972–1017Liu L, Zhang Y, Hughes T J R, Scott M A, Sederberg T W (2014) Volumetric T-spline Construction using Boolean Operations. Eng Comput 30(4):425–439Marco O, Sevilla R, Zhang Y, Ródenas J J, Tur M (2015) Exact 3D boundary representation in finite element analysis based on Cartesian grids independent of the geometry. Int J Numer Methods Eng 103:445–468Marco O, Ródenas J J, Fuenmayor FJ, Tur M (2017a) An extension of shape sensitivity analysis to an immersed boundary method based on cartesian grids. Computational Mechanics SubmittedMarco O, Ródenas J J, Navarro-Jiménez JM, Tur M (2017b) Robust h-adaptive meshing strategy for arbitrary cad geometries in a cartesian grid framework. Computers & Structures SubmittedMeagher D (1980) Octree Encoding: A New Technique for the Representation, Manipulation and Display of Arbitrary 3-D Objects by Computer. Tech. Rep. IPL-TR-80-11 I, Rensselaer Polytechnic InstituteMoita J S, Infante J, Mota C M, Mota C A (2000) Sensitivity analysis and optimal design of geometrically non-linear laminated plates and shells. Comput Struct 76(1–3):407–420Nadal E (2014) Cartesian Grid FEM (cgFEM): High Performance h-adaptive FE Analysis with Efficient Error Control. Application to Structural Shape Optimization. PhD Thesis. Universitat Politècnica de ValènciaNadal E, Ródenas J J, Albelda J, Tur M, Tarancón J E, Fuenmayor F J (2013) Efficient finite element methodology based on cartesian grids: application to structural shape optimization. Abstr Appl Anal 2013:1–19Najafi A R, Safdari M, Tortorelli D A, Geubelle P H (2015) A gradient-based shape optimization scheme using an interface-enriched generalized FEM. Comput Methods Appl Mech Eng 296:1–17Nguyen V P, Anitescu C, Bordas S P A, Rabczuk T (2015) Isogeometric analysis: An overview and computer implementation aspects. Math Comput Simul 117:89–116Nishita T, Sederberg TW, Kakimoto M (1990) Ray Tracing Trimmed Rational Surface Patches. SIGGRAPH Comput Graph 24(4):337–345Nocedal J, Wright SJ (2006) Numerical optimization, 2nd edn. Springer-Verlag, New YorkPandey P C, Bakshi P (1999) Analytical response sensitivity computation using hybrid finite elements. Comput Struct 71(5):525–534Parvizian J, Düster A, Rank E (2007) Finite Cell Method: h- and p- Extension for Embedded Domain Methods in Solid Mechanics. Comput Mech 41(1):121–133Peskin C S (1977) Numerical Analysis of Blood Flow in the Heart. J Comput Phys 25:220–252Poldneff M J, Rai I S, Arora J S (1993) Implementation of design sensitivity analysis for nonlinear structures. AIAA J 31(11):2137–2142Powell M (1983) Variable metric methods for constrained optimization. In: Bachem A, Grotschel M, Korte B (eds) Mathematical Programming: The State of the Art, Springer, Berlin, Heidelberg, pp 288–311Qian X (2010) Full analytical sensitivities in NURBS based isogeometric shape optimization. Comput Methods Appl Mech Eng 199(29–32):2059–2071Riehl S, Steinmann P (2014) An integrated approach to shape optimization and mesh adaptivity based on material residual forces. Comput Methods Appl Mech Eng 278:640–663Riehl S, Steinmann P (2016) On structural shape optimization using an embedding domain discretization technique. Int J Numer Methods Eng 109(9):1315–1343Ródenas J J, Tarancón J E, Albelda J, Roda A, Fuenmayor F J (2005) Hierarchical Properties in Elements Obtained by Subdivision: a Hierarquical h-adaptivity Program. In: Díez P, Wiberg N E (eds) Adaptive Modeling and Simulation, p 2005Ródenas J J, Corral C, Albelda J, Mas J, Adam C (2007a) Nested domain decomposition direct and iterative solvers based on a hierarchical h-adaptive finite element code. In: Runesson K, Díez P (eds) Adaptive Modeling and Simulation 2007, Internacional Center for Numerical Methods in Engineering (CIMNE), pp 206–209Ródenas J J, Tur M, Fuenmayor F J, Vercher A (2007b) Improvement of the superconvergent patch recovery technique by the use of constraint equations: the SPR-C technique. Int J Numer Methods Eng 70(6):705–727Ródenas J J, Bugeda G, Albelda J, Oñate E (2011) On the need for the use of error-controlled finite element analyses in structural shape optimization processes. Int J Numer Methods Eng 87(11):1105–1126Schillinger D, Ruess M (2015) The finite cell method: A review in the context of higher-order structural analysis of cad and image-based geometric models. Arch Comput Meth Eng 22(3):391– 455Sevilla R, Fernández-Méndez S, Huerta A (2011a) 3D-NURBS-enhanced Finite Element Method (NEFEM). Int J Numer Methods Eng 88(2):103–125Sevilla R, Fernández-Méndez S, Huerta A (2011b) Comparison of High-order Curved Finite Elements. Int J Numer Methods Eng 87(8):719–734Sevilla R, Fernández-Méndez S, Huerta A (2011c) NURBS-enhanced Finite Element Method (NEFEM): A Seamless Bridge Between CAD and FEM. Arch Comput Meth Eng 18(4):441–484Sweeney M, Bartels R (1986) Ray tracing free-form b-spline surfaces. IEEE Comput Graph Appl 6(2):41–49Toth D L (1985) On Ray Tracing Parametric Surfaces. SIGGRAPH Comput Graph 19(3):171–179Tur M, Albelda J, Nadal E, Ródenas J J (2014) Imposing dirichlet boundary conditions in hierarchical cartesian meshes by means of stabilized lagrange multipliers. Int J Numer Methods Eng 98(6):399–417Tur M, Albelda J, Marco O, Ródenas J J (2015) Stabilized Method to Impose Dirichlet Boundary Conditions using a Smooth Stress Field. Comput Methods Appl Mech Eng 296:352–375Yao T, Choi KK (1989) 3-d shape optimal design and automatic finite element regridding. Int J Numer Methods Eng 28(2):369–384Zhang L, Gerstenberger A, Wang X, Liu W K (2004) Immersed Finite Element Method. Comput Methods Appl Mech Eng 293(21):2051–2067Zhang Y, Wang W, Hughes T J R (2013) Conformal Solid T-spline Construction from Boundary T-spline Representations. Comput Mech 6(51):1051–1059Zienkiewicz O C, Zhu J Z (1987) A Simple Error Estimator and Adaptive Procedure for Practical Engineering Analysis. Int J Numer Methods Eng 24(2):337–35

BMP signaling components in embryonic transcriptomes of the hover fly Episyrphus balteatus (Syrphidae)

<p>Abstract</p> <p>Background</p> <p>In animals, signaling of Bone Morphogenetic Proteins (BMPs) is essential for dorsoventral (DV) patterning of the embryo, but how BMP signaling evolved with changes in embryonic DV differentiation is largely unclear. Based on the extensive knowledge of BMP signaling in <it>Drosophila melanogaster</it>, the morphological diversity of extraembryonic tissues in different fly species provides a comparative system to address this question. The closest relatives of <it>D. melanogaster </it>with clearly distinct DV differentiation are hover flies (Diptera: Syrphidae). The syrphid <it>Episyrphus balteatus </it>is a commercial bio-agent against aphids and has been established as a model organism for developmental studies and chemical ecology. The dorsal blastoderm of <it>E. balteatus </it>gives rise to two extraembryonic tissues (serosa and amnion), whereas in <it>D. melanogaster</it>, the dorsal blastoderm differentiates into a single extraembryonic epithelium (amnioserosa). Recent studies indicate that several BMP signaling components of <it>D. melanogaster</it>, including the BMP ligand Screw (Scw) and other extracellular regulators, evolved in the dipteran lineage through gene duplication and functional divergence. These findings raise the question of whether the complement of BMP signaling components changed with the origin of the amnioserosa.</p> <p>Results</p> <p>To search for BMP signaling components in <it>E. balteatus</it>, we generated and analyzed transcriptomes of freshly laid eggs (0-30 minutes) and late blastoderm to early germband extension stages (3-6 hours) using Roche/454 sequencing. We identified putative <it>E. balteatus </it>orthologues of 43% of all annotated <it>D. melanogaster </it>genes, including the genes of all BMP ligands and other BMP signaling components.</p> <p>Conclusion</p> <p>The diversification of several BMP signaling components in the dipteran linage of <it>D. melanogaster </it>preceded the origin of the amnioserosa.</p> <p>[Transcriptome sequence data from this study have been deposited at the NCBI Sequence Read Archive (SRP005289); individually assembled sequences have been deposited at GenBank (<ext-link ext-link-id="JN006969" ext-link-type="gen">JN006969</ext-link>-<ext-link ext-link-id="JN006986" ext-link-type="gen">JN006986</ext-link>).]</p

GW190814: Gravitational Waves from the Coalescence of a 23 Solar Mass Black Hole with a 2.6 Solar Mass Compact Object

We report the observation of a compact binary coalescence involving a 22.2–24.3 M⊙ black hole and a compact object with a mass of 2.50–2.67 M⊙ (all measurements quoted at the 90% credible level). The gravitational-wave signal, GW190814, was observed during LIGO's and Virgo's third observing run on 2019 August 14 at 21:10:39 UTC and has a signal-to-noise ratio of 25 in the three-detector network. The source was localized to 18.5 deg2 at a distance of Mpc; no electromagnetic counterpart has been confirmed to date. The source has the most unequal mass ratio yet measured with gravitational waves, , and its secondary component is either the lightest black hole or the heaviest neutron star ever discovered in a double compact-object system. The dimensionless spin of the primary black hole is tightly constrained to ≤0.07. Tests of general relativity reveal no measurable deviations from the theory, and its prediction of higher-multipole emission is confirmed at high confidence. We estimate a merger rate density of 1–23 Gpc−3 yr−1 for the new class of binary coalescence sources that GW190814 represents. Astrophysical models predict that binaries with mass ratios similar to this event can form through several channels, but are unlikely to have formed in globular clusters. However, the combination of mass ratio, component masses, and the inferred merger rate for this event challenges all current models of the formation and mass distribution of compact-object binaries

Erratum: "Searches for Gravitational Waves from Known Pulsars at Two Harmonics in 2015–2017 LIGO Data" (2019, ApJ, 879, 10)

This is a correction for 2019 ApJ 879 1

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO’s second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h95%0=3.47×10−25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering

Search for gravitational waves from Scorpius X-1 in the second Advanced LIGO observing run with an improved hidden Markov model

We present results from a semicoherent search for continuous gravitational waves from the low-mass x-ray binary Scorpius X-1, using a hidden Markov model (HMM) to track spin wandering. This search improves on previous HMM-based searches of LIGO data by using an improved frequency domain matched filter, the J-statistic, and by analyzing data from Advanced LIGO's second observing run. In the frequency range searched, from 60 to 650 Hz, we find no evidence of gravitational radiation. At 194.6 Hz, the most sensitive search frequency, we report an upper limit on gravitational wave strain (at 95% confidence) of h095%=3.47×10-25 when marginalizing over source inclination angle. This is the most sensitive search for Scorpius X-1, to date, that is specifically designed to be robust in the presence of spin wandering. © 2019 American Physical Society