-
On the logarithmic energy of solutions to the polynomial eigenvalue problem
Authors:
Diego Armentano,
Federico Carrasco,
Marcelo Fiori
Abstract:
In this paper, we compute the expected logarithmic energy of solutions to the polynomial eigenvalue problem for random matrices. We generalize some known results for the Shub-Smale polynomials, and the spherical ensemble. These two processes are the two extremal particular cases of the polynomial eigenvalue problem, and we prove that the logarithmic energy lies between these two cases. In particul…
▽ More
In this paper, we compute the expected logarithmic energy of solutions to the polynomial eigenvalue problem for random matrices. We generalize some known results for the Shub-Smale polynomials, and the spherical ensemble. These two processes are the two extremal particular cases of the polynomial eigenvalue problem, and we prove that the logarithmic energy lies between these two cases. In particular, the roots of the Shub-Smale polynomials are the ones with the lowest logarithmic energy of the family.
△ Less
Submitted 2 September, 2024; v1 submitted 20 August, 2024;
originally announced August 2024.
-
Delayed jet launching in binary neutron star mergers with realistic initial magnetic fields
Authors:
Ricard Aguilera-Miret,
Carlos Palenzuela,
Federico Carrasco,
Stephan Rosswog,
Daniele Viganò
Abstract:
We analyze a long-lived hyper-massive neutron star merger remnant (post-merger lifetime $>250$ ms) that has been obtained via large eddy simulations with a gradient subgrid-scale model. We find a clear helicoidal magnetic field structure that is governed by the toroidal component of the magnetic field. Although no jet emerges during the simulation time, we observe at late times a significant incre…
▽ More
We analyze a long-lived hyper-massive neutron star merger remnant (post-merger lifetime $>250$ ms) that has been obtained via large eddy simulations with a gradient subgrid-scale model. We find a clear helicoidal magnetic field structure that is governed by the toroidal component of the magnetic field. Although no jet emerges during the simulation time, we observe at late times a significant increase of the poloidal component of the magnetic field at all scales. We also compare with the results of several binary neutron star simulations with moderate resolution of $120$~m, that are evolved up to $50$~ms after the merger, which differ in terms of the initial topology and strength of the magnetic field. We find that the best choice is an isotropic small-scale magnetic field distribution that mimics the turbulent state that generically develops during the merger. This initial configuration reaches a closer agreement with our high-resolution simulation results than the purely dipolar large-scale fields that are commonly employed in these type of simulations. This provides a recipe to perform such simulations avoiding the computationally expensive grids required to faithfully capture the amplification of the magnetic field by Kelvin-Helmholtz instabilities.
△ Less
Submitted 29 July, 2024;
originally announced July 2024.
-
Generalizable Human Gaussians for Sparse View Synthesis
Authors:
Youngjoong Kwon,
Baole Fang,
Yixing Lu,
Haoye Dong,
Cheng Zhang,
Francisco Vicente Carrasco,
Albert Mosella-Montoro,
Jianjin Xu,
Shingo Takagi,
Daeil Kim,
Aayush Prakash,
Fernando De la Torre
Abstract:
Recent progress in neural rendering has brought forth pioneering methods, such as NeRF and Gaussian Splatting, which revolutionize view rendering across various domains like AR/VR, gaming, and content creation. While these methods excel at interpolating {\em within the training data}, the challenge of generalizing to new scenes and objects from very sparse views persists. Specifically, modeling 3D…
▽ More
Recent progress in neural rendering has brought forth pioneering methods, such as NeRF and Gaussian Splatting, which revolutionize view rendering across various domains like AR/VR, gaming, and content creation. While these methods excel at interpolating {\em within the training data}, the challenge of generalizing to new scenes and objects from very sparse views persists. Specifically, modeling 3D humans from sparse views presents formidable hurdles due to the inherent complexity of human geometry, resulting in inaccurate reconstructions of geometry and textures. To tackle this challenge, this paper leverages recent advancements in Gaussian Splatting and introduces a new method to learn generalizable human Gaussians that allows photorealistic and accurate view-rendering of a new human subject from a limited set of sparse views in a feed-forward manner. A pivotal innovation of our approach involves reformulating the learning of 3D Gaussian parameters into a regression process defined on the 2D UV space of a human template, which allows leveraging the strong geometry prior and the advantages of 2D convolutions. In addition, a multi-scaffold is proposed to effectively represent the offset details. Our method outperforms recent methods on both within-dataset generalization as well as cross-dataset generalization settings.
△ Less
Submitted 17 July, 2024;
originally announced July 2024.
-
Hamba: Single-view 3D Hand Reconstruction with Graph-guided Bi-Scanning Mamba
Authors:
Haoye Dong,
Aviral Chharia,
Wenbo Gou,
Francisco Vicente Carrasco,
Fernando De la Torre
Abstract:
3D Hand reconstruction from a single RGB image is challenging due to the articulated motion, self-occlusion, and interaction with objects. Existing SOTA methods employ attention-based transformers to learn the 3D hand pose and shape, but they fail to achieve robust and accurate performance due to insufficient modeling of joint spatial relations. To address this problem, we propose a novel graph-gu…
▽ More
3D Hand reconstruction from a single RGB image is challenging due to the articulated motion, self-occlusion, and interaction with objects. Existing SOTA methods employ attention-based transformers to learn the 3D hand pose and shape, but they fail to achieve robust and accurate performance due to insufficient modeling of joint spatial relations. To address this problem, we propose a novel graph-guided Mamba framework, named Hamba, which bridges graph learning and state space modeling. Our core idea is to reformulate Mamba's scanning into graph-guided bidirectional scanning for 3D reconstruction using a few effective tokens. This enables us to learn the joint relations and spatial sequences for enhancing the reconstruction performance. Specifically, we design a novel Graph-guided State Space (GSS) block that learns the graph-structured relations and spatial sequences of joints and uses 88.5% fewer tokens than attention-based methods. Additionally, we integrate the state space features and the global features using a fusion module. By utilizing the GSS block and the fusion module, Hamba effectively leverages the graph-guided state space modeling features and jointly considers global and local features to improve performance. Extensive experiments on several benchmarks and in-the-wild tests demonstrate that Hamba significantly outperforms existing SOTAs, achieving the PA-MPVPE of 5.3mm and F@15mm of 0.992 on FreiHAND. Hamba is currently Rank 1 in two challenging competition leaderboards on 3D hand reconstruction. The code will be available upon acceptance. [Website](https://humansensinglab.github.io/Hamba/).
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Graph Permutation Entropy: Extensions to the Continuous Case, A step towards Ordinal Deep Learning, and More
Authors:
Om Roy,
Avalon Campbell-Cousins,
John Stewart Fabila Carrasco,
Mario A Parra,
Javier Escudero
Abstract:
Nonlinear dynamics play an important role in the analysis of signals. A popular, readily interpretable nonlinear measure is Permutation Entropy. It has recently been extended for the analysis of graph signals, thus providing a framework for non-linear analysis of data sampled on irregular domains. Here, we introduce a continuous version of Permutation Entropy, extend it to the graph domain, and de…
▽ More
Nonlinear dynamics play an important role in the analysis of signals. A popular, readily interpretable nonlinear measure is Permutation Entropy. It has recently been extended for the analysis of graph signals, thus providing a framework for non-linear analysis of data sampled on irregular domains. Here, we introduce a continuous version of Permutation Entropy, extend it to the graph domain, and develop a ordinal activation function akin to the one of neural networks. This is a step towards Ordinal Deep Learning, a potentially effective and very recently posited concept. We also formally extend ordinal contrasts to the graph domain. Continuous versions of ordinal contrasts of length 3 are also introduced and their advantage is shown in experiments. We also integrate specific contrasts for the analysis of images and show that it generalizes well to the graph domain allowing a representation of images, represented as graph signals, in a plane similar to the entropy-complexity one. Applications to synthetic data, including fractal patterns and popular non-linear maps, and real-life MRI data show the validity of these novel extensions and potential benefits over the state of the art. By extending very recent concepts related to permutation entropy to the graph domain, we expect to accelerate the development of more graph-based entropy methods to enable nonlinear analysis of a broader kind of data and establishing relationships with emerging ideas in data science.
△ Less
Submitted 11 July, 2024; v1 submitted 10 July, 2024;
originally announced July 2024.
-
Taming 3DGS: High-Quality Radiance Fields with Limited Resources
Authors:
Saswat Subhajyoti Mallick,
Rahul Goel,
Bernhard Kerbl,
Francisco Vicente Carrasco,
Markus Steinberger,
Fernando De La Torre
Abstract:
3D Gaussian Splatting (3DGS) has transformed novel-view synthesis with its fast, interpretable, and high-fidelity rendering. However, its resource requirements limit its usability. Especially on constrained devices, training performance degrades quickly and often cannot complete due to excessive memory consumption of the model. The method converges with an indefinite number of Gaussians -- many of…
▽ More
3D Gaussian Splatting (3DGS) has transformed novel-view synthesis with its fast, interpretable, and high-fidelity rendering. However, its resource requirements limit its usability. Especially on constrained devices, training performance degrades quickly and often cannot complete due to excessive memory consumption of the model. The method converges with an indefinite number of Gaussians -- many of them redundant -- making rendering unnecessarily slow and preventing its usage in downstream tasks that expect fixed-size inputs. To address these issues, we tackle the challenges of training and rendering 3DGS models on a budget. We use a guided, purely constructive densification process that steers densification toward Gaussians that raise the reconstruction quality. Model size continuously increases in a controlled manner towards an exact budget, using score-based densification of Gaussians with training-time priors that measure their contribution. We further address training speed obstacles: following a careful analysis of 3DGS' original pipeline, we derive faster, numerically equivalent solutions for gradient computation and attribute updates, including an alternative parallelization for efficient backpropagation. We also propose quality-preserving approximations where suitable to reduce training time even further. Taken together, these enhancements yield a robust, scalable solution with reduced training times, lower compute and memory requirements, and high quality. Our evaluation shows that in a budgeted setting, we obtain competitive quality metrics with 3DGS while achieving a 4--5x reduction in both model size and training time. With more generous budgets, our measured quality surpasses theirs. These advances open the door for novel-view synthesis in constrained environments, e.g., mobile devices.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
A novel multivariate regression model for unbalanced binary data : a strong conjugacy under random effect approach
Authors:
Lizandra C. Fabio,
Vanessa Barros,
Cristian Villegas,
Jalmar M. F. Carrasco
Abstract:
In this paper, we deduce a new multivariate regression model designed to fit correlated binary data. The multivariate distribution is derived from a Bernoulli mixed model with a nonnormal random intercept on the marginal approach. The random effect distribution is assumed to be the generalized log-gamma (GLG) distribution by considering a particular parameter setting. The complement log-log functi…
▽ More
In this paper, we deduce a new multivariate regression model designed to fit correlated binary data. The multivariate distribution is derived from a Bernoulli mixed model with a nonnormal random intercept on the marginal approach. The random effect distribution is assumed to be the generalized log-gamma (GLG) distribution by considering a particular parameter setting. The complement log-log function is specified to lead to strong conjugacy between the response variable and random effect. The new discrete multivariate distribution, named MBerGLG distribution, has location and dispersion parameters. The MBerGLG distribution leads to the MBerGLG regression (MBerGLGR) model, providing an alternative approach to fitting both unbalanced and balanced correlated response binary data. Monte Carlo simulation studies show that its maximum likelihood estimators are unbiased, efficient, and consistent asymptotically. The randomized quantile residuals are performed to identify possible departures from the proposal model and the data and detect atypical subjects. Finally, two applications are presented in the data analysis section.
△ Less
Submitted 6 June, 2024;
originally announced June 2024.
-
The role of turbulence and winding in the development of large-scale, strong magnetic fields in long-lived remnants of binary neutron star mergers
Authors:
Ricard Aguilera-Miret,
Carlos Palenzuela,
Federico Carrasco,
Daniele Viganò
Abstract:
We perform a long and accurate Large-Eddy Simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution and the gradient subgrid-scale model allow us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-He…
▽ More
We perform a long and accurate Large-Eddy Simulation of a binary neutron star merger, following the newly formed remnant up to 110 milliseconds. The combination of high-order schemes, high-resolution and the gradient subgrid-scale model allow us to have among the highest effective resolutions ever achieved. Our results show that, although the magnetic fields are strongly amplified by the Kelvin-Helmholtz instability, they are coherent only over very short spatial scales until t \gtrsim 30 ms. Around that time, magnetic winding becomes more efficient leading to a linear growth of the toroidal component and slowly ordering the field to more axisymmetric, large scales. The poloidal component only starts to grow at small scales at much later times t \gtrsim 90 ms, in a way compatible with the magneto-rotational instability. No strong large-scale poloidal field or jet is produced in the timescales spanned by our simulation, although there is an helicoidal structure gradually developing at late times. We highlight that soon after the merger the topology is always strongly dominated by toroidal structures, with a complex distribution in the meridional plane and highly turbulent perturbations. Thus, starting with strong purely dipolar fields before the merger is largely inconsistent with the outcomes of a realistic evolution. Finally, we confirm the universality of the evolved topology, even when starting with very different magnetic fields confined to the outermost layers of each neutron star.
△ Less
Submitted 10 July, 2023;
originally announced July 2023.
-
Average conditioning of underdetermined polynomial systems
Authors:
Federico Carrasco
Abstract:
This article study the average conditioning for a random underdetermined polynomial system. The expected value of the moments of the condition number are compared to the moments of the condition number of random matrices. An expression for these moments is given by studying the kernel finding problem for random matrices. Furthermore, the second moment of the Frobenius condition number is computed.
This article study the average conditioning for a random underdetermined polynomial system. The expected value of the moments of the condition number are compared to the moments of the condition number of random matrices. An expression for these moments is given by studying the kernel finding problem for random matrices. Furthermore, the second moment of the Frobenius condition number is computed.
△ Less
Submitted 7 August, 2023; v1 submitted 10 May, 2023;
originally announced May 2023.
-
Towards Realistic Generative 3D Face Models
Authors:
Aashish Rai,
Hiresh Gupta,
Ayush Pandey,
Francisco Vicente Carrasco,
Shingo Jason Takagi,
Amaury Aubel,
Daeil Kim,
Aayush Prakash,
Fernando de la Torre
Abstract:
In recent years, there has been significant progress in 2D generative face models fueled by applications such as animation, synthetic data generation, and digital avatars. However, due to the absence of 3D information, these 2D models often struggle to accurately disentangle facial attributes like pose, expression, and illumination, limiting their editing capabilities. To address this limitation,…
▽ More
In recent years, there has been significant progress in 2D generative face models fueled by applications such as animation, synthetic data generation, and digital avatars. However, due to the absence of 3D information, these 2D models often struggle to accurately disentangle facial attributes like pose, expression, and illumination, limiting their editing capabilities. To address this limitation, this paper proposes a 3D controllable generative face model to produce high-quality albedo and precise 3D shape leveraging existing 2D generative models. By combining 2D face generative models with semantic face manipulation, this method enables editing of detailed 3D rendered faces. The proposed framework utilizes an alternating descent optimization approach over shape and albedo. Differentiable rendering is used to train high-quality shapes and albedo without 3D supervision. Moreover, this approach outperforms the state-of-the-art (SOTA) methods in the well-known NoW benchmark for shape reconstruction. It also outperforms the SOTA reconstruction models in recovering rendered faces' identities across novel poses by an average of 10%. Additionally, the paper demonstrates direct control of expressions in 3D faces by exploiting latent space leading to text-based editing of 3D faces.
△ Less
Submitted 26 October, 2023; v1 submitted 24 April, 2023;
originally announced April 2023.
-
Relativistic force-free models of the thermal X-ray emission in millisecond pulsars observed by NICER
Authors:
Federico Carrasco,
Joaquin Pelle,
Oscar Reula,
Daniele Viganò,
Carlos Palenzuela
Abstract:
Several important properties of rotation-powered millisecond pulsars (MSPs), such as their mass-radius ratio, equation of state and magnetic field topology, can be inferred from precise observations and modelling of their X-ray light curves. In the present study, we model the thermal X-ray signals originated in MSPs, all the way from numerically solving the surrounding magnetospheres up to the ray…
▽ More
Several important properties of rotation-powered millisecond pulsars (MSPs), such as their mass-radius ratio, equation of state and magnetic field topology, can be inferred from precise observations and modelling of their X-ray light curves. In the present study, we model the thermal X-ray signals originated in MSPs, all the way from numerically solving the surrounding magnetospheres up to the ray tracing of the emitted photons and the final computation of their light curves and spectra. Our modelled X-ray signals are then compared against a set of very accurate NICER observations of four target pulsars: PSR J0437-4715, PSR J1231-1411, PSR J2124-3358 and PSR J0030+0451. We find very good simultaneous fits for the light curve and spectral distribution in all these pulsars. The magnetosphere is solved by performing general relativistic force-free simulations of a rotating neutron star (NS) endowed with a simple centered dipolar magnetic field, for many different stellar compactness and pulsar misalignments. From these solutions, we derive an emissivity map over the surface of the star, which is based on the electric currents in the magnetosphere. In particular, the emission regions (ERs) are determined in this model by spacelike four-currents that reach the NS. We show that this assumption, together with the inclusion of the gravitational curvature on the force-free simulations, lead to non-standard ERs facing the closed-zone of the pulsar, in addition to other ERs within the polar caps. The combined X-ray signals from these two kinds of ERs (both antipodal) allow to approximate the non-trivial interpulses found in all the target MSPs light curves.
△ Less
Submitted 23 September, 2022;
originally announced September 2022.
-
Skylight: a new code for general-relativistic ray-tracing and radiative transfer in arbitrary spacetimes
Authors:
Joaquín Pelle,
Oscar Reula,
Federico Carrasco,
Carlos Bederian
Abstract:
To reproduce the observed spectra and light curves originated in the neighborhood of compact objects requires accurate relativistic ray-tracing codes. In this work, we present Skylight, a new numerical code for general-relativistic ray tracing and radiative transfer in arbitrary space-time geometries and coordinate systems. The code is capable of producing images, spectra, and light curves from as…
▽ More
To reproduce the observed spectra and light curves originated in the neighborhood of compact objects requires accurate relativistic ray-tracing codes. In this work, we present Skylight, a new numerical code for general-relativistic ray tracing and radiative transfer in arbitrary space-time geometries and coordinate systems. The code is capable of producing images, spectra, and light curves from astrophysical models of compact objects as seen by distant observers. We incorporate two different schemes, namely Monte Carlo radiative transfer, integrating geodesics from the astrophysical region to distant observers, and camera techniques with backwards integration from the observer to the emission region. The code is validated by successfully passing several test cases, among them: thin accretion disks and neutron star hot spot emission.
△ Less
Submitted 13 June, 2022;
originally announced June 2022.
-
Turbulent magnetic field amplification in binary neutron star mergers
Authors:
C. Palenzuela,
R. Aguilera-Miret,
F. Carrasco,
R. Ciolfi,
J. V. Kalinani,
W. Kastaun,
B. Miñano,
D. Viganò
Abstract:
Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insu…
▽ More
Magnetic fields are expected to play a key role in the dynamics and the ejection mechanisms that accompany the merger of two neutron stars. General relativistic magnetohydrodynamic (MHD) simulations offer a unique opportunity to unravel the details of the ongoing physical processes. Nevertheless, current numerical studies are severely limited by the fact that any affordable resolution remains insufficient to fully capture the small-scale dynamo, initially triggered by the Kelvin-Helmholtz instability, and later sourced by several MHD processes involving differential rotation. Here, we alleviate this limitation by using explicit large-eddy simulations, a technique where the unresolved dynamics occurring at the sub-grid scales (SGS) is modeled by extra terms, which are functions of the resolved fields and their derivatives. The combination of high-order numerical schemes, high resolutions, and the gradient SGS model allow us to capture the small-scale dynamos produced during the binary neutron star mergers. Here we follow the first 50 milliseconds after the merger and, for the first time, we find numerical convergence on the magnetic field amplification, in terms of integrated energy and spectral distribution over spatial scales. We also find that the average intensity of the magnetic field in the remnant saturates at $\sim 10^{16}$~G around $5$~ms after the merger. After $20-30$~ms, both toroidal and poloidal magnetic field components grow continuously, fed by the winding mechanism that provides a slow inverse cascade. We find no clear hints for magneto-rotational instabilities, and no significant impact of the magnetic field on the redistribution of angular momentum in the remnant in our simulations, probably due to the very turbulent and dynamical topology of the magnetic field at all stages, with small-scale components largely dominating over the large-scale ones.
△ Less
Submitted 3 July, 2022; v1 submitted 15 December, 2021;
originally announced December 2021.
-
Gamma-radiation sky maps from compact binaries
Authors:
Néstor Ortiz,
Federico Carrasco,
Stephen R. Green,
Luis Lehner,
Steven L. Liebling,
John Ryan Westernacher-Schneider
Abstract:
We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our…
▽ More
We study sky maps and light curves of gamma-ray emission from neutron stars in compact binaries, and in isolation. We briefly review some gamma-ray emission models, and reproduce sky maps from a standard isolated pulsar in the Separatrix Layer model. We consider isolated pulsars with several variations of a dipole magnetic field, including superpositions, and predict their gamma-ray emission. Our results provide new heuristics on what can and cannot be inferred about the magnetic field configuration of pulsars from high-energy observations. We find that typical double-peak light curves can be produced by pulsars with significant multipole structure beyond a single dipole. For binary systems, we also present a simple approximation that is useful for rapid explorations of binary magnetic field structure. Finally, we predict the gamma-ray emission pattern from a compact black hole-neutron star binary moments before merger by applying the Separatrix Layer model to data simulated in full general relativity; we find that face-on observers receive little emission, equatorial observers see one broad peak, and more generic observers typically see two peaks.
△ Less
Submitted 17 February, 2022; v1 submitted 14 July, 2021;
originally announced July 2021.
-
Magnetospheres of black hole-neutron star binaries
Authors:
Federico Carrasco,
Masaru Shibata,
Oscar Reula
Abstract:
We perform force-free simulations for a neutron star orbiting a black hole, aiming at clarifying the main magnetosphere properties of such binaries towards their innermost stable circular orbits. Several configurations are explored, varying the orbital separation, the individual spins and misalignment angle among the magnetic and orbital axes. We find significant electromagnetic luminosities,…
▽ More
We perform force-free simulations for a neutron star orbiting a black hole, aiming at clarifying the main magnetosphere properties of such binaries towards their innermost stable circular orbits. Several configurations are explored, varying the orbital separation, the individual spins and misalignment angle among the magnetic and orbital axes. We find significant electromagnetic luminosities, $L\sim 10^{42-46} \, [B_{\rm pole}/ 10^{12}{\rm G}]^2 \, {\rm erg/s}$ (depending on the specific setting), primarily powered by the orbital kinetic energy, being about one order of magnitude higher than those expected from unipolar induction. The systems typically develop current sheets that extend to long distances following a spiral arm structure. The intense curvature of the black hole produces extreme bending on a particular set of magnetic field lines as it moves along the orbit, leading to magnetic reconnections in the vicinity of the horizon. For the most symmetric scenario (aligned cases), these reconnection events can release large-scale plasmoids that carry the majority of the Poynting fluxes. On the other hand, for misaligned cases, a larger fraction of the luminosity is instead carried outwards by large-amplitude Alfv{é}n waves disturbances. We estimate possible precursor electromagnetic emissions based on our numerical solutions, finding radio signals as the most promising candidates to be detectable within distances of $\lesssim 200$\,Mpc by forthcoming facilities like the Square Kilometer Array.
△ Less
Submitted 16 June, 2021;
originally announced June 2021.
-
Diagnostic tools for a multivariate negative binomial model for fitting correlated data with overdispersion
Authors:
Lizandra Castilho Fabio,
Cristian Villegas,
Jalmar M. F. Carrasco,
Mário de Castro
Abstract:
We focus on the development of diagnostic tools and an R package called MNB for a multivariate negative binomial (MNB) regression model for detecting atypical and influential subjects. The MNB model is deduced from a Poisson mixed model in which the random intercept follows the generalized log-gamma (GLG) distribution. The MNB model for correlated count data leads to an MNB regression model that i…
▽ More
We focus on the development of diagnostic tools and an R package called MNB for a multivariate negative binomial (MNB) regression model for detecting atypical and influential subjects. The MNB model is deduced from a Poisson mixed model in which the random intercept follows the generalized log-gamma (GLG) distribution. The MNB model for correlated count data leads to an MNB regression model that inherits the features of a hierarchical model to accommodate the intraclass correlation and the occurrence of overdispersion simultaneously. The asymptotic consistency of the dispersion parameter estimator depends on the asymmetry of the GLG distribution. Inferential procedures for the MNB regression model are simple, although it can provide inconsistent estimates of the asymptotic variance when the correlation structure is misspecified. We propose the randomized quantile residual for checking the adequacy of the multivariate model, and derive global and local influence measures from the multivariate model to assess influential subjects. Finally, two applications are presented in the data analysis section. The code for installing the MNB package and the code used in the two examples is exhibited in the Appendix.
△ Less
Submitted 15 February, 2021;
originally announced February 2021.
-
Turbulent magnetic-field amplification in the first 10 milliseconds after a binary neutron star merger: comparing high-resolution and large eddy simulations
Authors:
Ricard Aguilera-Miret,
Daniele Viganò,
Federico Carrasco,
Borja Miñano,
Carlos Palenzuela
Abstract:
The detection of binary neutron star mergers represents one of the most important and complex astrophysical discoveries of the recent years. One of the unclear aspects of the problem is the turbulent magnetic field amplification, initially triggered by the Kelvin-Helmholtz instability at much smaller scales than any reachable numerical resolution nowadays. Here we present numerical simulations of…
▽ More
The detection of binary neutron star mergers represents one of the most important and complex astrophysical discoveries of the recent years. One of the unclear aspects of the problem is the turbulent magnetic field amplification, initially triggered by the Kelvin-Helmholtz instability at much smaller scales than any reachable numerical resolution nowadays. Here we present numerical simulations of the first ten milliseconds of a binary neutron star merger. First, we confirm in detail how the simulated amplification depends on the numerical resolution and is distributed on a broad range of scales, as expected from turbulent MHD theory. We find that an initial large-scale magnetic field of $10^{11}\,$G inside each star is amplified in the remnant to root-mean-square values above $10^{16}\,$G within the first $5$ milliseconds for our highest-resolution run. Then, we run large eddy simulations, exploring the performance of the subgrid-scale gradient model, already tested successfully in previous turbulent box simulations. We show that the addition of this model is especially important in the induction equation, since it leads to an amplification of the magnetic field comparable to a higher-resolution run, but with a greatly reduced computational cost. In the first 10 milliseconds, there is no clear hint for an ordered, large-scale magnetic field, which should indeed occur in longer timescales through magnetic winding and the magneto-rotational instability.
△ Less
Submitted 14 September, 2020;
originally announced September 2020.
-
GRMHD large eddy simulations with gradient subgrid-scale model
Authors:
Daniele Viganò,
Ricard Aguilera-Miret,
Federico Carrasco,
Borja Miñano,
Carlos Palenzuela
Abstract:
The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations w…
▽ More
The detection of binary neutron star mergers represents one of the most important astrophysical discoveries of the recent years. Due to the extreme matter and gravity conditions and the rich dynamics developed, it becomes a tremendous challenge to accurately simulate numerically all the scales present during the collision. Here we present how to study such systems by using large eddy simulations with a self-consistent subgrid-scale gradient model, that we generalized to the special relativistic case in a previous work and now extend to the general relativistic case. Adapted from nonrelativistic scenarios, the so-called gradient model allows to capture part of the effects of the hidden dynamics on the resolved scales, by means of a physically-agnostic, mathematically-based Taylor expansion of the nonlinear terms in the conservative evolution equations' fluxes. We assess the validity of this approach in bounding-box simulations of the magnetic Kelvin-Helmholtz instability. Several resolutions and a broad range of scenarios are considered in order to carefully test the performance of the model under three crucial aspects: (i) highly curved backgrounds, (ii) jumps on the fluid density profiles and (iii) strong shocks. The results suggest our extension of the gradient subgrid-scale model to general relativistic magnetohydrodynamics is a promising approach for studying binary neutron stars mergers, and potentially to other relevant astrophysical scenarios.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
Magnetosphere of an orbiting neutron star
Authors:
Federico Carrasco,
Masaru Shibata
Abstract:
We conduct force-free simulations of a single neutron star undergoing orbital motion in flat spacetime, mimicking the trajectory of the star about the center of mass on a compact binary system. Our attention is focused on the kinetic energy being extracted from the orbit by the acceleration of the magnetic dipole moment of the neutron star, and particularly, on how this energy gets distributed wit…
▽ More
We conduct force-free simulations of a single neutron star undergoing orbital motion in flat spacetime, mimicking the trajectory of the star about the center of mass on a compact binary system. Our attention is focused on the kinetic energy being extracted from the orbit by the acceleration of the magnetic dipole moment of the neutron star, and particularly, on how this energy gets distributed within its surrounding magnetosphere. A detailed study of the resulting magnetospheric configurations in our setting is presented, incorporating as well the effects due to neutron star spin and the misalignment of the magnetic and orbital axes. We find many features resembling those of pulsar magnetospheres for the orbiting neutron star --even in the absence of spin--, being of particular interest the development of a spiral current sheet that extends beyond the light cylinder. Then, we use recent advances in pulsar theory to estimate electromagnetic emissions produced at the reconnection regions of such current sheets.
△ Less
Submitted 19 February, 2020; v1 submitted 13 January, 2020;
originally announced January 2020.
-
Gradient sub-grid-scale model for relativistic MHD Large Eddy Simulations
Authors:
Federico Carrasco,
Daniele Viganò,
Carlos Palenzuela
Abstract:
MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at…
▽ More
MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at play, which is instead effectively modeled by a set of additional terms in the evolution equations, dubbed as sub-grid-scale model. Here we extend such approach, commonly used in non-relativistic/non-magnetic/incompressible fluid dynamics, applying the so-called gradient model to a general set of balance-law equations, that includes the relevant case in which a non-trivial inversion of conserved to primitive fields is needed. In particular, we focus on the relativistic compressible ideal MHD scenario, providing for the first time (and for any equation of state) all the additional sub-grid-scale terms. As an application, we consider box simulations of the relativistic Kelvin-Helmholtz instability, which is also the first mechanism responsible for the magnetic field amplification in binary neutron star mergers and cannot yet be fully captured by the finest-grid and longest simulations available. The performance of our model is numerically assessed by comparing it to the residuals arising from the filtering of high-resolution simulations. We find that the model can fit very well those residuals from resolutions a few times higher. Although the application shown here explicitly considers the Minkowski metric, it can be directly extended to general relativity, thus settling the basis to implement the gradient sub-grid model in a GRMHD binary merger. Our results suggest that this approach will be potentially able to unveil much better the small-scale dynamics achievable in full GRMHD simulations.
△ Less
Submitted 23 March, 2020; v1 submitted 4 August, 2019;
originally announced August 2019.
-
Astrophysical jets from boosted compact objects
Authors:
Ramiro Cayuso,
Federico Carrasco,
Barbara Sbarato,
Oscar Reula
Abstract:
We perform full 3D numerical simulations of compact objects, such as black holes or neutron stars, boosted through an ambient force-free plasma that posses a uniform magnetization. We study jet formation and energy extraction from the resulting stationary late time solutions.The implementation of appropriate boundary conditions has allowed us to explore a wide range of boost velocities, finding th…
▽ More
We perform full 3D numerical simulations of compact objects, such as black holes or neutron stars, boosted through an ambient force-free plasma that posses a uniform magnetization. We study jet formation and energy extraction from the resulting stationary late time solutions.The implementation of appropriate boundary conditions has allowed us to explore a wide range of boost velocities, finding the jet power scales as $γv^2$ (being $γ$ the Lorentz factor). We also explore other parameters of the problem like the orientation of the motion respect to the asymptotic magnetic field or the inclusion of black hole spin. Additionally, by comparing a black hole with a perfectly conducting sphere in flat spacetime, we manage to disentangle curvature effects from those produced by the perfect conducting surface. It is shown that when the stellar compactness is increased these two effects act in combination, further enhancing the luminosity produced by the neutron star.
△ Less
Submitted 1 May, 2019;
originally announced May 2019.
-
Triggering magnetar outbursts in 3D force-free simulations
Authors:
Federico Carrasco,
Daniele Viganò,
Carlos Palenzuela,
Jose A. Pons
Abstract:
In this letter, we present the first 3D force-free general relativity simulations of the magnetosphere dynamics related to the magnetar outburst/flare phenomenology. Starting from an initial dipole configuration, we adiabatically increase the helicity by twisting the footprints of a spot on the stellar surface and follow the succession of quasi-equilibrium states until a critical twist is reached.…
▽ More
In this letter, we present the first 3D force-free general relativity simulations of the magnetosphere dynamics related to the magnetar outburst/flare phenomenology. Starting from an initial dipole configuration, we adiabatically increase the helicity by twisting the footprints of a spot on the stellar surface and follow the succession of quasi-equilibrium states until a critical twist is reached. Twisting beyond that point triggers instabilities that results in the rapid expansion of magnetic field lines, followed by reconnection, as observed in previous axi-symmetric simulations. If the injection of magnetic helicity goes on, the process is recurrent, periodically releasing a similar amount of energy, of the order of a few % of the total magnetic energy. From our current distribution, we estimate the local temperature assuming that dissipation occurs mainly in the highly resistive outermost layer of the neutron star. We find that the temperature smoothly increases with injected twist, being larger for spots located in the tropical regions than in polar regions, and rather independent of their sizes. After the injection of helicity ceases, the magnetosphere relaxes to a new stable state, in which the persistent currents maintain the footprints area slightly hotter than before the onset of the instability.
△ Less
Submitted 25 January, 2019;
originally announced January 2019.
-
A Simflowny-based high-performance 3D code for the generalized induction equation
Authors:
Daniele Viganò,
David Martínez-Gómez,
José A. Pons,
Carlos Palenzuela,
Federico Carrasco,
Borja Miñano,
Antoni Arbona,
Carles Bona,
Joan Massó
Abstract:
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage…
▽ More
In the interior of neutron stars, the induction equation regulates the long-term evolution of the magnetic fields by means of resistivity, Hall dynamics and ambipolar diffusion. Despite the apparent simplicity and compactness of the equation, the dynamics it describes is not trivial and its understanding relies on accurate numerical simulations. While a few works in 2D have reached a mature stage and a consensus on the general dynamics at least for some simple initial data, only few attempts have been performed in 3D, due to the computational costs and the need for a proper numerical treatment of the intrinsic non-linearity of the equation. Here, we carefully analyze the general induction equation, studying its characteristic structure, and we present a new Cartesian 3D code, generated by the user-friendly, publicly available {\em Simflowny} platform. The code uses high-order numerical schemes for the time and spatial discretization, and relies on the highly-scalable {\em SAMRAI} architecture for the adaptive mesh refinement. We present the application of the code to several benchmark tests, showing the high order of convergence and accuracy achieved and the capabilities in terms of magnetic shock resolution and three-dimensionality. This paper paves the way for the applications to a realistic, 3D long-term evolution of neutron stars interior and, possibly, of other astrophysical sources.
△ Less
Submitted 11 December, 2018; v1 submitted 20 November, 2018;
originally announced November 2018.
-
Pulsar magnetospheres in General Relativity
Authors:
Federico Carrasco,
Carlos Palenzuela,
Oscar Reula
Abstract:
The main contribution to the pulsar power can be calculated by assuming a rotating magnetically-dominated magnetosphere described by the force-free approximation. Although this simple model has been used thoroughly to study pulsar magnetospheres in the flat spacetime regime, only few works have considered the relativistic corrections introduced by the curvature and frame-dragging effects induced b…
▽ More
The main contribution to the pulsar power can be calculated by assuming a rotating magnetically-dominated magnetosphere described by the force-free approximation. Although this simple model has been used thoroughly to study pulsar magnetospheres in the flat spacetime regime, only few works have considered the relativistic corrections introduced by the curvature and frame-dragging effects induced by a rotating neutron star. Here we revisit the problem and describe pulsar magnetospheres within full General Relativity, quantifying the corrections as a function of the angular velocity, the compactness of the star and the misalignment angle between the spin and the magnetic dipole. We provide analytical expressions for the pulsar luminosity by fitting our numerical results. Finally, we also analyze the effect of the relativistic corrections on the braking index, which indicates a slight increment in its value.
△ Less
Submitted 15 June, 2018; v1 submitted 10 May, 2018;
originally announced May 2018.
-
A novel scheme for simulating the force-free equations: boundary conditions and the evolution of solutions towards stationarity
Authors:
Federico Carrasco,
Oscar Reula
Abstract:
Force-Free Electrodynamics (FFE) describes a particular regime of magnetically dominated relativistic plasmas, which arises on several astrophysical scenarios of interest such as pulsars or active galactic nuclei. In this article, we present a full 3D numerical implementation of the FFE evolution around a Kerr black hole. The novelty of our approach is three-folded: i) We use the "multi-block" tec…
▽ More
Force-Free Electrodynamics (FFE) describes a particular regime of magnetically dominated relativistic plasmas, which arises on several astrophysical scenarios of interest such as pulsars or active galactic nuclei. In this article, we present a full 3D numerical implementation of the FFE evolution around a Kerr black hole. The novelty of our approach is three-folded: i) We use the "multi-block" technique to represent a domain with $S^2 \times \mathbb{R}^{+}$ topology within a stable finite-differences scheme. ii) We employ as evolution equations those arising from a covariant hyperbolization of the FFE system. iii) We implement stable and constraint-preserving boundary conditions to represent an outer region given by a uniform magnetic field aligned or misaligned respect to the symmetry axis.
We find stationary jet solutions which reach equilibrium --through boundary conditions-- with the outer numerical surface. This is so, even when the outer boundary is located very close to the central region (i.e. $r_{out}\sim 10M $). These numerical solutions reproduce most of the known results for analogue astrophysical settings.
△ Less
Submitted 29 March, 2017;
originally announced March 2017.
-
Covariant Hyperbolization of Force-free Electrodynamics
Authors:
Federico Carrasco,
Oscar Reula
Abstract:
Force-Free Flectrodynamics (FFE) is a non-linear system of equations modeling the evolution of the electromagnetic field, in the presence of a magnetically dominated relativistic plasma. This configuration arises on several astrophysical scenarios, which represent exciting laboratories to understand physics in extreme regimes. We show that this system, when restricted to the correct constraint sub…
▽ More
Force-Free Flectrodynamics (FFE) is a non-linear system of equations modeling the evolution of the electromagnetic field, in the presence of a magnetically dominated relativistic plasma. This configuration arises on several astrophysical scenarios, which represent exciting laboratories to understand physics in extreme regimes. We show that this system, when restricted to the correct constraint submanifold, is symmetric hyperbolic. In numerical applications is not feasible to keep the system in that submanifold, and so, it is necessary to analyze its structure first in the tangent space of that submanifold and then in a whole neighborhood of it. As already shown by Pfeiffer, a direct (or naive) formulation of this system (in the whole tangent space) results in a weakly hyperbolic system of evolution equations for which well-possednes for the initial value formulation does not follows.
Using the generalized symmetric hyperbolic formalism due to Geroch, we introduce here a covariant hyperbolization for the FFE system. In fact, in analogy to the usual Maxwell case, a complete family of hyperbolizers is found, both for the restricted system on the constraint submanifold as well as for a suitably extended system defined in a whole neighborhood of it. A particular symmetrizer among the family is then used to write down the pertaining evolution equations, in a generic ($3+1$)-decomposition on a background spacetime. Interestingly, it turns out that for a particular choice of the lapse and shift functions of the foliation, our symmetrized system reduces to the one found by Pfeiffer. Finally we analyze the characteristic structure of the resulting evolution system.
△ Less
Submitted 13 February, 2016; v1 submitted 4 February, 2016;
originally announced February 2016.
-
Nonlinear electrodynamics as a symmetric hyperbolic system
Authors:
Fernando Abalos,
Federico Carrasco,
Érico Goulart,
Oscar Reula
Abstract:
Nonlinear theories generalizing Maxwell's electromagnetism and arising from a Lagrangian formalism have dispersion relations in which propagation planes factor into null planes corresponding to two effective metrics which depend on the point-wise values of the electromagnetic field. These effective Lorentzian metrics share the null (generically two) directions of the electromagnetic field. We show…
▽ More
Nonlinear theories generalizing Maxwell's electromagnetism and arising from a Lagrangian formalism have dispersion relations in which propagation planes factor into null planes corresponding to two effective metrics which depend on the point-wise values of the electromagnetic field. These effective Lorentzian metrics share the null (generically two) directions of the electromagnetic field. We show that, the theory is symmetric hyperbolic if and only if the cones these metrics give rise to have a non-empty intersection. Namely that there exist families of symmetrizers in the sense of Geroch which are positive definite for all covectors in the interior of the cones intersection. Thus, for these theories, the initial value problem is well-posed. We illustrate the power of this approach with several nonlinear models of physical interest such as Born-Infeld, Gauss-Bonnet and Euler-Heisenberg.
△ Less
Submitted 13 July, 2015; v1 submitted 8 July, 2015;
originally announced July 2015.
-
A Holographic Path to the Turbulent Side of Gravity
Authors:
Stephen R. Green,
Federico Carrasco,
Luis Lehner
Abstract:
We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black b…
▽ More
We study the dynamics of a 2+1 dimensional relativistic viscous conformal fluid in Minkowski spacetime. Such fluid solutions arise as duals, under the "gravity/fluid correspondence", to 3+1 dimensional asymptotically anti-de Sitter (AAdS) black brane solutions to the Einstein equation. We examine stability properties of shear flows, which correspond to hydrodynamic quasinormal modes of the black brane. We find that, for sufficiently high Reynolds number, the solution undergoes an inverse turbulent cascade to long wavelength modes. We then map this fluid solution, via the gravity/fluid duality, into a bulk metric. This suggests a new and interesting feature of the behavior of perturbed AAdS black holes and black branes, which is not readily captured by a standard quasinormal mode analysis. Namely, for sufficiently large perturbed black objects (with long-lived quasinormal modes), nonlinear effects transfer energy from short to long wavelength modes via a turbulent cascade within the metric perturbation. As long wavelength modes have slower decay, this lengthens the overall lifetime of the perturbation. We also discuss various implications of this behavior, including expectations for higher dimensions, and the possibility of predicting turbulence in more general gravitational scenarios.
△ Less
Submitted 7 October, 2013; v1 submitted 30 September, 2013;
originally announced September 2013.
-
Multi-solitons with vector mesons on the two-sphere
Authors:
Federico Carrasco,
Oscar Reula
Abstract:
Recent studies have suggested a strong connection between the static solutions of the 3D Skyrme model and those corresponding to its low-dimensional analog (baby-Skyrme model) on a two-sphere. We have found almost identical solutions considering an alternative two-dimensional model in which a vector meson field is introduced and coupled to the system, instead of the usual Skyrme term. It has been…
▽ More
Recent studies have suggested a strong connection between the static solutions of the 3D Skyrme model and those corresponding to its low-dimensional analog (baby-Skyrme model) on a two-sphere. We have found almost identical solutions considering an alternative two-dimensional model in which a vector meson field is introduced and coupled to the system, instead of the usual Skyrme term. It has been known that including this vector meson field in three dimensions stabilizes the non-linear sigma model without the need of a term quartic on derivatives of the pion fields (Skyrme term). In the present work, we have numerically searched for static multi-solitonic solutions of this alternative stabilization, for the case in which the base-space is a two-sphere. Moreover, we analyze the stability of these solutions under small perturbations in a fully dynamical setting. We have also considered the inclusion of a particular potential term into the Lagrangian, and explored the low/high-density phases of solitons for different ranges of the parameter space, achieving solitons localized enough that allow us to compare with planar (two-dimensional) studies.
△ Less
Submitted 16 July, 2013;
originally announced July 2013.
-
Errors-in-variables beta regression models
Authors:
Jalmar M. F. Carrasco,
Silvia L. P. Ferrari,
Reinaldo B. Arellano-Valle
Abstract:
Beta regression models provide an adequate approach for modeling continuous outcomes limited to the interval (0,1). This paper deals with an extension of beta regression models that allow for explanatory variables to be measured with error. The structural approach, in which the covariates measured with error are assumed to be random variables, is employed. Three estimation methods are presented, n…
▽ More
Beta regression models provide an adequate approach for modeling continuous outcomes limited to the interval (0,1). This paper deals with an extension of beta regression models that allow for explanatory variables to be measured with error. The structural approach, in which the covariates measured with error are assumed to be random variables, is employed. Three estimation methods are presented, namely maximum likelihood, maximum pseudo-likelihood and regression calibration. Monte Carlo simulations are used to evaluate the performance of the proposed estimators and the naïve estimator. Also, a residual analysis for beta regression models with measurement errors is proposed. The results are illustrated in a real data set.
△ Less
Submitted 10 April, 2013; v1 submitted 4 December, 2012;
originally announced December 2012.
-
Turbulent flows for relativistic conformal fluids in 2+1 dimensions
Authors:
Federico Carrasco,
Luis Lehner,
Robert C. Myers,
Oscar Reula,
Ajay Singh
Abstract:
We demonstrate that relativistic conformal hydrodynamics in 2+1 dimensions displays a turbulent behaviour which cascades energy to longer wavelengths on both flat and spherical manifolds. Our motivation for this study is to understand the implications for gravitational solutions through the AdS/CFT correspondence. The observed behaviour implies gravitational perturbations of the corresponding blac…
▽ More
We demonstrate that relativistic conformal hydrodynamics in 2+1 dimensions displays a turbulent behaviour which cascades energy to longer wavelengths on both flat and spherical manifolds. Our motivation for this study is to understand the implications for gravitational solutions through the AdS/CFT correspondence. The observed behaviour implies gravitational perturbations of the corresponding black brane/black hole spacetimes (for sufficiently large scales/temperatures) will display a similar cascade towards longer wavelengths.
△ Less
Submitted 28 October, 2012; v1 submitted 24 October, 2012;
originally announced October 2012.
-
A New Generalized Kumaraswamy Distribution
Authors:
Jalmar M. F. Carrasco,
Silvia L. P. Ferrari,
Gauss M. Cordeiro
Abstract:
A new five-parameter continuous distribution which generalizes the Kumaraswamy and the beta distributions as well as some other well-known distributions is proposed and studied. The model has as special cases new four- and three-parameter distributions on the standard unit interval. Moments, mean deviations, Rényi's entropy and the moments of order statistics are obtained for the new generalized K…
▽ More
A new five-parameter continuous distribution which generalizes the Kumaraswamy and the beta distributions as well as some other well-known distributions is proposed and studied. The model has as special cases new four- and three-parameter distributions on the standard unit interval. Moments, mean deviations, Rényi's entropy and the moments of order statistics are obtained for the new generalized Kumaraswamy distribution. The score function is given and estimation is performed by maximum likelihood. Hypothesis testing is also discussed. A data set is used to illustrate an application of the proposed distribution.
△ Less
Submitted 6 April, 2010;
originally announced April 2010.