Average Distance Research Articles

Application of Unsupervised Multi-Omic Factor Analysis to Uncover Patterns of Variation and Molecular Processes Linked to Cardiovascular Disease.

Disease mechanisms are usually complex and governed by the interaction of several distinct molecular processes. Complex, multidimensional datasets are a valuable resource to generate more insights into those processes, but the analysis of such datasets can be challenging due to the high dimensionality resulting, for example, from different disease conditions, timepoints, and omics capturing the process at different resolutions. Here, we showcase an approach to analyze and explore such a complex multiomics dataset in an unsupervised way by applying multi-omics factor analysis (MOFA) to a dataset generated from blood samples that capture the immune response in acute and chronic coronary syndromes. The dataset consists of several assays at differing resolutions, including sample-level cytokine data, plasma-proteomics and neutrophil prime-seq, and single-cell RNA-seq (scRNA-seq)data. Further complexity is added by having several different time points measured per patient and several patient subgroups. The analysis workflow outlines how to integrate and analyze the data in several steps: (1) Data pre-processing and harmonization, (2) Estimation of the MOFA model, (3) Downstream analysis. Step 1 outlines how to process the features of the different data types, filter out low-quality features, and normalize them to harmonize their distributions for further analysis. Step 2 shows how to apply the MOFA model and explore the major sources of variance within the dataset across all omics and features. Step 3 presents several strategies for the downstream analysis of the captured patterns, linking them to the disease conditions and potential molecular processes governing those conditions. Overall, we present a workflow for unsupervised data exploration of complex multi-omics datasets to enable the identification of major axes of variation composed of differing molecular features that can also be applied to other contexts and multi-omics datasets (including other assays as presented in the exemplary use case).

Why are inner planets not inclined?

AbstractPoincaré’s work more than one century ago, or Laskar’s numerical simulations from the 1990’s on, have irrevocably impaired the long-held belief that the Solar System should be stable. But mathematical mechanisms explaining this instability have remained mysterious. In 1968, Arnold conjectured the existence of “Arnold diffusion” in celestial mechanics. We prove Arnold’s conjecture in the planetary spatial 4-body problem as well as in the corresponding hierarchical problem (where the bodies are increasingly separated), and show that this diffusion leads, on a long time interval, to some large-scale instability. Along the diffusive orbits, the mutual inclination of the two inner planets is close to $\pi /2$ π / 2 , which hints at why even marginal stability in planetary systems may exist only when inner planets are not inclined. More precisely, consider the normalised angular momentum of the second planet, obtained by rescaling the angular momentum by the square root of its semimajor axis and by an adequate mass factor (its direction and norm give the plane of revolution and the eccentricity of the second planet). It is a vector of the unit 3-ball. We show that any finite sequence in this ball may be realised, up to an arbitrary precision, as a sequence of values of the normalised angular momentum in the 4-body problem. For example, the second planet may flip from prograde nearly horizontal revolutions to retrograde ones. As a consequence of the proof, the non-recurrent set of any finite-order secular normal form accumulates on circular motions – a weak form of a celebrated conjecture of Herman.

The interaction between vegetation patchiness and tidal flows in a shortleaf seagrass meadow

AbstractSeagrasses are critical coastal habitats that provide numerous ecosystem services. The inherent patchiness within meadows exerts a significant influence on the flow‐vegetation interaction, which, in turn, affects the array of services provided by these environments. We present field observations of vegetation distribution and hydrodynamic measurements within and surrounding an approximately 2 m diameter gap (bare sediment) in a fragmented seagrass meadow. We show that variability in mean flows and turbulence is correlated with meadow structure at small (O (100 m)) and large (O (102 m)) spatial scales. Our observations reveal that bare gaps within seagrass meadows lead to faster flows and lower bed elevations. Despite slower flow speeds above the dense seagrass adjacent to the gap, the rates of dissipation of turbulent energy above the vegetation are typically around an order of magnitude larger than above the bare bed. Associated with this enhanced dissipation of turbulent energy, we observed a dominance of down‐deceleration events promoting fluid exchange and mixing and driving mass flux into the canopy. Analysis of directional variograms demonstrates that the major continuity axis within NDVI (normalized difference vegetation index) imagery (used as a proxy for vegetation biomass) coincides with the major axis of flow on both gap and meadow scales. Conversely, along the axis of minor flow variance, vegetation remains dense and exhibits greater uniformity. These findings indicate a feedback mechanism between seagrass meadow patchiness and spatial structure through flow modification, which may be beneficial for plant development.

Automated detection and mapping of crystal tilt using thermal diffuse scattering in transmission electron microscopy

Quantitative interpretation of transmission electron microscopy (TEM) data of crystalline specimens often requires the accurate knowledge of the local crystal orientation. A method is presented which exploits momentum-resolved scanning TEM (STEM) data to determine the local mistilt from a major zone axis. It is based on a geometric analysis of Kikuchi bands within a single diffraction pattern, yielding the center of the Laue circle. Whereas the approach is not limited to convergent illumination, it is here developed using unit-cell averaged diffraction patterns corresponding to high-resolution STEM settings. In simulation studies, an accuracy of approximately 0.1 mrad is found. The method is implemented in automated software and applied to crystallographic tilt and in-plane rotation mapping in two experimental cases. In particular, orientation maps of high-Mn steel and an epitaxially grown La0.7Sr0.3MnO3-SrTiO3 interface are presented. The results confirm the estimates of the simulation study and indicate that tilt mapping can be performed consistently over a wide field of view with diameters well above 100 nm at unit cell real space sampling.

Investigation of the incremental benefits of eccentric collisions in kinetic deflection of potentially hazardous asteroids

In asteroid momentum deflection missions, the presence of ejecta leads to a phenomenon where the system’s momentum appears “amplified” after the impact. This paper makes use of this phenomenon and demonstrates through computational simulations that targeting a point off the geometric center of an asteroid can further enhance the collisional benefit after impact. Due to uncertainties in the attitude of the asteroid and the momentum transfer coefficient (β,γ), this study employs a Monte Carlo approach to address these uncertainties. The results indicate that the strategy proposed in this paper can increase the post-collision deflection distance of the asteroid relative to Earth by an average of 81.05%, while also reducing the standard deviation by an order of magnitude, significantly lowering the uncertainty of the deflection mission. Furthermore, the results show that for certain asteroids particularly sensitive to changes in velocity Δv, blindly targeting their geometric center could result in a 48% probability of reducing the minimum distance to Earth. However, the striking strategy developed in this study can avoid this negative outcome. Finally, based on the computational results, a statistical formula is derived to predict the relative gain of the two strategies, concluding that for asteroids with smaller semi-major axes a, and the interception angle α at impact is greater, the benefits of employing the approach discussed in this paper are greater.

Seismic behavior of innovative reinforced concrete composite shear wall with PEC boundary columns

Considering the inadequate integrity between boundary elements and wall connections in shear walls with PEC columns. A new type of PEC shear wall is proposed to enhance the connection between the PEC columns and the wall by introducing T-section connectors to improve the bearing capacity as well as the deformation capacity. In this study, the seismic behavior of composite shear walls was investigated in terms of the failure process, hysteresis, energy consumption, etc. The test results showed that the specimens had good stiffness and energy dissipation capacity. The combination of PEC columns with RC shear walls under the action of the T-section connector was effective. An increase in the width-to-thickness ratio can lead to shear bonding failure. In the PEC-RCSW structure with a shear-to-span ratio of 1.35, the bearing capacity and deformation capacity of the T-section connectors and PEC columns are slightly higher in the major axis direction than in the minor axis direction. Based on the test results, the finite element model was established and further research was conducted on the shear span ratio and T-section connector size. And provided engineering application suggestions for the innovative composite shear wall shear wall. Finally, flexural and shear capacity calculation formulae of the innovative composite shear wall were proposed.

Single-shot polarization detection with a highly scattering system

Polarization detection plays a significant role in optics. However, the current detection methods usually involve mechanically rotating components, multiple measurement steps, complicated optical design, and precise microfabrication process. To address this issue, we propose a single-shot method to detect the polarization state of light based on a highly scattering system, which is constituted by a spatial light modulator and a highly scattering medium. When the incident light beam shaped by a superimposed wavefront is incident on a highly scattering medium, the foci represented the six components at horizontal, vertical, diagonal, antidiagonal, right circularly polarized, and left circularly polarized directions will appear behind the highly scattering medium simultaneously. By measuring the intensities of these six foci, all the Stokes parameters can be extracted. Taking advantage of the measured Stokes parameters, the orientation angle of major axis, the ellipticity, and the handedness of the polarization ellipse of incident light beam can be determined. Various light beams with different polarization states are detected to demonstrate the viability of the method. The experimental results and theoretical values are in a good agreement. Compared to the existing methods, this approach is fast, free of complicated fabrication, and independent of mechanical movement. The proposed method is expected to promote the development of real-time and broadband polarimetry.

A New Bubble Image Model Based on the Recognition of Bubble Flow

AbstractIn this study, a new ellipse‐fitting algorithm is proposed to achieve the reconstruction of bubble shapes in bubbly flow captured by a high‐speed camera in the gas–liquid two‐phase column reactor. Bubble flow patterns and geometric parameters in the experimental images are recognized and identified successfully, represented by means of the topological parameters. Three logical steps are carried out in detail. First, the area threshold and the circularity factors are established to identify the bubbles whether belonging to a single bubble or not. The overlapping bubbles in images can be separated from single bubbles based on a watershed segmentation algorithm. Second, a single bubble image and an overlapping bubble image are combined into one image. After that, statistical analysis for the size distributions and ellipse area bubbles is performed for further analysis and discussion. The advantage of this algorithm is that it can make use of a set of major and minor axes of an ellipse to capture the ellipse parameters more effectively. Simulation results are well agreed with experimental measurements. Moreover, it can be used to detect many ellipse‐like bubbles that are dispersed in high‐speed camera images, indicating that it is a better strategy for the recognition and identification of bubbly turbulent flow accurately.

Flexural behaviour and design of hybrid cold-formed steel built-up I-section beams

The current study investigates the flexural performance and design of hybrid cold-formed steel built-up I-section beams, which are increasingly popular in the construction industry. These beams combine the benefits of both an "I" section and a closed-box section in their cross-section profile. They are constructed from two identical plain channels placed back-to-back with space between them and top and bottom cover plates, fastened together by screws. A nonlinear finite element (FE) model was developed to consider initial geometric imperfections, material and geometric non-linearities. The model was validated against the test results of cold-formed steel (CFS) built-up beams from the companion paper and literature. Once validation was successful, a comprehensive parametric study using finite element modelling was conducted to generate data on hybrid cold-formed steel built-up I-section beams across a broader range of cross-sectional slenderness, hybrid ratio, aspect ratio of the cross section and thickness of the cover plates. The parametric results were used to analyse the influence of key parameters on the moment carrying capacity and buckling modes. The applicability of the current Direct Strength Method in the AISI S100–16 specifications was evaluated based on these results, revealing that it is unconservative. Additionally, modified design procedures and equations were suggested for hybrid CFS built-up I-section beams vulnerable to local buckling in major axis bending, which were then evaluated through reliability analysis.

Signatures of gas flows – II. Connecting the kinematics of the multiphase circumgalactic medium to galaxy rotation

ABSTRACT The multiphase circumgalactic medium (CGM) hosts critical processes that affect galaxy evolution such as accretion and outflows. We searched for evidence of these phenomena by using the EW co-rotation fraction ($f_{\rm EWcorot}$) to study the kinematic connection between the multiphase CGM and host galaxy rotation. We examined CGM absorption from Hubble Space Telescope /Cosmic Origins Spectrograph (including, but not limited to, Si ii, C ii, Si iii, C iii, and O vi) within $21\le D\le ~276$ kpc of 27 galaxies. We find the median $f_{\rm EWcorot}$ for all ions is consistent within errors and the $f_{\rm EWcorot}$ increases with increasing N$({{{\rm H} \rm{\small I}}})$. The $f_{\rm EWcorot}$ of lower ionization gas decreases with increasing $D/R_{\rm vir}$, while O vi and H i are consistent with being flat. The $f_{\rm EWcorot}$ varies minimally as a function of azimuthal angle and is similar for all ions at a fixed azimuthal angle. The larger number of O vi detections enabled us to investigate where the majority of co-rotating gas is found. Highly co-rotating O vi primarily resides along the galaxies’ major axis. Looking at the $f_{\rm EWcorot}$ as a function of ionization potential (${{\rm d}{({f_{\rm EWcorot}})}}/{{\rm d}{(\rm eV)}}$), we find a stronger co-rotation signature for lower ionization gas. There are suggestions of a connection between the CGM metallicity and major axis co-rotation where low-ionization gas with higher $f_{\rm EWcorot}$ exhibits lower metallicity and may trace large-scale filamentary inflows. Higher ionization gas with higher $f_{\rm EWcorot}$ exhibits higher metallicity and may instead trace co-planar recycled gas accretion. Our results stress the importance of comparing absorption originating from a range of ionization phases to differentiate between various gas flow scenarios.

Mechanical behavior of dumbbell-type metamaterials constructed with elliptic curves

In this paper, based on the traditional concave ortho-hexagonal honeycomb (CHH), its curvilinear modified design was introduced into the elliptic curve configuration, and a novel dumbbell-shaped honeycomb (DH) was proposed. All specimens were prepared by Liquid Crystal Display (LCD) light-cured 3D printing. The aim of this study is to investigate the deformation and damage behavior and mechanical properties of DH and CHH under in-plane quasi-static uniaxial compression and tension. The effects of variation of the cell parameters semi-major axis length a and semi-minor axis length b on the deformation, damage behavior and mechanical properties of DH have been investigated. The experimental results show that the DH incorporating the elliptical curvilinear design has a reduced internal compressive strength under compression, but the negative Poisson’s ratio property is enhanced, which is 17% higher compared to the CHH, and the tensile strength under tension is reduced compared to the CHH, but the tensile ultimate strain is increased by 25% compared to the CHH. This will help to improve the problem of stress concentration and negative Poisson’s ratio in the process of deformation of conventional structures. Differences in the a and b parameters in the DH structure will have different effects on its deformation-fracture behavior and mechanical properties.

Ellipsoidal micelle formation of quaternary ammonium salt-based gemini surfactants: structural analysis through small-angle X-ray scattering.

Gemini surfactants exhibit better adsorption and aggregation properties than those of monomeric surfactants. However, to enhance the functional properties of gemini surfactants, the effect of spacer structures on their aggregation behavior must be elucidated. Small-angle X-ray scattering (SAXS) has been performed to study the aggregate structures of monomeric surfactants, but its application has not been expanded to the analysis of gemini surfactants. Therefore, in this study, we performed a structural analysis of aggregates formed by quaternary ammonium salt-based gemini surfactants with different spacer structures in aqueous solutions through viscosity, dynamic light scattering (DLS), and SAXS measurements. We also investigated the effects of the spacer structure and surfactant concentration on the aggregation behavior. The DLS results indicated that gemini surfactants with spacers containing cyclic structures, such as diethylene and triethylene chains, formed small micelles (several nanometers in size) at the limiting concentration for dissolution in water. In contrast, gemini surfactants with nitrogen and oxygen atoms at the center of the spacer formed ellipsoidal micelles at low concentrations, as shown by SAXS results. The core and overall radii of the minor and major axes of the ellipsoidal micelles decreased with increasing surfactant concentration and were larger for spacers with two ethylene chains connected to central nitrogen and oxygen atoms than for spacers with rigid diethylene and ethylene chains connected to a central nitrogen atom.

SPGAN Optimized by Piranha Foraging Optimization for Thyroid Nodule Classification in Ultrasound Images.

In this research work, Semantic-Preserved Generative Adversarial Network optimized by Piranha Foraging Optimization for Thyroid Nodule Classification in Ultrasound Images (SPGAN-PFO-TNC-UI) is proposed. Initially, ultrasound images are gathered from the DDTI dataset. Then the input image is sent to the pre-processing step. During pre-processing stage, the Multi-Window Savitzky-Golay Filter (MWSGF) is employed to reduce the noise and improve the quality of the ultrasound (US) images. The pre-processed output is supplied to the Generalized Intuitionistic Fuzzy C-Means Clustering (GIFCMC). Here, the ultrasound image's Region of Interest (ROI) is segmented. The segmentation output is supplied to the Fully Numerical Laplace Transform (FNLT) to extract the features, such as geometric features like solidity, orientation, roundness, main axis length, minor axis length, bounding box, convex area, and morphological features, like area, perimeter, aspect ratio, and AP ratio. The Semantic-Preserved Generative Adversarial Network (SPGAN) separates the image as benign or malignant nodules. Generally, SPGAN does not express any optimization adaptation methodologies for determining the best parameters to ensure the accurate classification of thyroid nodules. Therefore, the Piranha Foraging Optimization (PFO) algorithm is proposed to improve the SPGAN classifier and accurately identify the thyroid nodules. The metrics, like F-score, accuracy, error rate, precision, sensitivity, specificity, ROC, computing time is examined. The proposed SPGAN-PFO-TNC-UI method attains 30.54%, 21.30%, 27.40%, and 18.92% higher precision and 26.97%, 20.41%, 15.09%, and 18.27% lower error rate compared with existing techniques, like Thyroid detection and classification using DNN with Hybrid Meta-Heuristic and LSTM (TD-DL-HMH-LSTM), Quantum-Inspired convolutional neural networks for optimized thyroid nodule categorization (QCNN-OTNC), Thyroid nodules classification under Follow the Regularized Leader Optimization based Deep Neural Networks (CTN-FRL-DNN), Automatic classification of ultrasound thyroids images using vision transformers and generative adversarial networks (ACUTI-VT-GAN) respectively.

Neoichnology of Rhinella dorbignyi (Bufonidae) burrows: improving the recognition and interpretation of toad burrows

A neoichnological characterization of toad burrows from organic-rich soils is presented in order to establish key ichnological signatures for their recognition in the geological record. Dorbigny’s Toad (Rhinella dorbignyi) produces J-shaped burrows with a single elliptical opening resembling the ichnogenus Macanopsis. The burrow varies from 15 to 20 cm in length and has a constant diameter along the shaft, enlarged in the chamber. The shaft is elliptical with a diameter of 2.5-3.0 cm in the major axis and 1.0–1.5 cm in the minor axis. The chamber is semicircular, with a distinct floor, and an average diameter of 5 cm in the major axis. A thick (3–5 mm), compacted, smooth, and uniform inorganic lining made of soil covers the burrow throughout its length. Surficial features, such as horizontal and vertical millimeter-scale striations and ridges, occur in the burrow surface and are produced by keratinized granules and spikes on the skin of R. dorbignyi. Those surficial features are generated during the toad’s movements inside the burrow as well as anchoring during ambush predation behavior. The main distinct features observed on the burrow surface are well-developed circular pit impressions produced by the keratinized tubercles present in the toad’s ankles during climbing activity. These characteristics encompass useful ichnotaxobases for recognition and distinguishing toad burrows from J-shaped, Macanopsis-like burrows produced by other soil-dwelling animals. The burrow architecture of R. dorbignyi suggests that J-shaped burrows are efficient structures for ambush predatory behavior in terrestrial ecosystems, as observed in spiders and scorpions.

Expression and significance of cystine transporter SLC7A11/xCT in early colorectal cancer specimens from endoscopic submucosal dissection.

The SLC7A11/xCT cystine transporter is intricately linked with ferroptosis. By mediating intracellular cystine flux, it regulates oxidative stress within neoplastic cells, thereby curtailing ferroptosis and influencing the emergence of colorectal cancer. This study aimed to gauge the SLC7A11/xCT expression across various tumorigenic stages in early colorectal adenocarcinoma tissues, shedding light on its specific role in the genesis of these early malignancies. Sixty specimens that underwent endoscopic submucosal dissection (ESD) resection with pathologic diagnosis of colorectal adenocarcinoma were collected. SLC7A11/xCT expression was pinpointed using immunohistochemistry, and correlations with the patients' clinical-pathological features were drawn. Additionally, a comprehensive bioinformatics assessment was undertaken to discern differential SLC7A11/xCT expressions across a spectrum of cancers. Immunohistochemical assessments unveiled a pronounced cytoplasmic SLC7A11/xCT expression, manifesting as a brownish-yellow hue, particularly in nascent colorectal cancer samples. Its expression was discernibly correlated with both patient gender and adenocarcinoma differentiation grade (P<0.05). Nevertheless, factors such as patient age, tumor localization, infiltration depth, diameter, adjacent adenoma histology, its major axis, and dysplasia degree bore no statistical significance with SLC7A11/xCT levels (P>0.05). Bioinformatics insights pointed to an upregulated SLC7A11/xCT expression across diverse malignancies, inclusive of colon adenocarcinoma, esophageal cancer, acute myeloid leukemia, lung squamous cell carcinoma, colorectal cancer, and endometrial cancer (P<0.05). Elevated SLC7A11/xCT expression marks early colorectal adenocarcinoma, with the intensity of this expression being intertwined with the patient's gender and the tumor's differentiation grade. It is postulated that colorectal cancer cells might amplify SLC7A11/xCT to stymie ferroptosis, thus fostering neoplastic proliferation, metastasis, and cellular stemness.

Multi-factor coupled forest fire model based on cellular automata

The risk of forest fires is substantial due to uneven precipitation distributions and abnormal climate change. This study employs cellular automata principles to analyze forest fire behavior, taking into account meteorological elements, combustible material types, and terrain slopes. The Wang Zhengfei model is utilized to compute fire spread speed, and a multifactor coupled forest fire model is developed. Comparisons with experimental data show a mean calculated fire spread speed of 0.69 m/min, which is consistent with the experimental results. Using the forest fire in Anning city, Yunnan Province, as a case study with a mean burned area of 2281 ha, the burned area, rate of change in burned area, and burning area demonstrated an increasing trend, with fluctuating states in the rate of change of the burning area. Employing the controlled variable method to examine forest fire spreading patterns under varying factors such as wind speed, vegetation type, and maximum slope reveals that under wind influence, the fire site adopts an elliptical shape with the downwind direction as the major axis. Quantitatively, when the wind speed increases from 2 m/s to 10 m/s, the burned area expands by a factor of 1.37. The ratio of the combustible material configuration coefficient to the burned area remains consistent across the different vegetation types, and the burned area increases by a factor of 1.92 when the maximum slope increases from 5° to 25°.

Quantitative classification evaluation model for tight sandstone reservoirs based on machine learning

Tight sandstone reservoirs are a primary focus of research on the geological exploration of petroleum. However, many reservoir classification criteria are of limited applicability due to the inherent strong heterogeneity and complex micropore structure of tight sandstone reservoirs. This investigation focused on the Chang 8 tight reservoir situated in the Jiyuan region of the Ordos Basin. High-pressure mercury intrusion experiments, casting thin sections, and scanning electron microscopy experiments were conducted. Image recognition technology was used to extract the pore shape parameters of each sample. Based on the above, through grey relational analysis (GRA), analytic hierarchy process (AHP), entropy weight method (EWM) and comprehensive weight method, the relationship index Q1 between initial productivity and high pressure mercury injection parameters and the relationship index Q2 between initial productivity and pore shape parameters are obtained by fitting. Then a dual-coupled comprehensive quantitative classification prediction model for tight sandstone reservoirs was developed based on pore structure and shape parameters. A quantitative classification study was conducted on the target reservoir, analyzing the correlation between reservoir quality and pore structure and shape parameters, leading to the proposal of favourable exploration areas. The research results showed that when Q1 ≥ 0.5 and Q2 ≥ 0.5, the reservoir was classified as type I. When Q1 > 0.7 and Q2 > 0.57, it was classified as type I1, indicating a high-yield reservoir. When 0.32 < Q1 < 0.47 and 0.44 < Q2 < 0.56, was classified as type II. When 0.1 < Q1 < 0.32 and 0.3 < Q2 < 0.44, it was classified as type III. Type I reservoirs exhibit a zigzag pattern in the northwest part of the study area. Thus, the northwest should be prioritized in actual exploration and development. Additionally, the initial productivity of tight sandstone reservoirs showed a positive correlation with the porosity, permeability, sorting coefficient, coefficient of variation, and median radius. Conversely, it demonstrated a negative correlation with the median pressure and displacement pressure. The perimeters of pores, their circularity, and the length of the major axis showed a positive correlation with the porosity, permeability, sorting coefficient, coefficient of variation, and median radius. On the other hand, they exhibited a negative correlation with the median pressure and displacement pressure. This study quantitatively constructed a new classification and evaluation system for tight sandstone reservoirs from the perspective of microscopic pore structure, achieving an overall model accuracy of 93.3%. This model effectively predicts and evaluates tight sandstone reservoirs. It provides new guidance for identifying favorable areas in the study region and other tight sandstone reservoirs.

Influence of planets on debris discs in star clusters – II. The impact of stellar density

ABSTRACT We present numerical simulations of planetary systems in star clusters with different initial stellar densities, to investigate the impact of the density on debris disc dynamics. We use lps+ to combine N-body codes nbody6++gpu and rebound for simulations. We simulate debris discs with and without a Jupiter-mass planet at 50 au, in star clusters with $N=$ 1k–64k stars. The spatial range of the remaining planetary systems decreases with increasing N. As cluster density increases, the planet’s influence range first increases and then decreases. For debris particles escaping from planetary systems, the probability of their direct ejection from the star cluster decreases as their initial semimajor axis ($a_0$) or the cluster density increases. The eccentricity and inclination of surviving particles increase as cluster density increases. The presence of a planet leads to lower eccentricities and inclinations of surviving particles. The radial density distribution of the remaining discs decays exponentially in sparse clusters. We derive a general expression of the gravitational encounter rate. Our results are unable to directly explain the scarcity of debris discs in star clusters. Nevertheless, given that many planetary systems have multiple planets, the mechanism of the planet-cluster combined gravitational influence on the disc remains appealing as a potential explanation.

Flexible Bridge-Conic Strategies for Moon-to-Moon Analytical Transfer Design

The multi-burn transfer framework based on the moon-to-moon analytical transfer method allows rapid patching of low-energy trajectories in the vicinity of different moons in a planetary system. However, the semimajor axis and eccentricity of intermediate arcs of one transfer mission are fixed in the existent multi-burn strategies. Additional propellant consumption is required to enable the transfers at all departure epochs, and a nonanalytical formulation still remains for transferring to spatial periodic orbits. To address these difficulties, new multi-burn strategies with variable intermediate arcs are proposed to meet the transfer conditions for different types of arrival orbits. Specifically, one or two intermediate arcs with one variable parameter are introduced and can be determined analytically. In numerical simulations, the proposed flexible bridge-conic strategies are applied to design transfers between libration point orbits in both the Jovian and Saturnian systems. Final results validate the effectiveness of the proposed strategies and indicate the important advantage in propellant saving compared with the existent multi-burn strategies.

Ultraconfined oblate hard particles between hybrid penetrable walls.

We have investigated, by Monte Carlo simulation, the orientational structure of very thin films of a discotic liquid crystal (DLC) confined between hybrid walls of controllable penetrability, as a function of wall separation L_{z}. Our purpose was to clarify whether, as predicted by continuum theory, the preferred orientation of the DLC is uniform, changes linearly, or changes discontinuously, when L_{z} and the anchoring strengths at either wall are changed. The model consists of oblate hard Gaussian overlap (HGO) particles: each wall sees a particle as a disk of zero thickness and diameter D less than or equal to that of the actual particle σ_{0}, embedded inside the particle and located halfway along, and perpendicular to, its minor axis. This provides a particle-level mechanism to control the anchoring properties of the walls, from planar (edge-on) for D∼0 to homeotropic (face-on) for D∼σ_{0}, which can be done independently at either wall. As in our earlier work [C. Anquetil-Deck et al., J. Phys. Chem. B 124, 7709 (2020)1520-610610.1021/acs.jpcb.0c05027], which was restricted to L_{z}=6σ_{0}, depending on the values of D_{s}≡D/σ_{0} at the top (D_{s}^{t}) and bottom (D_{s}^{b}) walls, we find domains in (D_{s}^{b},D_{s}^{t}) space in which particle alignment is uniform planar (UP), uniform homeotropic (UH), or varies linearly from planar at one wall to homeotropic at the other (L), but no bistable or tristable regions are identified between these domains. Most importantly, there appears never to occur an abrupt change of the LC orientation when the walls strongly favor different anchorings, in general agreement with the scenario proposed by Velasco and co-workers [D. de las Heras et al., Phys. Rev. E 79, 011712 (2009)1539-375510.1103/PhysRevE.79.011712], but in contrast to the behavior of equivalent calamitic systems [F. Barmes et al., Phys. Rev. E 69, 061705 (2004)1539-375510.1103/PhysRevE.69.061705; Phys. Rev. E, 71, 021705 (2005)1539-375510.1103/PhysRevE.71.021705; C. Anquetil-Deck et al., Phys. Rev. E 86, 041707 (2012)1539-375510.1103/PhysRevE.86.041707]. However, for the thinnest films investigated (L_{z}=2σ_{0}), the system is unable to accommodate a rotation of the preferred particle orientation from one wall to the other and adopts instead a tilted configuration, similar to that reported earlier for Gay-Berne films in symmetric confinement [T. Gruhn et al., Thin Solid Films 330, 46 (1998)0040-609010.1016/S0040-6090(98)00799-8; Mol. Phys. 93, 681 (1998)10.1080/002689798169014] but which, as far as we know, has been missed in most earlier work.