To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Interactions play a significant role in the formation and evolution of galaxies in the Universe. The galaxy systems, NGC 7252 and NGC 5291, are two nearby interacting systems that are hosting tidal dwarf galaxies (TDGs) and star-forming knots. The present work aims (a) to determine the attenuation-corrected star formation rate (SFR) of the interacting system NGC 7252, (b) to compare the star formation in the NGC 7252 system with that of the NGC 5291 system, and (c) to explore the relation between surface densities of gas and SFR in these two systems. The study utilises high-resolution FUV and NUV imaging data from the ultraviolet imaging telescope on board AstroSat. Six star-forming regions, including the merger remnant, were identified in the NGC 7252 system. The SFR corrected for attenuation of the knots in the NGC 7252 system is determined using the continuum slope $\beta$ calculated from the FUV-NUV colour. It has been observed that the attenuation-corrected SFR values of the knots in this system fall within the range of SFR values determined for the NGC 5291 knots. The TDGs in both systems adhere to the same Kennicutt–Schmidt relation as regular spiral galaxies.
Obscuration in active galactic nuclei (AGN) provides valuable insights into the nature of the material surrounding the central engine. Compton-thick AGN (CTAGN), characterised by a column density of $N_{\mathrm{H}} \geq 1.5 \times 10^{24} \ \mathrm{cm}^{-2}$, are heavily obscured by substantial amounts of dust and gas. While X-ray observations are primarily used to determine this column density, our understanding of obscuration properties in the sub-mm regime, particularly for CTAGN, remains limited. In this study, we analyse archival data from the Atacama Large Millimetre/sub-millimetre Array (ALMA) for both CTAGN and non-CTAGN sources, as identified by the 70-month catalogue of the all-sky hard X-ray Swift/Burst Alert Monitor survey and other X-ray surveys. Integrated intensity maps (moment 0) of CO(3–2) emission reveal a concentrated distribution of dense gas around the nucleus. Utilising a constant CO-to-H2 conversion factor, $X_{\mathrm{CO}} = 2.2 \times 10^{20} \ \mathrm{cm}^{-2} \ (\mathrm{K\ km\ s}^{-1})^{-1}$, we find that the derived molecular hydrogen column densities, $N_{\mathrm{H_2}}$, are generally lower than the total hydrogen column densities, $N_{\mathrm{H}}$, obtained from X-ray observations. However, the $N_{\mathrm{H_2}}$ values derived in this work are slightly higher than those reported in previous studies due to the adoption of a higher CO-to-H2 conversion factor. This discrepancy between $N_{\mathrm{H}}$ and $N_{\mathrm{H_2}}$ is consistent with prior findings that X-ray-derived column densities are typically higher, except in the case of non-CTAGN, where $N_{\mathrm{H_2}}$ can exceed $N_{\mathrm{H}}$. Statistical analysis using Kendall and Spearman tests reveals a positive monotonic relationship between $N_{\mathrm{H}}$ and $N_{\mathrm{H_2}}$, although the correlation is not statistically significant. This suggests a complex interplay of factors influencing these properties. The optically thick nature of CO in dense regions may contribute to the observed discrepancies. Our results highlight the importance of adopting an accurate CO-to-H2 conversion factor to derive reliable column densities, which could serve as an alternative method for identifying CTAGN. Further investigations with more comprehensive data sets and refined methodologies are needed to better understand the relationship between sub-millimetre and X-ray properties in AGNs.
The effects of diffraction, reflection, and mutual coupling on the spectral smoothness of radio telescopes become increasingly important at low frequencies, where the observing wavelength may be significant compared with the antenna or array dimensions. These effects are important for both traditional parabolic antennas, which are prone to the ‘standing wave’ phenomenon caused by interference between direct and scattered wavefronts, and aperture arrays, such as the SKA-Low, MWA, HERA, and LOFAR which have more complicated scattering geometries and added dependence on pointing direction (scan angle). Electromagnetic modelling of these effects is computationally intensive and often only possible at coarse frequency resolution. Therefore, using the example of SKA-Low station configurations, we investigate the feasibility of parameterising scattering matrices and separating antenna and array contributions to telescope chromaticity. This allows deeper insights into the effect on spectral smoothness and frequency-dependent beam patterns of differing antenna configurations. Even for the complicated SKA-Low element design, band-limited delay-space techniques appear to produce similar results to brute-force electromagnetic models and allow for faster computation of station beam hypercubes (position, frequency, and polarisation-dependent point spread functions) at arbitrary spectral resolution. As such techniques could facilitate improvements in the design of low-frequency spectral-line surveys, we conduct a simulated Cosmic Dawn experiment using different observing techniques and station configurations.
Polar ring galaxies (PRGs) are a unique class of galaxies characterised by a ring of gas and stars orbiting nearly orthogonal to the main body. This study delves into the evolutionary trajectory of PRGs using the exemplary trio of NGC 3718, NGC 2685, and NGC 4262. We investigate the distinct features of PRGs by analysing their ring and host components to reveal their unique characteristics through spectral energy distribution (SED) fitting. Using CIGALE, we performed SED fitting to independently analyse the ring and host spatially resolved regions, marking the first decomposed SED analysis for PRGs, which examines stellar populations using high-resolution observations from AstroSat UVIT at a resolved scale. The UV-optical surface profiles provide an initial idea that distinct patterns in the galaxies, with differences in FUV and NUV, suggest three distinct stages of ring evolution in the selected galaxies. The study of resolved-scale stellar regions reveals that the ring regions are generally younger than their host galaxies, with the age disparity progressively decreasing along the evolutionary sequence from NGC 3718 to NGC 4262. Star formation rates (SFR) also exhibit a consistent pattern, with higher SFR in the ring of NGC 3718 compared to the others, and a progressive decrease through NGC 2685 and NGC 4262. Finally, the representation of the galaxies in the HI gas fraction versus the NUV–$\text r$ plane supports the idea that they are in three different evolutionary stages of PRG evolution, with NGC 3718 in the initial stage, NGC 2685 in the intermediate stage, and NGC 4262 representing the final stage. This study concludes that PRGs undergo various evolutionary stages, as evidenced by the observed features in the ring and host components. NGC 3718, NGC 2685, and NGC 4262 represent different stages of this evolution, highlighting the dynamic nature of PRGs and emphasising the importance of studying their evolutionary processes to gain insights into galactic formation and evolution.
Using clean numerical simulation (CNS) in which artificial numerical noise is negligible over a finite, sufficiently long interval of time, we provide evidence, for the first time, that artificial numerical noise in direct numerical simulation (DNS) of turbulence is approximately equivalent to thermal fluctuation and/or stochastic environmental noise. This confers physical significance on the artificial numerical noise of DNS of the Navier–Stokes equations. As a result, DNS on a fine mesh should correspond to turbulence under small internal/external physical disturbance, whereas DNS on a sparse mesh corresponds to turbulent flow under large physical disturbance. The key point is that all of them have physical meanings and so are correct in terms of their deterministic physics, even if their statistics are quite different. This is illustrated herein. Our paper provides a positive viewpoint regarding the presence of artificial numerical noise in DNS.
High-beta magnetised plasmas often exhibit anomalously structured temperature profiles, as seen from galaxy cluster observations and recent experiments. It is well known that when such plasmas are collisionless, temperature gradients along the magnetic field can excite whistler waves that efficiently scatter electrons to limit their heat transport. Only recently has it been shown that parallel temperature gradients can excite whistler waves also in collisional plasmas. Here, we develop a Wigner–Moyal theory for the collisional whistler instability starting from Braginskii-like fluid equations in a slab geometry. This formalism is necessary because, for a large region in parameter space, the fastest-growing whistler waves have wavelengths comparable to the background temperature gradients. We find additional damping terms in the expression for the instability growth rate involving inhomogeneous Nernst advection and resistivity. They (i) enable whistler waves to re-arrange the electron temperature profile via growth, propagation and subsequent dissipation, and (ii) allow non-constant temperature profiles to exist stably. For high-beta plasmas, the marginally stable solutions take the form of a temperature staircase along the magnetic field lines. The electron heat flux can also be suppressed by the Ettingshausen effect when the whistler intensity profile is sufficiently peaked and oriented opposite the background temperature gradient. This mechanism allows cold fronts without magnetic draping, might reduce parallel heat losses in inertial fusion experiments and generally demonstrates that whistler waves can regulate transport even in the collisional limit.
A charged particle in a suitably strong magnetic field spirals along the field lines while slowly drifting transversely. This note provides a brief derivation of an effective Lagrangian formulation for the guiding-centre approximation that captures this dynamics without resolving the gyro motion. It also explains how the effective Lagrangian may, for special magnetic fields, admit a ‘quasi-symmetry’ which can give rise to a conserved quantity helpful for plasma confinement in fields lacking a geometric isometry. The aim of this note is to offer a pedagogical introduction and some perspectives on this well-established subject.
We provide a fundamentally new perspective on subcritical turbulence in plasmas, based on coherent structures, which are obtained and characterised via direct numerical solution. The domains where these coherent states exist appear to be closely connected to the those where related turbulent states can exist, so there may be a deep connection between the stability of these coherent structures and the domain where sustained turbulence is possible. In contrast to previous descriptions of turbulence in terms of a stochastic collection of linear waves, we present a fundamentally nonlinear representation based on more general classes of translating oscillatory nonlinear solutions. In turbulent tokamak plasmas, the transport can often be completely suppressed by introducing a background shear flow, whose amplitude is an important control parameter. As this parameter is decreased below a critical value, radially localised structures appear, becoming larger and more complex, in both gyrokinetic simulations and a simpler fluid model of the plasma. For the fluid model, we directly solve for a particular class of nonlinear solutions, relative periodic orbits, and determine their stability, thus explaining why these isolated structures appear in initial-value simulations. The increase of complexity as the flow shear is reduced is explained by a series of Hopf bifurcations of these nonlinear solutions, which we quantify via stability analysis. In gyrokinetic simulations, we are able to indirectly determine the underlying relative periodic orbits by imposing symmetry conditions on the simulations.
Landau damping is one of the cornerstones of plasma physics. Based on the initial-value approach adopted by Landau in his original derivation of Landau damping, we examine the solutions of the linear Vlasov–Poisson system for different equilibrium distribution functions $f_0(v)$, going beyond the traditional focus on the root with largest imaginary part and investigating the full set of roots that the dispersion relation of the system generally admits. Specifically, we provide analytical insights into the number and the structure of the roots for entire and meromorphic functions $f_0(v)$, such as Maxwellian and $\kappa$ distributions, we discuss the potential issues related to the redefinition of $\partial{f}_0/\partial{v}$ as a complex variable function and we show how different sigmoids affect the root structure associated with non-meromorphic cut-off distribution functions. Finally, based on the comparison of the several root structures considered, we wonder if the multiple roots might hint at a deeper understanding of the Landau damping phenomenon.
The change in direction of the wavevector and group velocity experienced by a wave refracted at the interface of an anisotropic medium in uniform linear motion are determined analytically. These transmission conditions, which are shown to be consistent with the generalised Snell’s law written in the laboratory frame, are then used to examine the effect of motion on waves incident on a magnetised plasma. For an incident wave in the plane perpendicular to the magnetic field the motion is observed to lead to non-negligible deviation of the low-frequency X-mode, as well as to non-symmetrical total reflection angles. These effects are shown to be further complicated when the magnetic field is in the plane formed by the incident wavevector and the medium’s velocity, as the anisotropy now competes with the motion-induced drag. Although obtained in simplified configurations, these results suggest that accounting for motion when modelling plasma wave trajectories could be important under certain conditions, calling for a more detailed quantification of the effect of motion in actual diagnostics and plasma control schemes.
We present a quantum algorithm based on repeated measurement to solve initial-value problems for nonlinear ordinary differential equations (ODEs), which may be generated from partial differential equations in plasma physics. We map a dynamical system to a Hamiltonian form, where the Hamiltonian matrix is a function of dynamical variables. To advance in time, we measure expectation values from the previous time step and evaluate the Hamiltonian function classically, which introduces stochasticity into the dynamics. We then perform standard quantum Hamiltonian simulation over a short time, using the evaluated constant Hamiltonian matrix. This approach requires evolving an ensemble of quantum states, which are consumed each step to measure the required observables. We apply this approach to the classic logistic and Lorenz systems, in both integrable and chaotic regimes. Our analysis shows that the solutions’ accuracy is influenced by both the stochastic sampling rate and the nature of the dynamical system.
Dusty plasmas typically contain various species of dust particles, though most studies have focused on homogeneous systems. This paper investigates the propagation of dust acoustic waves in an inhomogeneous dusty plasma with an interface, analysing how plasma inhomogeneity influences wave behaviour. Using scattering and reductive perturbation methods, we show that both transmitted and reflected waves depend strongly on the mass ratio between regions. Dust acoustic waves cannot propagate through a dust lattice when the wavelength is smaller than the lattice constant. At a discontinuous interface, at least one transmitted solitary wave is generated, with its amplitude determined by the mass ratio, while at most one reflected solitary wave can exist. These results underscore the critical role of the mass ratio in wave propagation and suggest a method for estimating dust particle masses and properties by analysing the incident, transmitted and reflected waves.
We study the onset of spontaneous dynamics in the follower force model of an active filament, wherein a slender elastic filament in a viscous liquid is clamped normal to a wall at one end and subjected to a tangential compressive force at the other. Clarke et al. (Phys. Rev. Fluids, vol. 9, 2024, 073101) recently conducted a thorough investigation of this model using methods of computational dynamical systems; inter alia, they showed that the filament first loses stability via a supercritical double-Hopf bifurcation, with periodic ‘planar-beating’ states (unstable) and ‘whirling’ states (stable) simultaneously emerging at the critical follower-force value. We complement their numerical study by carrying out a weakly nonlinear analysis close to this unconventional bifurcation, using the method of multiple scales. The main outcome is an ‘amplitude equation’ governing the slow modulation of small-magnitude oscillations of the filament in that regime. Analysis of this reduced-order model provides insights into the onset of spontaneous dynamics, including the creation of the nonlinear whirling states from particular superpositions of linear planar-beating modes as well as the selection of whirling over planar beating in three-dimensional scenarios.
This paper presents a novel approach for simulating plasma instabilities in tokamak plasmas using the piecewise field-aligned finite element method in combination with the particle-in-cell method. Our method traditionally aligns the computational grid, but defines the basis functions in piecewise field-aligned coordinates to avoid grid deformation while naturally representing the field-aligned mode structures. This scheme is formulated and implemented numerically. It also applies to the unstructured triangular meshes in principle. We have conducted linear benchmark tests, which agree well with previous results and traditional schemes. Furthermore, multiple-$n$ simulations are also carried out as a proof of principle, demonstrating the efficiency of this scheme in nonlinear turbulence simulations within the framework of the finite element method.
We analyse distributions of the spatial scales of coherent intermittent structures – current sheets – obtained from fully kinetic, two-dimensional simulations of relativistic turbulence in a collisionless pair plasma using unsupervised machine-learning data dissection. We find that the distribution functions of sheet length $\ell$ (longest scale of the analysed structure in the direction perpendicular to the dominant guide field) and curvature $r_c$ (radius of a circle fitted to the structures) can be well-approximated by power-law distributions, indicating self-similarity of the structures. The distribution for the sheet width $w$ (shortest scale of the structure) peaks at the kinetic scales and decays exponentially at larger values. The data shows little or no correlation between $w$ and $\ell$, as expected from theoretical considerations. The typical $r_c$ depends linearly on $\ell$, which indicates that the sheets all have a similar curvature relative to their sizes. We find a weak correlation between $r_c$ and $w$. Our results can be used to inform realistic magnetohydrodynamic subgrid models for plasma turbulence in high-energy astrophysics.
A point force acting on a Brinkman fluid in confinement is always counterbalanced by the force on the porous medium, the force on the walls and the stress at open boundaries. We discuss the distribution of those forces in different geometries: a long pipe, a medium with a single no-slip planar boundary, a porous sphere with an open boundary and a porous sphere with a no-slip wall. We determine the forces using the Lorentz reciprocal theorem and additionally validate the results with explicit analytical flow solutions. We discuss the relevance of our findings for cellular processes such as cytoplasmic streaming and centrosome positioning.
Particle-laden horizontal turbulent pipe flow is studied experimentally in the two-way coupling regime with a focus on delineating the effects of particle-to-fluid density ratio $\rho _{p}/\rho _{f}=1$ and 1.05 on the fluid and particle statistics. Particle volume fraction $\phi _{v}$ up to $1\,\%$ and viscous Stokes numbers ranging from $St^+ \approx 1.2$ to $St^+ \approx 3.8$ are investigated at friction Reynolds number $Re_\tau \approx 195$ using time-resolved two-dimensional particle image and tracking velocimetry. Substantial differences are observed between the statistics of neutrally buoyant (i.e. $\rho _{p}/\rho _{f}=1$) and denser (i.e. $\rho _{p}/\rho _{f}=1.05$) settling particles (with settling velocities 0.12–0.32 times the friction velocity), which, at most instances, show opposing trends compared to unladen pipe flow statistics. Neutrally buoyant particles show a slightly increased overall drag and suppressed turbulent stresses, but elevated particle–fluid interaction drag and results in elongated turbulent structures compared to the unladen flow, whereas $\rho _{p}/\rho _{f}=1.05$ particles exhibit a slight overall drag reduction even with increased radial turbulent stresses, and shorter streamwise structures compared to the unladen flow. These differences are enhanced with increasing $St^+$ and $\phi _v$, and can be attributed to the small but non-negligible settling velocity of denser particles, which also leads to differing statistics in the upper and lower pipe halves.
We perform simulations of a two-fluid–structure interaction problem involving liquid–gas flow past a fully submerged stationary circular cylinder. Interactions between the liquid–gas interface with finite surface tension and flow disturbances arising from the cylinder induce a variety of interfacial phenomena and wake structures. We map different interface regimes in a parameter space defined by the Bond number $Bo \in [100, 5000]$ and the submergence depth $h/D \in [1, 2.5]$ of the cylinder while keeping the Reynolds (Re) and Weber (We) numbers fixed at 150 and 1000, respectively. The emerging interface features are classified into three distinct regimes: interfacial waves generated by Strouhal vortices, the entrainment of multi-scale gas bubbles and the reduced deformation state. In the interfacial wave regime, we demonstrate that the frequency of transverse interface fluctuations at a specific streamwise location is identical to the vortex shedding frequency. Additionally, the wavelength of interfacial waves is determined by the size of vortex pairs consisting of alternating Strouhal vortices. In the gas entrainment regime at $ Bo = 1000$, our bubble-size distributions reveal that the entrained bubbles have sizes ranging from one to two orders of magnitude smaller than the cylinder. These multi-scale bubbles are formed primarily through plunging and surfing breakers at $h/D = 2.5$. In contrast, at $h/D = 1$, smaller bubbles initially emerge from the breakup of a gas finger. Over time, some of these bubbles grow in size through coalescence cascades. The influence of $ Re \in [50, 150]$ and $ We \in [700, 1100]$ on gas entrainment is quantified in terms of mean bubble size and count. Lastly, we demonstrate how the deformability of the liquid–gas interface drives the hydrodynamic lift force acting on the cylinder. The net downward lift materializes only in the gas entrainment and reduced deformation regimes due to the broken symmetry of the front stagnation point. While our study focuses on two-dimensional simulations, we also provide insights into the three-dimensional gas entrainment mechanism for one of the extreme cases at $h/D = 1$.
Raman fiber lasers, known for their capacity to provide both high-power and precise wavelength emissions, are gaining attraction across a spectrum of applications, including fiber optic communications, sensing, spectroscopy and imaging. However, the scalability of Raman laser power is impeded by the constraints of pump brightness and the deleterious effects of second-order Raman scattering. In this research, we have undertaken a comprehensive experimental and simulation-based investigation into the impact of pump brightness on the output characteristics within an amplifier framework. Our innovative approach integrates high-brightness pumping with multi-mode graded-index fibers. Notably, we have pioneered the introduction of multi-wavelength seed light to facilitate four-wave mixing, thereby effectively mitigating higher-order Raman scattering. This novel strategy has culminated in the achievement of a 4 kW Raman laser output in an all-fiber configuration, representing the highest output power reported so far.