-
Ncorpi$\mathcal{O}$N : A $\mathcal{O}(N)$ software for N-body integration in collisional and fragmenting systems
Authors:
Jérémy Couturier,
Alice C. Quillen,
Miki Nakajima
Abstract:
Ncorpi$\mathcal{O}$N is a $N$-body software developed for the time-efficient integration of collisional and fragmenting systems of planetesimals or moonlets orbiting a central mass. It features a fragmentation model, based on crater scaling and ejecta models, able to realistically simulate a violent impact. The user of Ncorpi$\mathcal{O}$N can choose between four different built-in modules to comp…
▽ More
Ncorpi$\mathcal{O}$N is a $N$-body software developed for the time-efficient integration of collisional and fragmenting systems of planetesimals or moonlets orbiting a central mass. It features a fragmentation model, based on crater scaling and ejecta models, able to realistically simulate a violent impact. The user of Ncorpi$\mathcal{O}$N can choose between four different built-in modules to compute self-gravity and detect collisions. One of these makes use of a mesh-based algorithm to treat mutual interactions in $\mathcal{O}(N)$ time. Another module, much more efficient than the standard Barnes-Hut tree code, is a $\mathcal{O}(N)$ tree-based algorithm called FalcON. It relies on fast multipole expansion for gravity computation and we adapted it to collision detection as well. Computation time is reduced by building the tree structure using a three-dimensional Hilbert curve. For the same precision in mutual gravity computation, Ncorpi$\mathcal{O}$N is found to be up to 25 times faster than the famous software REBOUND. Ncorpi$\mathcal{O}$N is written entirely in the C language and only needs a C compiler to run. A python add-on, that requires only basic python libraries, produces animations of the simulations from the output files. The name Ncorpi$\mathcal{O}$N, reminding of a scorpion, comes from the French $N$-corps, meaning $N$-body, and from the mathematical notation $\mathcal{O}(N)$, due to the running time of the software being almost linear in the total number $N$ of moonlets. Ncorpi$\mathcal{O}$N is designed for the study of accreting or fragmenting disks of planetesimal or moonlets. It detects collisions and computes mutual gravity faster than REBOUND, and unlike other $N$-body integrators, it can resolve a collision by fragmentation. The fast multipole expansions are implemented up to order six to allow for a high precision in mutual gravity computation.
△ Less
Submitted 26 March, 2024; v1 submitted 31 October, 2023;
originally announced October 2023.
-
Wind erosion and transport on planetesimals
Authors:
Alice C. Quillen,
Stephen Luniewski,
Adam E. Rubinstein,
Jeremy Couturier,
Rachel Glade,
Miki Nakajima
Abstract:
We consider the possibility that aeolian (wind blown) processes occur on small, 1 to 100~km diameter, planetesimals when they were embedded in the protosolar nebula. Drag from a headwind within a protostellar disk is sufficiently large to loft cm and smaller sized particles off the surface of a 10 km diameter asteroid in the inner solar system (at a few AU), and micron sized particles off the surf…
▽ More
We consider the possibility that aeolian (wind blown) processes occur on small, 1 to 100~km diameter, planetesimals when they were embedded in the protosolar nebula. Drag from a headwind within a protostellar disk is sufficiently large to loft cm and smaller sized particles off the surface of a 10 km diameter asteroid in the inner solar system (at a few AU), and micron sized particles off the surface of a 10 km diameter object in the Transneptunian region. The headwind is sufficiently strong to overcome surface cohesion in the inner solar system, but not in the outer solar system. However, in the outer solar system, surface particles can be redistributed or escape due to impacts from particles that are in the protosolar disk's wind. Based on scaling crater ejecta, we estimate that impacts from particles in the headwind will lead to erosion of mass rather than accretion for planetesimals below about 6 km in diameter. The erosion limit is independent of material strength but proportional to the wind velocity. We explore the sensitivity of splash particle trajectories to particle size, headwind velocity and Reynolds number. Winds from a protostellar disk could account for Kuiper Belt Object (486958) Arrokoth's smooth undulating terrain but only during an epoch of high particle flux and low wind velocity. These conditions could have been present during and just after coalescence of Arrokoth's building blocks.
△ Less
Submitted 10 January, 2024; v1 submitted 14 October, 2023;
originally announced October 2023.
-
Subsurface pulse, crater and ejecta asymmetry from oblique impacts into granular media
Authors:
Bingcheng Suo,
A. C. Quillen,
Max Neiderbach,
Luke O'Brient,
Abobakar Sediq Miakhel,
Nathan Skerrett,
Jérémy Couturier,
Victor Lherm,
Jiaxin Wang,
Hesam Askari,
Esteban Wright,
Paul Sánchez
Abstract:
We carry out experiments of 104 m/s velocity oblique impacts into a granular medium (sand). Impact craters have nearly round rims even at a grazing angle of about $10^\circ$, however, the strength of seismic pulses excited by the impact is dependent upon impact angle, and the ratio between uprange and downrange velocity peaks can be as large as 5, particularly at shallow depths. Crater slope, an o…
▽ More
We carry out experiments of 104 m/s velocity oblique impacts into a granular medium (sand). Impact craters have nearly round rims even at a grazing angle of about $10^\circ$, however, the strength of seismic pulses excited by the impact is dependent upon impact angle, and the ratio between uprange and downrange velocity peaks can be as large as 5, particularly at shallow depths. Crater slope, an offset between crater center and impact site, crater volume, azimuthal variation in ejection angle, seismic pulse shapes and subsurface flow direction are also sensitive to impact angle, but to a much lower degree than subsurface pulse strength. Uprange and downrange pulse peak amplitudes can be estimated from the horizontal and vertical components of the momentum imparted to the medium from the projectile
△ Less
Submitted 22 September, 2023; v1 submitted 3 August, 2023;
originally announced August 2023.
-
Dynamics of co-orbital exoplanets in a first order resonance chain with tidal dissipation
Authors:
Jérémy Couturier,
Philippe Robutel,
Alexandre C. M. Correia
Abstract:
Co-orbital planets (in a $1:1$ mean motion resonance) can be formed within a Laplace resonance chain. Here, we develop a secular model to study the dynamics of the resonance chain $p:p:p+1$, where the co-orbital pair is in a first-order mean motion resonance with the outermost third planet. Our model takes into account tidal dissipation through the use of a Hamiltonian version of the constant time…
▽ More
Co-orbital planets (in a $1:1$ mean motion resonance) can be formed within a Laplace resonance chain. Here, we develop a secular model to study the dynamics of the resonance chain $p:p:p+1$, where the co-orbital pair is in a first-order mean motion resonance with the outermost third planet. Our model takes into account tidal dissipation through the use of a Hamiltonian version of the constant time-lag model, which extends the Hamiltonian formalism of the point-mass case. We show the existence of several families of equilibria, and how these equilibria extend to the complete system. In one family, which we call the main branch, a secular resonance between the libration frequency of the co-orbitals and the precession frequency of the pericentres has unexpected dynamical consequences when tidal dissipation is added. We report the existence of two distinct mechanisms that make co-orbital planets much more stable within the $p:p:p+1$ resonance chain rather than outside it. The first one is due to negative real parts of the eigenvalues of the linearised system with tides, in the region of the secular resonance mentioned above. The second one comes from non-linear contributions of the vector field and it is due to eccentricity damping. These two stabilising mechanisms increase the chances of a still-to-come detection of exoplanets in the co-orbital configuration.
△ Less
Submitted 17 April, 2022;
originally announced April 2022.
-
An analytical model for tidal evolution in co-orbital systems $\,$I. Application to exoplanets
Authors:
Jérémy Couturier,
Philippe Robutel,
Alexandre C. M. Correia
Abstract:
Close-in co-orbital planets (in a 1:1 mean motion resonance) can experience strong tidal interactions with the central star. Here, we develop an analytical model adapted to the study of the tidal evolution of those systems. We use a Hamiltonian version of the constant time-lag tidal model, which extends the Hamiltonian formalism developed for the point-mass case. We show that co-orbital systems un…
▽ More
Close-in co-orbital planets (in a 1:1 mean motion resonance) can experience strong tidal interactions with the central star. Here, we develop an analytical model adapted to the study of the tidal evolution of those systems. We use a Hamiltonian version of the constant time-lag tidal model, which extends the Hamiltonian formalism developed for the point-mass case. We show that co-orbital systems undergoing tidal dissipation either favour the Lagrange or the anti-Lagrange configurations, depending on the system parameters. However, for all range of parameters and initial conditions, both configurations become unstable, although the timescale for the destruction of the system can be larger than the lifetime of the star. We provide an easy-to-use criterion to determine if an already known close-in exoplanet may have an undetected co-orbital companion.
△ Less
Submitted 18 June, 2021;
originally announced June 2021.
-
Squares of Low Maximum Degree
Authors:
Manfred Cochefert,
Jean-François Couturier,
Petr A. Golovach,
Dieter Kratsch,
Daniël Paulusma,
Anthony Stewart
Abstract:
A graph H is a square root of a graph G if G can be obtained from H by adding an edge between any two vertices in H that are of distance 2. The Square Root problem is that of deciding whether a given graph admits a square root. This problem is only known to be NP-complete for chordal graphs and polynomial-time solvable for non-trivial minor-closed graph classes and a very limited number of other g…
▽ More
A graph H is a square root of a graph G if G can be obtained from H by adding an edge between any two vertices in H that are of distance 2. The Square Root problem is that of deciding whether a given graph admits a square root. This problem is only known to be NP-complete for chordal graphs and polynomial-time solvable for non-trivial minor-closed graph classes and a very limited number of other graph classes. We prove that Square Root is O(n)-time solvable for graphs of maximum degree 5 and O(n^4)-time solvable for graphs of maximum degree at most 6.
△ Less
Submitted 27 August, 2016; v1 submitted 22 August, 2016;
originally announced August 2016.
-
Faster algorithms to enumerate hypergraph transversals
Authors:
Manfred Cochefert,
Jean-Francois Couturier,
Serge Gaspers,
Dieter Kratsch
Abstract:
A transversal of a hypergraph is a set of vertices intersecting each hyperedge. We design and analyze new exponential-time algorithms to enumerate all inclusion-minimal transversals of a hypergraph. For each fixed k>2, our algorithms for hypergraphs of rank k, where the rank is the maximum size of a hyperedge, outperform the previous best. This also implies improved upper bounds on the maximum num…
▽ More
A transversal of a hypergraph is a set of vertices intersecting each hyperedge. We design and analyze new exponential-time algorithms to enumerate all inclusion-minimal transversals of a hypergraph. For each fixed k>2, our algorithms for hypergraphs of rank k, where the rank is the maximum size of a hyperedge, outperform the previous best. This also implies improved upper bounds on the maximum number of minimal transversals in n-vertex hypergraphs of rank k>2. Our main algorithm is a branching algorithm whose running time is analyzed with Measure and Conquer. It enumerates all minimal transversals of hypergraphs of rank 3 on n vertices in time O(1.6755^n). Our algorithm for hypergraphs of rank 4 is based on iterative compression. Our enumeration algorithms improve upon the best known algorithms for counting minimum transversals in hypergraphs of rank k for k>2 and for computing a minimum transversal in hypergraphs of rank k for k>5.
△ Less
Submitted 17 October, 2015;
originally announced October 2015.
-
Parameterized Algorithms for Finding Square Roots
Authors:
Manfred Cochefert,
Jean-François Couturier,
Petr A. Golovach,
Dieter Kratsch,
Daniël Paulusma
Abstract:
We show that the following two problems are fixed-parameter tractable with parameter k: testing whether a connected n-vertex graph with m edges has a square root with at most n-1+k edges and testing whether such a graph has a square root with at least m-k edges. Our first result implies that squares of graphs obtained from trees by adding at most k edges can be recognized in polynomial time for ev…
▽ More
We show that the following two problems are fixed-parameter tractable with parameter k: testing whether a connected n-vertex graph with m edges has a square root with at most n-1+k edges and testing whether such a graph has a square root with at least m-k edges. Our first result implies that squares of graphs obtained from trees by adding at most k edges can be recognized in polynomial time for every fixed k>=0; previously this result was known only for k=0. Our second result is equivalent to stating that deciding whether a graph can be modified into a square root of itself by at most k edge deletions is fixed-parameter tractable with parameter k.
△ Less
Submitted 21 October, 2013;
originally announced October 2013.
-
Identification of a Proliferation Gene Cluster Associated with HPV E6/E7 Expression Level and Viral DNA Load in Invasive Cervical Carcinoma
Authors:
Christophe Rosty,
Michal Sheffer,
Dafna Tsafrir,
Nicolas Stransky,
Ilan Tsafrir,
Martine Peter,
Patricia de Cremoux,
Anne de La Rochefordiere,
Remy Salmon,
Thierry Dorval,
Jean Paul Thiery,
Jerome Couturier,
Francois Radvanyi,
Eytan Domany,
Xavier Sastre-Garau
Abstract:
Specific HPV DNA sequences are associated with more than 90% of invasive carcinomas of the uterine cervix. Viral E6 and E7 oncogenes are key mediators in cell transformation by disrupting TP53 and RB pathways. To investigate molecular mechanisms involved in the progression of invasive cervical carcinoma, we performed a gene expression study on cases selected according to viral and clinical param…
▽ More
Specific HPV DNA sequences are associated with more than 90% of invasive carcinomas of the uterine cervix. Viral E6 and E7 oncogenes are key mediators in cell transformation by disrupting TP53 and RB pathways. To investigate molecular mechanisms involved in the progression of invasive cervical carcinoma, we performed a gene expression study on cases selected according to viral and clinical parameters. Using Coupled Two-Way Clustering and Sorting Points Into Neighbourhoods methods, we identified a Cervical Cancer Proliferation Cluster composed of 163 highly correlated transcripts, many of which corresponded to E2F pathway genes controlling cell proliferation, whereas no primary TP53 targets were present in this cluster. The average expression level of the genes of this cluster was higher in tumours with an early relapse than in tumours with a favourable course (P=0.026). Moreover, we found that E6/E7 mRNA expression level was positively correlated with the expression level of the cluster genes and with viral DNA load. These findings suggest that HPV E6/E7 expression level plays a key role in the progression of invasive carcinoma of the uterine cervix via the deregulation of cellular genes controlling tumour cell proliferation. HPV expression level may thus correspond to a biological marker useful for prognosis assessment and specific therapy of the disease.
△ Less
Submitted 15 November, 2005;
originally announced November 2005.