Search

The 21 cm transition from neutral Hydrogen promises to be the best observational probe of the Epoch of Reionisation (EoR). This has led to the construction of low-frequency radio interferometric arrays, such as the Hydrogen Epoch of Reionization Array (HERA), aimed at systematically mapping this emission for the first time. Precision calibration, however, is a requirement in 21 cm radio observations. Due to the spatial compactness of HERA, the array is prone to the effects of mutual coupling, which inevitably lead to non-smooth calibration errors that contaminate the data. When unsmooth gains are used in calibration, intrinsically spectrally-smooth foreground emission begins to contaminate the data in a way that can prohibit a clean detection of the cosmological EoR signal. In this paper, we show that the effects of mutual coupling on calibration quality can be reduced by applying custom time-domain filters to the data prior to calibration. We find that more robust calibration solutions are derived when filtering in this way, which reduces the observed foreground power leakage. Specifically, we find a reduction of foreground power leakage by 2 orders of magnitude at k=0.5.

We implemented symbolic regression techniques to identify suitable analytical functions that map various properties of neutron stars (NSs), obtained by solving the Tolman-Oppenheimer-Volkoff (TOV) equations, to a few key parameters of the equation of state (EoS). These symbolic regression models (SRMs) are then employed to perform Bayesian inference with a comprehensive dataset from nuclear physics experiments and astrophysical observations. The posterior distributions of EoS parameters obtained from Bayesian inference using SRMs closely match those obtained directly from the solutions of TOV equations. Our SRM-based approach is approximately 100 times faster, enabling efficient Bayesian analyses across different combinations of data to explore their sensitivity to various EoS parameters within a reasonably short time.

Interferometric experiments designed to detect the highly redshifted 21-cm signal from neutral hydrogen are producing increasingly stringent constraints on the 21-cm power spectrum, but some k-modes remain systematics-dominated. Mutual coupling is a major systematic that must be overcome in order to detect the 21-cm signal, and simulations that reproduce effects seen in the data can guide strategies for mitigating mutual coupling. In this paper, we analyse 12 nights of data from the Hydrogen Epoch of Reionization Array and compare the data against simulations that include a computationally efficient and physically motivated semi-analytic treatment of mutual coupling. We find that simulated coupling features qualitatively agree with coupling features in the data; however, coupling features in the data are brighter than the simulated features, indicating the presence of additional coupling mechanisms not captured by our model. We explore the use of fringe-rate filters as mutual coupling mitigation tools and use our simulations to investigate the effects of mutual coupling on a simulated cosmological 21-cm power spectrum in a "worst case" scenario where the foregrounds are particularly bright. We find that mutual coupling contaminates a large portion of the "EoR Window", and the contamination is several orders-of-magnitude larger than our simulated cosmic signal across a wide range of cosmological Fourier modes. While our fiducial fringe-rate filtering strategy reduces mutual coupling by roughly a factor of 100 in power, a non-negligible amount of coupling cannot be excised with fringe-rate filters, so more sophisticated mitigation strategies are required., Comment: 19 pages, 12 figures, submitted to MNRAS

A Bayesian method is used in this extensive work to generate a large set of minimally constrained equations of state (EOSs) for matters in neutron stars (NS). These EOSs are analyzed for their correlations with key NS properties, such as the tidal deformability, radius, and maximum mass, within the mass range of $1.2-2M_\odot$. The observed connections between the pressure of $\beta$-equilibrated matter and the properties of neutron stars at different densities offer significant insights into the behavior of NS matter in a nearly model-independent manner. The study also examines the influence of various factors on the correlation of symmetry energy parameters, such as slope and curvature parameters at saturation density ($\rho_0=0.16 ~\text{fm}^{-3}$) with the tidal deformability and radius of neutron stars. This study investigates the robustness of the observed correlations by considering the distributions and interdependence of symmetry energy parameters. Furthermore, the utilization of Principal Component Analysis (PCA) is employed to unveil the complicated relationship between various nuclear matter parameters and properties of neutron stars. This analysis highlights the importance of employing multivariate analysis techniques in order to comprehend the variety in tidal deformability and radius observed across distinct masses of NS. This comprehensive study aims to establish a connection between the parameters of nuclear matter and the properties of neutron stars, providing significant insights into the behavior of NS matter across different circumstances., Comment: PhD Thesis [Examiner:- Prof. Mariana Dutra ]

We have attempted to mitigate the challenge of connecting the neutron star (NS) properties with the nuclear matter parameters that describe equations of state (EoSs). The efforts to correlate various neutron star properties with individual nuclear matter parameters have been inconclusive. A Principal Component Analysis is employed as a tool to uncover the connection between multiple nuclear matter parameters and the tidal deformability as well as the radius of neutron stars within the mass range of $1.2-1.8M_\odot$. The essential EOSs for neutron star matter at low densities have been derived using both uncorrelated uniform distributions and minimally constrained joint posterior distributions of nuclear matter parameters. For higher densities ($\rho > 0.32$fm$^{-3}$), the EOSs have been established through a suitable parameterization of the speed of sound, which consistently maintains causality and gradually approaches the conformal limit. Our analysis reveals that in order to account for over 90\% of the variability in NS properties, it is crucial to consider two or more principal components, emphasizing the significance of employing multivariate analysis. To explain the variability in tidal deformability needs a greater number of principal components compared to those for the radius at a given NS mass. The contributions from iso-vector nuclear matter parameters to the tidal deformability and radius of NS decrease by $\sim$ 25\% with the increase in mass of NS from 1.2$M_\odot$ to 1.8$M_\odot$., Comment: 12 Pages, 7 Figures and 2 Tables (Accepted in Physical Review D)

The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets., Comment: Under review by MNRAS; 28 pages, 16 figures

The impacts of various symmetry energy parameters on the properties of neutron stars (NSs) have been recently investigated, and the outcomes are at variance, as summarized in Table III of Phys. Rev. D 106, 063005 (2022). We have systematically analyzed the correlations of slope and curvature parameters of symmetry energy at the saturation density ($\rho_0=0.16 \text{fm}^{-3}$) with the tidal deformability and stellar radius of non-spinning neutron stars in the mass range of $1.2 - 1.6 M_\odot$ using a large set of minimally constrained equations of state (EoSs). The EoSs at low densities correspond to the nucleonic matter and are constrained by empirical ranges of a few low-order nuclear matter parameters from the finite nuclei data and the pure neutron matter EoS from chiral effective field theory. The EoSs at high densities ($\rho > 1.5 - 2\rho_0$) are obtained by a parametric form for the speed of sound that satisfies the causality condition. Several factors affecting the correlations between the NS properties and the individual symmetry energy parameters usually encountered in the literature are considered. These correlations are quite sensitive to the choice of the distributions of symmetry energy parameters and their interdependence. But, variations of NS properties with the pressure of $\beta -$ equilibrated matter at twice the saturation density remain quite robust which maybe due to the fact that the pressure depends on the combination of multiple nuclear matter parameters that describe the symmetric nuclear matter as well as the density dependence of the symmetry energy. Our results are practically insensitive to the behavior of EoS at high densities., Comment: 14 pages, 8 figures, 4 tables (Published in PRC)

With an effective chiral model, we investigate how the mesonic cross couplings $\sigma-\rho$ and $\omega-\rho$ affect the density content of the symmetry energy and its higher-order slope parameters. Earlier mentioned cross-couplings are crucial to controlling the density content of symmetry energy. For this purpose, we did a case study for different values of the symmetry energy $J_1$, defined at density 0.1 fm$^{-3}$ in the range (23.4 - 25.2) for a fixed value of the slope of the symmetry energy $L_0 = 60$ MeV at saturation density and investigate its effect on the higher-order coefficients and their influence on the underlying equation of state. We found that the model with $J_1= 24.6$ MeV is more favorable with the pure neutron matter (PNM) constraints obtained from $\chi$EFT calculations. In addition, we show that all of our models predict a monotonically increasing speed of sound up to four times the saturation density. The speed of sound decreases/saturates above that point and approaches the conformal limit approximately $\sqrt{1/3}~c$ at the center of the maximum mass star., Comment: The manuscript contain 9 figures and 2 tables

Motivated by the recent observation of high-mass pulsars ($M \simeq 2 M_{\odot}$), we employ the $\sigma$-cut potential on the equation of state (EOS) of high-density matter and the properties of neutron stars within the relativistic mean-field (RMF) model using TM1$^{*}$ parameter set. The $\sigma$-cut potential is known to reduce the contributions of the $\sigma$ field, resulting in a stiffer EOS at high densities and hence leading to larger neutron star masses without affecting the properties of nuclear matter at normal saturation density. We also analyzed the effect of the same on pure neutron matter and also on the neutron star matter with and without hyperonic core and compared it with the available theoretical, experimental, and observational data. The corresponding tidal deformability ($\Lambda_{1.4}$) is also calculated. With the choice of meson-hyperon coupling fixed to hypernuclear potentials, we obtain $\approx 10~\%$ increase in mass by employing the $\sigma$-cut potential for $f_{s} = 0.6$. Our results are in good agreement with various experimental constraints and observational data, particularly with the GW170817 data., Comment: 7 Pages, 6 Figures and 1 Table (Accepted in Phys. Rev. C)

This paper presents the High-Performance computing efforts with FPGA for the accelerated pulsar/transient search for the SKA. Case studies are presented from within SKA and pathfinder telescopes highlighting future opportunities. It reviews the scenario that has shifted from offline processing of the radio telescope data to digitizing several hundreds/thousands of antenna outputs over huge bandwidths, forming several 100s of beams, and processing the data in the SKA real-time pulsar search pipelines. A brief account of the different architectures of the accelerators, primarily the new generation Field Programmable Gate Array-based accelerators, showing their critical roles to achieve high-performance computing and in handling the enormous data volume problems of the SKA is presented here. It also presents the power-performance efficiency of this emerging technology and presents potential future scenarios., Comment: Accepted for publication in the JoAA special issue on SKA (2023)

We apply Bayesian approach to construct a large number of minimally constrained equations of state (EOSs) and study their correlations with a few selected properties of a neutron star (NS). Our set of minimal constraints includes a few basic properties of saturated nuclear matter and low-density pure neutron matter EOS which is obtained from a precise next-to-next-to-next-to-leading-order (N$^{3}$LO) calculation in chiral effective field theory. The tidal deformability and radius of NS with mass $1-2 M_\odot$ are found to be strongly correlated with the pressure of $\beta$-equilibrated matter at densities higher than the saturation density ($\rho_0 = 0.16$ fm$^{-3}$) in a nearly model-independent manner. These correlations are employed to parametrize the pressure for $\beta$-equilibrated matter, around 2$\rho_0$, as a function of neutron star mass and the corresponding tidal deformability. The maximum mass of neutron star is also found to be strongly correlated with the pressure of $\beta$-equilibrated matter at densities $\sim 4.5\rho_0$., Comment: 19 pages, 7 figures, Accepted in Physical Review D. arXiv admin note: text overlap with arXiv:2110.15776

The nuclear matter parameters (NMPs), those underlie in the construction of the equation of state (EoS) of neutron star matter, are not directly accessible. The Bayesian approach is applied to reconstruct the posterior distributions of NMPs from the EoS of neutron star matter. The constraints on lower-order parameters as imposed by the finite nuclei observables are incorporated through appropriately chosen prior distributions. The calculations are performed with two sets of pseudo data on the EoS whose true models are known. The median values of second or higher order NMPs show sizeable deviations from their true values and associated uncertainties are also larger. The sources of these uncertainties are intrinsic in nature, identified as (i) the correlations among various NMPs and (ii) the variations in the EoS of symmetric nuclear matter, symmetry energy, and the neutron-proton asymmetry in such a way that the neutron star matter EoS remain almost unaffected., Comment: 11 pages, 7 figures, Published in PRC

Aiming to fill a crucial gap in our observational knowledge of the early Universe, experiments around the world continue to attempt to verify the claimed detection of the redshifted 21-cm signal from Cosmic Dawn by the EDGES experiment. This sky-averaged or 'global' signal from neutral hydrogen should be detectable at low radio frequencies (50-200 MHz), but is difficult to measure due to bright foreground emission and difficulties in reaching the required levels of instrumental-calibration precision. In this paper we outline our progress toward using a novel new method to measure the global redshifted 21-cm signal. Motivated by the need to use alternative methods with very different systematic errors to EDGES for an independent result, we employ an array of closely-spaced antennas to measure the global sky signal interferometrically, rather than using the conventional approach with a single antenna. We use simulations to demonstrate our newly-developed methods and show that, for an idealised instrument, a 21-cm signal could theoretically be extracted from the visibilities of an array of closely-spaced dipoles. We verify that our signal-extraction methods work on real data using observations made with a Square Kilometre-Array-like prototype; the Engineering Development Array-2. Finally, we use the lessons learned in both our simulations and observations to lay out a clear plan for future work, which will ultimately lead to a new global redshifted 21-cm instrument: the All-Sky SignAl Short-Spacing INterferometer (ASSASSIN)., Comment: 18 pages, 13 figures

We present new HI data of the dwarf galaxy KK 69, obtained with the Giant Metrewave Radio Telescope (GMRT) with a signal-to-noise ratio that almost double previous observations. We carried out a Gaussian spectral decomposition and stacking methods to identify the cold neutral medium (CNM) and the warm neutral medium (WNM) of the HI gas. We found that 30% of the total HI gas, which corresponds to a mass of approx. 10^7 Mo, is in the CNM phase. The distribution of the HI in KK 69 is not symmetric. Our GMRT HI intensity map of KK 69 overlaid onto a Hubble Space Telescope image reveals an offset of approx. 4 kpc between the HI high-density region and the stellar body, indicating it may be a dwarf transitional galaxy. The offset, along with the potential truncation of the HI body, are evidence of interaction with the central group spiral galaxy NGC 2683, indicating the HI gas is being stripped from KK 69. Additionally, we detected extended HI emission of a dwarf galaxy member of the group as well as a possible new galaxy located near the north-eastern part of the NGC 2683 HI disk., Comment: 13 pages, 7 figures, accepted for publication in: Astrophysics and Space Science

Ba(Fe$_{1/2}$Nb$_{1/2}$)O$_3$ (BFN) ceramics are considered to be promising for technological applications owing to their high dielectric constant over a wide range of temperatures. However, there exists considerable discrepancy over the structural details. We address this discrepancy through a combined x-ray diffraction at room temperature and neutron powder diffraction measurements in the range from 5K up to room temperature, supplemented by a comparative analysis of the earlier reported structures. Our study reveals a cubic structure with space group Pm$\bar{3}$m at all measured temperatures. Further, the x-ray near edge structure and the extended x-ray absorption fine structure studies on the local environment of the Fe ions is consistent with the cubic symmetry. An appropriate value of $U$ for DFT+$U$ calculations is obtained by comparison with x-ray absorption spectroscopy, which agrees well with the earlier reported electronic properties., Comment: 15 pages, 9 figures

With 30 antennas and a maximum baseline length of 25 km, the Giant Metrewave Radio Telescope (GMRT) is the premier low-frequency radio interferometer today. We have carried out a study of possible expansions of the GMRT, via adding new antennas and installing focal plane arrays (FPAs), to improve its point-source sensitivity, surface brightness sensitivity, angular resolution, field of view, and U-V coverage. We have carried out array configuration studies, aimed at minimizing the number of new GMRT antennas required to obtain a well-behaved synthesized beam over a wide range of angular resolutions for full-synthesis observations. This was done via two approaches, tomographic projection and random sampling, to identify the optimal locations for the new antennas. We report results for the optimal locations of the antennas of an expanded array (the "EGMRT"), consisting of the existing 30 GMRT antennas, 30 new antennas at short distances, $\leq 2.5$ km from the array centre, and 26 new antennas at long distances, $\approx 5-25$ km from the array centre. The collecting area and the field of view of the proposed EGMRT array would be larger by factors of, respectively, $\approx 3$ and $\approx 30$, than those of the GMRT. Indeed, the EGMRT continuum sensitivity and survey speed with 550-850 MHz FPAs installed on the 45 antennas within a distance of $\approx 2.5$ km of the array centre would be far better than those of any existing interferometer, and comparable to the sensitivity and survey speed of Phase-1 of the Square Kilometre Array., Comment: 16 pages, 13 figures; accepted for publication in MNRAS

Rare earth Manganite Nickelite double perovskites are prepared by solgel assisted combustion route, where, R varies from La, Pr, Nd, Sm, Gd, Tb, Dy, Y, and Ho. The samples have been systematically investigated using powder xray diffraction, Raman spectroscopy, ultraviolet visible spectroscopy, magnetization, and synchrotron based xray absorption spectroscopy measurements. All compounds in the family crystallize in the monoclinic structure and the monoclinic distortion enhances with decreasing trivalent rare earth radii. The magnetic ordering temperature, Tc, decreases from 270 K for La to 80 K for Ho samples as the rare earth radii decrease from 0.116nm for La to 0.102nm in case of Ho. An additional anomaly is observed in samples containing Nd, Sm, Tb, and Dy at lower temperatures, which originates from the 3d to 4f coupling between Mn and Ni and Nd, Sm, Tb and Dy magnetic moments. Further, high saturation magnetization is achieved for all samples indicating that they are atomically ordered and have less anti site disorders. Upon decreasing the size of rare earth ions, the local structure shows an expansion of NiO6 octahedra and almost unchanged MnO6 octahedra. Xray absorption near edge spectroscopy reveals that majority of Ni and Mn ions are in positive 2 and positive 4 valence states in all the samples. Raman spectra of RNMO show a softening of phonon modes resulting in the elongation of Ni to O and Mn to O bond lengths. Finally, a correlation between lattice parameters, structural distortion, octahedral tilting, superexchange angle, and electronic band gap, Curie temperature, and the rare earth ionic radius is established., Comment: 29

Local and long range structure, optical and photoluminescence properties of sol-gel synthesized Ce1-xNixO2 nanostructures have been studied. The crystal structure, lattice strain and crystallite size have been analyzed. A decrease in lattice parameter may be attributed to substitution of Ce with smaller Ni ion. UV-Vis measurement is used for studying the effect of Ni substitution on bandgap and disorder. The bandgap decreases with Ni substitution and disorder increases. The PL spectra show five major peaks attributed to various defect states. The PL emission decreases with Ni substitution owing to increase in defects which acts as emission quenching centers. The lattice disorder and defects have been studied using Raman spectroscopy. Raman measurement shows that oxygen vacancies related defects are increasing with Ni substitution which causes changes in optical and PL properties. Local structure measurements show that Ni substitution leads to oxygen vacancies which does change host lattice structure notably. Ce4+ to Ce3+ conversion increases with Ni substitution.

In this article, we investigated the effect of Cr substitution in place of Fe on the structural, magnetic and transport properties of $\mathrm{Co_2FeSi}$ alloy. A comprehensive structural analysis is done using X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS) spectroscopy. Quaternary Heusler compounds $\mathrm{Co_2Fe_{1-x}Cr_xSi}$ with Cr content (x = 0.1, 0.3, 0.5) were found to crystallize in cubic structure. The synchrotron based EXAFS studies reveal that the anti-site disorder increases with the increase in Cr concentration. The saturation magnetization values in all the alloys are found to be less than those expected from the Slater-Pauling rule, which may be due to the some inherent disorder. A detailed resistivity analysis in the temperature range of 5-300 K is done, taking into account different scattering mechanisms. The residual resistivity ratio is found to decrease with increasing Cr concentration. A disorder induced resistivity minimum due to weak localization effect is seen for x = 0.5. The resistivity measurements also indicate that the half-metallic character survives upto 100 K for x = 0.1, whereas the alloys with x= 0.3 and 0.5 show signature of half- metallic nature even at higher temperatures. First principles calculation done with a more robust exchange correlation functional (namely HSE-06) confirms the half metallicity in the entire concentration range. Theoretically simulated band gap and magnetic moments compliment the experimental findings and are compared wherever possible. All these properties make $\mathrm{Co_2Fe_{1-x}Cr_xSi}$ a promising material for spintronics application., Comment: 14 pages, 9 figures

Among multifunctional double perovskite oxides, La2NiMnO6 has recently drawn significant attention due to its importance both in terms of understanding of fundamental physics and potential for device applications. The relative alteration in Ni:Mn ratio strongly influences the structural, magnetic and electrical properties of La2NiMnO6. In the present study, La2Ni1.5Mn0.5O6 sample with Ni:Mn = 3:1 ratio has been prepared using sol-gel method and modifications of the above physical properties from that of a stiochiometric sample of La2NiMnO6 are discussed. Rietveld analysis of X-ray diffraction data shows that La2Ni1.5Mn0.5O6 samples belong to a major monoclinic structural phase (P21/n), with partially disordered arrangement of Ni and Mn ions. Magnetic characterization and X-ray absorption near edge structure analysis reveal ferromagnetic ordering around 260 K between Ni2+ and Mn4+ ions while spin glass like behavior is revealed at ~45 K. Semiconducting nature of the sample was confirmed from temperature dependent resistivity measurements as well as from the I-V characteristics measured at room temperature., Comment: 10

We present high-resolution Giant Metrewave Radio Telescope (GMRT) HI observations and deep Canada-France-Hawaii Telescope (CFHT) optical imaging of two galaxy groups: NGC 4725/47 and NGC 3166/9. These data are part of a multi-wavelength unbiased survey of the gas-rich dwarf galaxy populations in three nearby interacting galaxy groups. The NGC 4725/47 group hosts two tidal knots and one dIrr. Both tidal knots are located within a prominent HI tidal tail, appear to have sufficient mass (M_gas~10^8 M_sol) to evolve into long-lived tidal dwarf galaxies (TDGs) and are fairly young in age. The NGC 3166/9 group contains a TDG candidate, AGC 208457, at least three dIrrs and four HI knots. Deep CFHT imaging confirms that the optical component of AGC 208457 is bluer -- with a 0.28 mag g-r colour -- and a few Gyr younger than its purported parent galaxies. Combining the results for these groups with those from the NGC 871/6/7 group reported earlier, we find that the HI properties, estimated stellar ages and baryonic content of the gas-rich dwarfs clearly distinguish tidal features from their classical counterparts. We optimistically identify four potentially long-lived tidal objects associated to three separate pairs of interacting galaxies, implying that TDGs are not readily produced during interaction events as suggested by some recent simulations. The tidal objects examined in this survey also appear to have a wider variety of properties than TDGs of similar mass formed in current simulations of interacting galaxies, which could be the result of pre- or post-formation environmental influences., Comment: 18 pages, 14 figures, accepted for publication in MNRAS

We present GMRT HI observations and deep CFHT MegaCam optical images of the gas-rich interacting galaxy group NGC 871/NGC 876/NGC 877 (hereafter NGC 871/6/7). Our high-resolution data sets provide a census of the HI and stellar properties of the detected gas-rich group members. In addition to a handful of spiral, irregular and dwarf galaxies, this group harbours an intriguing HI feature, AGC 749170, that has a gas mass of ~10^9.3 M_sol, a dynamical-to-gas mass ratio of ~1 (assuming the cloud is rotating and in dynamical equilibrium) and no optical counterpart in previous imaging. Our observations have revealed a faint feature in the CFHT g'- and r'-bands; if it is physically associated with AGC 749170, the latter has M/L_g > 1000 M_sol/L_sol as well as a higher metallicity (estimated using photometric colours) and a significantly younger stellar population than the other low-mass gas-rich group members. These properties, as well as its spectral and spatial location with respect to its suspected parent galaxies, strongly indicate a tidal origin for AGC 749170. Overall, the HI properties of AGC 749170 resemble those of other optically dark/dim clouds that have been found in groups. These clouds could represent a class of relatively long-lived HI-rich tidal remnants that survive in intermediate-density environments., Comment: 11 pages, 6 figures, accepted for publication in MNRAS

We study conditions under which carbon clusters of different sizes form and stabilize. {We describe an approach to equilibrium by simulating tenuous carbon gas dynamics to long times.} First, we use reactive molecular dynamics simulations to describe the nucleation of long chains, large clusters, and complex cage structures in carbon and hydrogen rich interstellar gas phases. We study how temperature, particle density, presence of hydrogen, and carbon inflow affect the nucleation of molecular moieties with different characteristics, in accordance with astrophysical conditions. We extend the simulations to densities which are orders of magnitude lower than current laboratory densities, to temperatures relevant to circumstellar environments of planetary nebulae, and to longtime (microsecond) formation timescales. We correlate cluster size distributions from dynamical simulations with thermodynamic equilibrium intuitions, where at low temperatures and gas densities, entropy plays a significant role., Comment: 5 figures (accepted for ApJ)