-
Extending the Scope of Inference About Predictive Ability to Machine Learning Methods
Authors:
Juan Carlos Escanciano,
Ricardo Parra
Abstract:
Though out-of-sample forecast evaluation is systematically employed with modern machine learning methods and there exists a well-established classic inference theory for predictive ability, see, e.g., West (1996, Asymptotic Inference About Predictive Ability, Econometrica, 64, 1067-1084), such theory is not directly applicable to modern machine learners such as the Lasso in the high dimensional se…
▽ More
Though out-of-sample forecast evaluation is systematically employed with modern machine learning methods and there exists a well-established classic inference theory for predictive ability, see, e.g., West (1996, Asymptotic Inference About Predictive Ability, Econometrica, 64, 1067-1084), such theory is not directly applicable to modern machine learners such as the Lasso in the high dimensional setting. We investigate under which conditions such extensions are possible. Two key properties for standard out-of-sample asymptotic inference to be valid with machine learning are (i) a zero-mean condition for the score of the prediction loss function; and (ii) a fast rate of convergence for the machine learner. Monte Carlo simulations confirm our theoretical findings. We recommend a small out-of-sample vs in-sample size ratio for accurate finite sample inferences with machine learning. We illustrate the wide applicability of our results with a new out-of-sample test for the Martingale Difference Hypothesis (MDH). We obtain the asymptotic null distribution of our test and use it to evaluate the MDH of some major exchange rates at daily and higher frequencies.
△ Less
Submitted 16 April, 2024; v1 submitted 20 February, 2024;
originally announced February 2024.
-
Machine Learning based Pointing Models for Radio/Sub-millimeter Telescopes
Authors:
Bendik Nyheim,
Signe Riemer-Sørensen,
Rodrigo Parra,
Claudia Cicone
Abstract:
Radio, sub-millimiter and millimeter ground-based telescopes are powerful instruments for studying the gas and dust-rich regions of the Universe that are invisible at optical wavelengths, but the pointing accuracy is crucial for obtaining high-quality data. Pointing errors are small deviations of the telescope's orientation from its desired direction. The telescopes use linear regression pointing…
▽ More
Radio, sub-millimiter and millimeter ground-based telescopes are powerful instruments for studying the gas and dust-rich regions of the Universe that are invisible at optical wavelengths, but the pointing accuracy is crucial for obtaining high-quality data. Pointing errors are small deviations of the telescope's orientation from its desired direction. The telescopes use linear regression pointing models to correct for these errors, taking into account various factors such as weather conditions, telescope mechanical structure, and the target's position in the sky. However, residual pointing errors can still occur due to factors that are hard to model accurately, such as thermal and gravitational deformation and environmental conditions like humidity and wind. Here we present a proof-of-concept for reducing pointing error for the Atacama Pathfinder EXperiment (APEX) telescope in the high-altitude Atacama Desert in Chile based on machine learning. Using historic pointing data from 2022, we trained eXtreme Gradient Boosting (XGBoost) models that reduced the root-mean-square errors (RMSE) for azimuth and elevation (horizontal and vertical angle) pointing corrections by 4.3% and 9.5%, respectively, on hold-out test data. Our results will inform operations of current and future facilities such as the next-generation Atacama Large Aperture Submillimeter Telescope (AtLAST).
△ Less
Submitted 13 February, 2024;
originally announced February 2024.
-
Discrete-time staged progression epidemic models
Authors:
Luis Sanz-Lorenzo,
Rafael Bravo de la Parra
Abstract:
In the Staged Progression (SP) epidemic models, infected individuals are classified into a suitable number of states. The goal of these models is to describe as closely as possible the effect of differences in infectiousness exhibited by individuals going through the different stages. The main objective of this work is to study, from the methodological point of view, the behavior of solutions of t…
▽ More
In the Staged Progression (SP) epidemic models, infected individuals are classified into a suitable number of states. The goal of these models is to describe as closely as possible the effect of differences in infectiousness exhibited by individuals going through the different stages. The main objective of this work is to study, from the methodological point of view, the behavior of solutions of the discrete time SP models without reinfection and with a general incidence function. Besides calculating $\mathcal{R}_{0}$, we find bounds for the epidemic final size, characterize the asymptotic behavior of the infected classes, give results about the final monotonicity of the infected classes, and obtain results regarding the initial dynamics of the prevalence of the disease. Moreover, we incorporate into the model the probability distribution of the number of contacts in order to make the model amenable to study its effect in the dynamics of the disease.
△ Less
Submitted 7 February, 2024;
originally announced February 2024.
-
Non-linear population discrete models with two time scales: re-scaling of part of the slow process
Authors:
Luis Sanz,
Rafael Bravo de la Parra,
Marcos Marvá,
Eva Sánchez
Abstract:
In this work we present a reduction result for discrete time systems with two time scales. In order to be valid, previous results in the field require some strong hypotheses that are difficult to check in practical applications. Roughly speaking, the iterates of a map as well as their differentials must converge uniformly on compact sets. Here, we eliminate the hypothesis of uniform convergence of…
▽ More
In this work we present a reduction result for discrete time systems with two time scales. In order to be valid, previous results in the field require some strong hypotheses that are difficult to check in practical applications. Roughly speaking, the iterates of a map as well as their differentials must converge uniformly on compact sets. Here, we eliminate the hypothesis of uniform convergence of the differentials at no significant cost in the conclusions of the result. This new result is then used to extend to nonlinear cases the reduction of some population discrete models involving processes acting at different time scales. In practical cases, some processes that occur at a fast time scale are often only measured at slow time intervals, notably mortality. For a general class of linear models that include such kind of processes, it has been shown that a more realistic approach requires the re-scaling of those processes to be considered at the fast time scale. We develop the same type of re-scaling in some nonlinear models and prove the corresponding reduction results. We also provide an application to a particular model of a structured population in a two-patch environment.
△ Less
Submitted 7 February, 2024;
originally announced February 2024.
-
A discrete model of competing species sharing a parasite
Authors:
Rafael Bravo de La Parra,
Luis Sanz
Abstract:
In this work we develop a discrete model of competing species affected by a common parasite. We analyze the influence of the fast development of the shared disease on the community dynamics. The model is presented under the form of a two time scales discrete system with four variables. Thus, it becomes analytically tractable with the help of the appropriate reduction method. The 2-dimensional redu…
▽ More
In this work we develop a discrete model of competing species affected by a common parasite. We analyze the influence of the fast development of the shared disease on the community dynamics. The model is presented under the form of a two time scales discrete system with four variables. Thus, it becomes analytically tractable with the help of the appropriate reduction method. The 2-dimensional reduced system, that has the same the asymptotic behaviour of the full model, is a generalization of the Leslie-Gower competition model. It has the unfrequent property in this kind of models of including multiple equilibrium attractors of mixed type. The analysis of the reduced system shows that parasites can completely alter the outcome of competition depending on the parasite's basic reproductive number R0. In some cases, initial conditions decide among several exclusion or coexistence scenarios.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Discrete epidemic models with two time scales
Authors:
Luis Sanz-Lorenzo,
Rafael Bravo de la Parra
Abstract:
The main aim of the work is to present a general class of two time scales discrete-time epidemic models. In the proposed framework the disease dynamics is considered to act on a slower time scale than a second different process that could represent movements between spatial locations, changes of individual activities or behaviours, or others. To include a sufficiently general disease model, we fir…
▽ More
The main aim of the work is to present a general class of two time scales discrete-time epidemic models. In the proposed framework the disease dynamics is considered to act on a slower time scale than a second different process that could represent movements between spatial locations, changes of individual activities or behaviours, or others. To include a sufficiently general disease model, we first build up from first principles a discrete-time Susceptible-Exposed-Infectious-Recovered-Susceptible (SEIRS) model and characterize the eradication or endemicity of the disease with the help of its basic reproduction number R0. Then, we propose a general full model that includes sequentially the two processes at different time scales, and proceed to its analysis through a reduced model. The basic reproduction number R0 of the reduced system gives a good approximation of the R0 of the full model since it serves at analyzing its asymptotic behaviour. As an illustration of the proposed general framework, it is shown that there exist conditions under which a locally endemic disease, considering isolated patches in a metapopulation, can be eradicated globally by establishing the appropriate movements between patches.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
Stochastic matrix metapopulation models with fast migration: re-scaling survival to the fast scale
Authors:
Luis Sanz,
Rafael Bravo de la Parra
Abstract:
In this work we address the analysis of discrete-time models of structured metapopulations subject to environmental stochasticity. Previous works on these models made use of the fact that migrations between the patches can be considered fast with respect to demography (maturation, survival, reproduction) in the population. It was assumed that, within each time step of the model, there are many fas…
▽ More
In this work we address the analysis of discrete-time models of structured metapopulations subject to environmental stochasticity. Previous works on these models made use of the fact that migrations between the patches can be considered fast with respect to demography (maturation, survival, reproduction) in the population. It was assumed that, within each time step of the model, there are many fast migration steps followed by one slow demographic event. This assumption allowed one to apply approximate reduction techniques that eased the model analysis. It is however a questionable issue in some cases since, in particular, individuals can die at any moment of the time step. We propose new non-equivalent models in which we re-scale survival to consider its effect on the fast scale. We propose a more general formulation of the approximate reduction techniques so that they also apply to the proposed new models. We prove that the main asymptotic elements in this kind of stochastic models, the Stochastic Growth Rate (SGR) and the Scaled Logarithmic Variance (SLV), can be related between the original and the reduced systems, so that the analysis of the latter allows us to ascertain the population fate in the first. Then we go on to considering some cases where we illustrate the reduction technique and show the differences between both modelling options. In some cases using one option represents exponential growth, whereas the other yields extinction.
△ Less
Submitted 6 February, 2024;
originally announced February 2024.
-
The coloured mutation class of $\mathbb{A}_n$ quivers
Authors:
Viviana Gubitosi,
Claudio Qureshi,
Rafael Parra
Abstract:
In this paper we give an explicit and pure combinatorial description of the $m$-coloured quivers that appears in the $m$-coloured mutation class of a quiver of type $\mathbb{A}_n$. The $m$-coloured mutation defined by Buan and Thomas in \cite{BT} generalizes the well-known quiver mutation of Fomin and Zelevinsky \cite{FZ}. In particular, our description generalizes a result of Buan and Vatne, \cit…
▽ More
In this paper we give an explicit and pure combinatorial description of the $m$-coloured quivers that appears in the $m$-coloured mutation class of a quiver of type $\mathbb{A}_n$. The $m$-coloured mutation defined by Buan and Thomas in \cite{BT} generalizes the well-known quiver mutation of Fomin and Zelevinsky \cite{FZ}. In particular, our description generalizes a result of Buan and Vatne, \cite{BV}, which we recover when $m=1$.
△ Less
Submitted 9 June, 2023;
originally announced June 2023.
-
Finiteness properties and homological dimensions of Skew Group rings
Authors:
Viviana Gubitosi,
Rafael Parra
Abstract:
Let $G$ be a finite group acting on a ring $R$ and $H$ a subgroup of $G$. In this paper we compare some homological dimensions over the skew group rings $RG$ and $RH$. Moreover, under the assumption that $RG$ is a separable extension over $RH$, we show that the skew group rings $RG$ and $RH$ share some properties such as being $n$-Gorenstein, $n$-perfect, $n$-coherent, $(n,d)$, Ding-Chen or IF-rin…
▽ More
Let $G$ be a finite group acting on a ring $R$ and $H$ a subgroup of $G$. In this paper we compare some homological dimensions over the skew group rings $RG$ and $RH$. Moreover, under the assumption that $RG$ is a separable extension over $RH$, we show that the skew group rings $RG$ and $RH$ share some properties such as being $n$-Gorenstein, $n$-perfect, $n$-coherent, $(n,d)$, Ding-Chen or IF-rings.
△ Less
Submitted 22 July, 2022; v1 submitted 21 July, 2022;
originally announced July 2022.
-
On the K-theory of $\mathbb{Z}$-categories
Authors:
Eugenia Ellis,
Rafael Parra
Abstract:
We establish connections between the concepts of Noetherian, regular coherent, and regular n-coherent categories for Z-linear categories with finitely many objects and the corresponding notions for unital rings. These connections enable us to obtain a negative K-theory vanishing result, a fundamental theorem, and a homotopy invariance result for the K-theory of Z-linear categories.
We establish connections between the concepts of Noetherian, regular coherent, and regular n-coherent categories for Z-linear categories with finitely many objects and the corresponding notions for unital rings. These connections enable us to obtain a negative K-theory vanishing result, a fundamental theorem, and a homotopy invariance result for the K-theory of Z-linear categories.
△ Less
Submitted 31 October, 2023; v1 submitted 11 July, 2022;
originally announced July 2022.
-
Emission from HCN and CH$_3$OH in comets. Onsala 20-m observations and radiative transfer modelling
Authors:
P. Bergman,
M. S. Lerner,
A. O. H. Olofsson,
E. Wirström,
J. H. Black,
P. Bjerkeli,
R. Parra,
K. Torstensson
Abstract:
The aim of this work is to characterize HCN and CH$_3$OH emission from recent comets.
We used the Onsala 20-m telescope to search for millimetre transitions of HCN towards a sample of 11 recent and mostly bright comets in the period December 2016 to November 2019. Also CH$_3$OH was searched for in two comets. The HCN sample includes the interstellar comet 2I/Borisov. For the short-period comet 4…
▽ More
The aim of this work is to characterize HCN and CH$_3$OH emission from recent comets.
We used the Onsala 20-m telescope to search for millimetre transitions of HCN towards a sample of 11 recent and mostly bright comets in the period December 2016 to November 2019. Also CH$_3$OH was searched for in two comets. The HCN sample includes the interstellar comet 2I/Borisov. For the short-period comet 46P/Wirtanen we could monitor the variation of HCN emission over a time span of about one month. We performed radiative transfer modelling for the observed molecular emission by also including time-dependent effects due to the outgassing of molecules.
HCN was detected in 6 comets. Two of these are short-period comets and 4 of them are long-period. Six methanol transitions were detected in 46P/Wirtanen, enabling us to determine the gas kinetic temperature. From the observations, we determined the molecular production rates using time-dependent radiative transfer modelling. For 5 comets, we could determine that the HCN mixing ratios lie near 0.1% using contemporary water production rates, $Q_\mathrm{H_2O}$, taken from other studies. This HCN mixing ratio was also found typical in our monitoring observations of 46P/Wirtanen but here we notice deviations, on a daily time scale, up to 0.2% which could indicate short-time changes in the outgassing activity. From our radiative transfer modelling of cometary comae, we found that time-dependent effects on the HCN level populations are of the order 5-15% when $Q_\mathrm{H_2O}$ is around $2\times 10^{28}\,\mathrm{mol\, s^{-1}}$. The effects may be relatively stronger for comets with lower $Q_\mathrm{H_2O}$. The exact details of the time-dependent effects depend on the amount of neutral and electron collisions, radiative pumping, and molecular parameters such as the spontaneous rate coefficient.
△ Less
Submitted 7 February, 2022;
originally announced February 2022.
-
CONCERTO at APEX: installation and technical commissioning
Authors:
A. Monfardini,
A. Beelen,
A. Benoit,
J. Bounmy,
M. Calvo,
A. Catalano,
J. Goupy,
G. Lagache,
P. Ade,
E. Barria,
M. Bethermin,
O. Bourrion,
G. Bres,
C. De Breuck,
F. -X. Desert,
G. Duvauchelle,
A. Fasano,
T. Fenouillet,
J. Garcia,
G. Garde,
C. Hoarau,
W. Hu,
J. -C. Lambert,
F. Levy-Bertrand,
A. Lundgren
, et al. (19 additional authors not shown)
Abstract:
We describe the deployment and first tests on Sky of CONCERTO, a large field-of-view (18.6arc-min) spectral-imaging instrument. The instrument operates in the range 130-310GHz from the APEX 12-meters telescope located at 5100m a.s.l. on the Chajnantor plateau. Spectra with R=1-300 are obtained using a fast (2.5Hz mechanical frequency) Fourier Transform Spectrometer (FTS), coupled to a continuous d…
▽ More
We describe the deployment and first tests on Sky of CONCERTO, a large field-of-view (18.6arc-min) spectral-imaging instrument. The instrument operates in the range 130-310GHz from the APEX 12-meters telescope located at 5100m a.s.l. on the Chajnantor plateau. Spectra with R=1-300 are obtained using a fast (2.5Hz mechanical frequency) Fourier Transform Spectrometer (FTS), coupled to a continuous dilution cryostat with a base temperature of 60mK. Two 2152-pixels arrays of Lumped Element Kinetic Inductance Detectors (LEKID) are installed in the cryostat that also contains the cold optics and the front-end electronics. CONCERTO, installed in April 2021, generates more than 20k spectra per second during observations. We describe the final development phases, the installation and the first results obtained on Sky.
△ Less
Submitted 26 June, 2021;
originally announced June 2021.
-
The SEDIGISM survey: first data release and overview of the Galactic structure
Authors:
F. Schuller,
J. S. Urquhart,
T. Csengeri,
D. Colombo,
A. Duarte-Cabral,
M. Mattern,
A. Ginsburg,
A. R. Pettitt,
F. Wyrowski,
L. Anderson,
F. Azagra,
P. Barnes,
M. Beltran,
H. Beuther,
S. Billington,
L. Bronfman,
R. Cesaroni,
C. Dobbs,
D. Eden,
M. -Y. Lee,
S. -N. X. Medina,
K. M. Menten,
T. Moore,
F. M. Montenegro-Montes,
S. Ragan
, et al. (35 additional authors not shown)
Abstract:
The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical…
▽ More
The SEDIGISM (Structure, Excitation and Dynamics of the Inner Galactic Interstellar Medium) survey used the APEX telescope to map 84 deg^2 of the Galactic plane between l = -60 deg and l = +31 deg in several molecular transitions, including 13CO(2-1) and C18O(2-1), thus probing the moderately dense (~10^3 cm^-3) component of the interstellar medium. With an angular resolution of 30'' and a typical 1-sigma sensitivity of 0.8-1.0 K at 0.25 km/s velocity resolution, it gives access to a wide range of structures, from individual star-forming clumps to giant molecular clouds and complexes. The coverage includes a good fraction of the first and fourth Galactic quadrants, allowing us to constrain the large scale distribution of cold molecular gas in the inner Galaxy. In this paper we provide an updated overview of the full survey and the data reduction procedures used. We also assess the quality of these data and describe the data products that are being made publicly available as part of this first data release (DR1). We present integrated maps and position-velocity maps of the molecular gas and use these to investigate the correlation between the molecular gas and the large scale structural features of the Milky Way such as the spiral arms, Galactic bar and Galactic centre. We find that approximately 60 per cent of the molecular gas is associated with the spiral arms and these appear as strong intensity peaks in the derived Galactocentric distribution. We also find strong peaks in intensity at specific longitudes that correspond to the Galactic centre and well known star forming complexes, revealing that the 13CO emission is concentrated in a small number of complexes rather than evenly distributed along spiral arms.
△ Less
Submitted 2 December, 2020;
originally announced December 2020.
-
The SEDIGISM survey: Molecular clouds in the inner Galaxy
Authors:
A. Duarte-Cabral,
D. Colombo,
J. S. Urquhart,
A. Ginsburg,
D. Russeil,
F. Schuller,
L. D. Anderson,
P. J. Barnes,
M. T. Beltran,
H. Beuther,
S. Bontemps,
L. Bronfman,
T. Csengeri,
C. L. Dobbs,
D. Eden,
A. Giannetti,
J. Kauffmann,
M. Mattern,
S. -N. X. Medina,
K. M. Menten,
M. -Y. Lee,
A. R. Pettitt,
M. Riener,
A. J. Rigby,
A. Trafficante
, et al. (35 additional authors not shown)
Abstract:
We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the…
▽ More
We use the 13CO(2-1) emission from the SEDIGISM high-resolution spectral-line survey of the inner Galaxy, to extract the molecular cloud population with a large dynamic range in spatial scales, using the SCIMES algorithm. This work compiles a cloud catalogue with a total of 10663 molecular clouds, 10300 of which we were able to assign distances and compute physical properties. We study some of the global properties of clouds using a science sample, consisting of 6664 well resolved sources and for which the distance estimates are reliable. In particular, we compare the scaling relations retrieved from SEDIGISM to those of other surveys, and we explore the properties of clouds with and without high-mass star formation. Our results suggest that there is no single global property of a cloud that determines its ability to form massive stars, although we find combined trends of increasing mass, size, surface density and velocity dispersion for the sub-sample of clouds with ongoing high-mass star formation. We then isolate the most extreme clouds in the SEDIGISM sample (i.e. clouds in the tails of the distributions) to look at their overall Galactic distribution, in search for hints of environmental effects. We find that, for most properties, the Galactic distribution of the most extreme clouds is only marginally different to that of the global cloud population. The Galactic distribution of the largest clouds, the turbulent clouds and the high-mass star-forming clouds are those that deviate most significantly from the global cloud population. We also find that the least dynamically active clouds (with low velocity dispersion or low virial parameter) are situated further afield, mostly in the least populated areas. However, we suspect that part of these trends may be affected by some observational biases, and thus require further follow up work in order to be confirmed.
△ Less
Submitted 2 December, 2020;
originally announced December 2020.
-
K-theory of n-coherent rings
Authors:
Eugenia Ellis,
Rafael Parra
Abstract:
Let $R$ be a strong $n$-coherent ring such that each finitely $n$-presented $R$-module has finite projective dimension. We consider $\mathcal{FP}_{n}(R)$ the full subcategory of $R$-Mod of finitely $n$-presented modules. We prove that $\mathcal{FP}_{n}(R)$ is an exact category, $K_{i}(R) = K_{i}(\mathcal{FP}_{n}(R))$ for every $i\geq 0$ and obtain an expression of $\operatorname{Nil}_{i}(R)$.
Let $R$ be a strong $n$-coherent ring such that each finitely $n$-presented $R$-module has finite projective dimension. We consider $\mathcal{FP}_{n}(R)$ the full subcategory of $R$-Mod of finitely $n$-presented modules. We prove that $\mathcal{FP}_{n}(R)$ is an exact category, $K_{i}(R) = K_{i}(\mathcal{FP}_{n}(R))$ for every $i\geq 0$ and obtain an expression of $\operatorname{Nil}_{i}(R)$.
△ Less
Submitted 6 November, 2020; v1 submitted 20 August, 2020;
originally announced August 2020.
-
Particle Ratios in a Multi Component Non-Ideal Hadron Resonance Gas
Authors:
Rameez Ahmad Parra,
Saeed Uddin,
Waseem Bashir,
Inam-ul Bashir
Abstract:
We have considered formation of a multi component non-ideal hot and dense gas of hadronic resonances in the ultra-relativistic heavy-ion collisions. In the statistical thermal model approach the equation of state (EoS) of the non interacting ideal hadron resonance gas (IHRG) does not incorporate either the attractive part or the short range repulsive part of the baryonic interaction. On the other…
▽ More
We have considered formation of a multi component non-ideal hot and dense gas of hadronic resonances in the ultra-relativistic heavy-ion collisions. In the statistical thermal model approach the equation of state (EoS) of the non interacting ideal hadron resonance gas (IHRG) does not incorporate either the attractive part or the short range repulsive part of the baryonic interaction. On the other hand in the non-ideal hadron resonance gas (NIHRG) model we can incorporate these interactions using the Van der Waals (VDW) type approach. Studies have been made to see its effect on the critical parameters of the quark-hadron phase transition. However, it can also lead to modifications in the calculated relative particle yields. In this paper we have attempted to understand the effect of such Van der Waals type interactions on the relative particle yields and also studied their dependences on the system thermal parameters, such as the temperature and baryon chemical potential (μ_B). We have also taken into account the decay contributions of the heavier resonances. These results on particle ratios are compared with the corresponding results obtained from the point-like i.e. non interacting IHRG model. It is found that the particle ratios get modified by incorporating the Van der Waals type interactions, especially in a baryon rich system which is expected to be formed at lower RHIC energies, SPS energies and in the forthcoming CBM experiments due to high degree of nuclear stopping in these experiments.
△ Less
Submitted 20 August, 2022; v1 submitted 17 December, 2019;
originally announced December 2019.
-
Continuum and Spectral Line Radiation from a Random Clumpy Medium
Authors:
John E. Conway,
Moshe Elitzur,
Rodrigo Parra
Abstract:
We present a formalism for continuum and line emission from random clumpy media together with its application to problems of current interest, including CO spectral lines from ensembles of clouds and radio emission from HII regions, supernovae and star-forming regions. For line emission we find that the effects of clump opacity on observed line ratios can be indistinguishable from variations of in…
▽ More
We present a formalism for continuum and line emission from random clumpy media together with its application to problems of current interest, including CO spectral lines from ensembles of clouds and radio emission from HII regions, supernovae and star-forming regions. For line emission we find that the effects of clump opacity on observed line ratios can be indistinguishable from variations of intrinsic line strengths, adding to the difficulties in determining abundances from line observations. Our formalism is applicable to arbitrary distributions of cloud properties, provided the cloud volume filling factor is small; numerical simulations show it to hold up to filling factors of about 10%. We show that irrespective of the complexity of the cloud ensemble, the radiative effect of clumpiness can be parametrized at each frequency by a single multiplicative correction to the overall optical depth; this multiplier is derived from appropriate averaging over individual cloud properties. Our main finding is that cloud shapes have only a negligible effect on radiation propagation in clumpy media; the results of calculations employing point-like clouds are practically indistinguishable from those for finite-size clouds with arbitrary geometrical shapes.
△ Less
Submitted 22 August, 2018;
originally announced August 2018.
-
Inferring repeat protein energetics from evolutionary information
Authors:
Rocío Espada,
R. Gonzalo Parra,
Thierry Mora,
Aleksandra M. Walczak,
Diego U. Ferreiro
Abstract:
Natural protein sequences contain a record of their history. A common constraint in a given protein family is the ability to fold to specific structures, and it has been shown possible to infer the main native ensemble by analyzing covariations in extant sequences. Still, many natural proteins that fold into the same structural topology show different stabilization energies, and these are often re…
▽ More
Natural protein sequences contain a record of their history. A common constraint in a given protein family is the ability to fold to specific structures, and it has been shown possible to infer the main native ensemble by analyzing covariations in extant sequences. Still, many natural proteins that fold into the same structural topology show different stabilization energies, and these are often related to their physiological behavior. We propose a description for the energetic variation given by sequence modifications in repeat proteins, systems for which the overall problem is simplified by their inherent symmetry. We explicitly account for single amino acid and pair-wise interactions and treat higher order correlations with a single term. We show that the resulting force field can be interpreted with structural detail. We trace the variations in the energetic scores of natural proteins and relate them to their experimental characterization. The resulting energetic force field allows the prediction of the folding free energy change for several mutants, and can be used to generate synthetic sequences that are statistically indistinguishable from the natural counterparts.
△ Less
Submitted 15 March, 2017; v1 submitted 9 March, 2017;
originally announced March 2017.
-
SEDIGISM: Structure, excitation, and dynamics of the inner Galactic interstellar medium
Authors:
F. Schuller,
T. Csengeri,
J. S. Urquhart,
A. Duarte-Cabral,
P. J. Barnes,
A. Giannetti,
A. K. Hernandez,
S. Leurini,
M. Mattern,
S. -N. X. Medina,
C. Agurto,
F. Azagra,
L. D. Anderson,
M. T. Beltrán,
H. Beuther,
S. Bontemps,
L. Bronfman,
C. L. Dobbs,
M. Dumke,
R. Finger,
A. Ginsburg,
E. Gonzalez,
T. Henning,
J. Kauffmann,
F. Mac-Auliffe
, et al. (19 additional authors not shown)
Abstract:
The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combin…
▽ More
The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".
△ Less
Submitted 17 January, 2017;
originally announced January 2017.
-
Protein Repeats from First Principles
Authors:
Pablo Turjanski,
R. Gonzalo Parra,
Rocío Espada,
Verónica Becher,
Diego U. Ferreiro
Abstract:
Some natural proteins display recurrent structural patterns. Despite being highly similar at the tertiary structure level, repetitions within a single repeat protein can be extremely variable at the sequence level. We propose a mathematical definition of a repeat and investigate the occurrences of these in different protein families. We found that long stretches of perfect repetitions are infreque…
▽ More
Some natural proteins display recurrent structural patterns. Despite being highly similar at the tertiary structure level, repetitions within a single repeat protein can be extremely variable at the sequence level. We propose a mathematical definition of a repeat and investigate the occurrences of these in different protein families. We found that long stretches of perfect repetitions are infrequent in individual natural proteins, even for those which are known to fold into structures of recurrent structural motifs. We found that natural repeat proteins are indeed repetitive in their families, exhibiting abundant stretches of 6 amino acids or longer that are perfect repetitions in the reference family. We provide a systematic quantification for this repetitiveness, and show that this form of repetitiveness is not exclusive of repeat proteins, but also occurs in globular domains. A by-product of this work is a fast classifier of proteins into families, which yields likelihood value about a given protein belonging to a given family.
△ Less
Submitted 8 October, 2015;
originally announced October 2015.
-
Capturing coevolutionary signals in repeat proteins
Authors:
Rocío Espada,
R. Gonzalo Parra,
Thierry Mora,
Aleksandra M. Walczak,
Diego Ferreiro
Abstract:
The analysis of correlations of amino acid occurrences in globular proteins has led to the development of statistical tools that can identify native contacts -- portions of the chains that come to close distance in folded structural ensembles. Here we introduce a statistical coupling analysis for repeat proteins -- natural systems for which the identification of domains remains challenging. We sho…
▽ More
The analysis of correlations of amino acid occurrences in globular proteins has led to the development of statistical tools that can identify native contacts -- portions of the chains that come to close distance in folded structural ensembles. Here we introduce a statistical coupling analysis for repeat proteins -- natural systems for which the identification of domains remains challenging. We show that the inherent translational symmetry of repeat protein sequences introduces a strong bias in the pair correlations at precisely the length scale of the repeat-unit. Equalizing for this bias reveals true co-evolutionary signals from which local native-contacts can be identified. Importantly, parameter values obtained for all other interactions are not significantly affected by the equalization. We quantify the robustness of the procedure and assign confidence levels to the interactions, identifying the minimum number of sequences needed to extract evolutionary information in several repeat protein families. The overall procedure can be used to reconstruct the interactions at long distances, identifying the characteristics of the strongest couplings in each family, and can be applied to any system that appears translationally symmetric.
△ Less
Submitted 25 July, 2014;
originally announced July 2014.
-
A time scales approach to coinfection by opportunistic diseases
Authors:
Marcos Marvá,
Ezio Venturino,
Rafael Bravo de la Parra
Abstract:
Traditional biomedical approaches treat diseases in isolation, but the importance of synergistic disease interactions is now recognized. As a first step we present and analyze a simple coinfection model for two diseases affecting simultaneously a population. The host population is affected by the \emph{primary disease}, a long-term infection whose dynamics is described by a SIS model with demograp…
▽ More
Traditional biomedical approaches treat diseases in isolation, but the importance of synergistic disease interactions is now recognized. As a first step we present and analyze a simple coinfection model for two diseases affecting simultaneously a population. The host population is affected by the \emph{primary disease}, a long-term infection whose dynamics is described by a SIS model with demography, which facilitates individuals acquiring a second disease, \emph{secondary (or \emph{opportunistic}) disease}. The secondary disease is instead a short-term infection affecting only the primary-infected individuals. Its dynamics is also represented by a SIS model with no demography. To distinguish between short and long-term infection the complete model is written as a two time scales system. The primary disease acts at the slow time scale while the secondary disease does at the fast one, allowing a dimension reduction of the system and making its analysis tractable. We show that an opportunistic disease outbreak might change drastically the outcome of the primary epidemic process, although it does among the outcomes allowed by the primary disease. We have found situations in which either acting on the opportunistic disease transmission or recovery rates or controlling the susceptible and infected population size allow to eradicate/promote disease endemicity.
△ Less
Submitted 20 April, 2015; v1 submitted 14 March, 2014;
originally announced March 2014.
-
Detecting Repetitions and Periodicities in Proteins by Tiling the Structural Space
Authors:
R. Gonzalo Parra,
Rocío Espada,
Ignacio E. Sánchez,
Manfred J. Sippl,
Diego U. Ferreiro
Abstract:
The notion of energy landscapes provides conceptual tools for understanding the complexities of protein folding and function. Energy Landscape Theory indicates that it is much easier to find sequences that satisfy the "Principle of Minimal Frustration" when the folded structure is symmetric (Wolynes, P. G. Symmetry and the Energy Landscapes of Biomolecules. Proc. Natl. Acad. Sci. U.S.A. 1996, 93,…
▽ More
The notion of energy landscapes provides conceptual tools for understanding the complexities of protein folding and function. Energy Landscape Theory indicates that it is much easier to find sequences that satisfy the "Principle of Minimal Frustration" when the folded structure is symmetric (Wolynes, P. G. Symmetry and the Energy Landscapes of Biomolecules. Proc. Natl. Acad. Sci. U.S.A. 1996, 93, 14249-14255). Similarly, repeats and structural mosaics may be fundamentally related to landscapes with multiple embedded funnels. Here we present analytical tools to detect and compare structural repetitions in protein molecules. By an exhaustive analysis of the distribution of structural repeats using a robust metric we define those portions of a protein molecule that best describe the overall structure as a tessellation of basic units. The patterns produced by such tessellations provide intuitive representations of the repeating regions and their association towards higher order arrangements. We find that some protein architectures can be described as nearly periodic, while in others clear separations between repetitions exist. Since the method is independent of amino acid sequence information we can identify structural units that can be encoded by a variety of distinct amino acid sequences.
△ Less
Submitted 12 June, 2013;
originally announced June 2013.
-
APEX sub-mm monitoring of gamma-ray blazars
Authors:
S. Larsson,
L. Fuhrmann,
A. Weiss,
E. Angelakis,
T. P. Krichbaum,
I. Nestoras,
J. A. Zensus,
M. Axelsson,
D. Nilsson,
F. Ryde,
L. Hjalmarsdotter,
J. Larsson,
A. Lundgren,
F. Mac-Auliffe,
R. Parra,
G. Siringo
Abstract:
So far, no systematic long-term blazar monitoring programs and detailed variability studies exist at sub-mm wavelengths. Here, we present a new sub-mm blazar monitoring program using the APEX 12-m telescope. A sample of about 40 gamma-ray blazars has been monitored since 2007/2008 with the LABOCA bolometer camera at 345 GHz. First light curves, preliminary variability results and a first compariso…
▽ More
So far, no systematic long-term blazar monitoring programs and detailed variability studies exist at sub-mm wavelengths. Here, we present a new sub-mm blazar monitoring program using the APEX 12-m telescope. A sample of about 40 gamma-ray blazars has been monitored since 2007/2008 with the LABOCA bolometer camera at 345 GHz. First light curves, preliminary variability results and a first comparison with the longer cm/mm bands (F-GAMMA program) are presented, demonstrating the extreme variability characteristics of blazars at such short wavelengths.
△ Less
Submitted 17 June, 2012;
originally announced June 2012.
-
Rapid variability of the compact radio sources in Arp220. Evidence for a population of microblazars?
Authors:
F. Batejat,
J. E. Conway,
A. Rushton,
R. Parra,
P. J. Diamond,
C. J. Lonsdale,
C. J. Lonsdale
Abstract:
Context. The two nuclei of the starburst galaxy Arp220 contain multiple compact radio sources previously identified as radio supernovae or supernova remnants. Aims. In order to search for an embedded radio AGN, or other possible exotic objects, we have carried out a program of VLBI monitoring at 6 cm over three epochs each separated by four months. Methods. Combining the new data with existing dat…
▽ More
Context. The two nuclei of the starburst galaxy Arp220 contain multiple compact radio sources previously identified as radio supernovae or supernova remnants. Aims. In order to search for an embedded radio AGN, or other possible exotic objects, we have carried out a program of VLBI monitoring at 6 cm over three epochs each separated by four months. Methods. Combining the new data with existing data at 6 cm and 18 cm (spanning 4 and 12 years respectively) we are able to characterise source flux density variability on a range of time-scales. Additionally we analyse the variability of sources in shape and position. Results. We detect rapid (< 4 months) variability in three sources (W7, W26 and W29). These sources show possible superluminal motion (> 4c) of jet-like features near rapidly varying almost stationary components. These enigmatic sources might be associated with an AGN or a highly beamed microquasar (i.e. microblazar). Other hypotheses include that the apparent variability is intrinsic and is produced by neutron star powered central components within a supernova remnant, by a sequence of several supernovae within super star clusters, or is extrinsic and is produced by Galactic interstellar scintillation of very compact non-varying objects. Conclusions. A microquasar/microblazar origin seems to be the best explanation for the nature of the variable sources in Arp220.
△ Less
Submitted 14 May, 2012;
originally announced May 2012.
-
Resolution of the Compact Radio Continuum Sources in Arp220
Authors:
Fabien Batejat,
John E. Conway,
Rossa Hurley,
Rodrigo Parra,
Philip J. Diamond,
Colin J. Lonsdale,
Carol J. Lonsdale
Abstract:
We present 2 cm and 3.6 cm wavelength very long baseline interferometry images of the compact radio continuum sources in the nearby ultra-luminous infrared galaxy Arp220. Based on their radio spectra and variability properties, we confirm these sources to be a mixture of supernovae (SNe) and supernova remnants (SNRs). Of the 17 detected sources we resolve 7 at both wavelengths. The SNe generally o…
▽ More
We present 2 cm and 3.6 cm wavelength very long baseline interferometry images of the compact radio continuum sources in the nearby ultra-luminous infrared galaxy Arp220. Based on their radio spectra and variability properties, we confirm these sources to be a mixture of supernovae (SNe) and supernova remnants (SNRs). Of the 17 detected sources we resolve 7 at both wavelengths. The SNe generally only have upper size limits. In contrast all the SNRs are resolved with diameters {\geq} 0.27 pc. This size limit is consistent with them having just entered their Sedov phase while embedded in an interstellar medium (ISM) of density 10^4 cm^{-3} . These objects lie on the diameter-luminosity correlation for SNRs (and so also on the diameter-surface brightness relation) and extend these correlations to very small sources. The data are consistent with the relation L {\propto} D^{-9/4}. Revised equipartition arguments adjusted to a magnetic field to relativistic particle energy density ratio of 1% combined with a reasonable synchrotron-emitting volume filling factor of 10% give estimated magnetic field strengths in the SNR shells of ~ 15-50 mG. The SNR shell magnetic fields are unlikely to come from compression of ambient ISM fields and must instead be internally generated. We set an upper limit of 7 mG for the ISM magnetic field. The estimated energy in relativistic particles, 2%-20% of the explosion kinetic energy, is consistent with estimates from models that fit the IR-radio correlation in compact starburst galaxies.
△ Less
Submitted 29 September, 2011;
originally announced September 2011.
-
The Jacobian cocycle and equidistribution towards the Green current
Authors:
Manuel Rodrigo Parra
Abstract:
Given a holomorphic self-map of complex projective space of degree larger than one, we prove that there exists a finite collection of totally invariant algebraic sets with the following property: given any positive closed (1,1)-current of mass 1 with no mass on any element of this family, the sequence of normalized pull-backs of the current converges to the Green current. Under suitable geometric…
▽ More
Given a holomorphic self-map of complex projective space of degree larger than one, we prove that there exists a finite collection of totally invariant algebraic sets with the following property: given any positive closed (1,1)-current of mass 1 with no mass on any element of this family, the sequence of normalized pull-backs of the current converges to the Green current. Under suitable geometric conditions on the collection of totally invariant algebraic sets, we prove a sharper equidistribution result.
△ Less
Submitted 23 March, 2011;
originally announced March 2011.
-
Lelong numbers on projective varieties
Authors:
Manuel Rodrigo Parra
Abstract:
Given a positive closed (1,1)-current $T$ defined on the regular locus of a projective variety $X$ with bounded mass near the singular part of $X$ and $Y$ an irreducible algebraic subset of $X$, we present uniform estimates for the locus inside $Y$ where the Lelong numbers of $T$ are larger than the generic Lelong number of $T$ along $Y$.
Given a positive closed (1,1)-current $T$ defined on the regular locus of a projective variety $X$ with bounded mass near the singular part of $X$ and $Y$ an irreducible algebraic subset of $X$, we present uniform estimates for the locus inside $Y$ where the Lelong numbers of $T$ are larger than the generic Lelong number of $T$ along $Y$.
△ Less
Submitted 23 November, 2010;
originally announced November 2010.
-
Compact source resolution and rapid variability in Arp220
Authors:
Fabien Batejat,
John E. Conway,
Philip J. Diamond,
Rodrigo Parra,
Colin J. Lonsdale,
Carol J. Lonsdale
Abstract:
We present multi-epoch global VLBI observations at 2 cm, 3.6 cm and 6 cm of the compact radio sources in Arp220. We resolve many sources and estimate sizes, expansion velocities and source classes. We find most source properties are consistent with them being either radio supernovae (SNe) or supernova remnants (SNRs). We extend the luminosity-diameter relation for SNRs to very small sources and ar…
▽ More
We present multi-epoch global VLBI observations at 2 cm, 3.6 cm and 6 cm of the compact radio sources in Arp220. We resolve many sources and estimate sizes, expansion velocities and source classes. We find most source properties are consistent with them being either radio supernovae (SNe) or supernova remnants (SNRs). We extend the luminosity-diameter relation for SNRs to very small sources and argue this supports models where shell magnetic fields are internally amplified. We also detect one probable SN/SNR transition object candidate and three highly variable sources with possible superluminal motion (approximately 4c) of jet-like features near rapidly varying almost stationary components. These enigmatic sources, which show similarities to the recently discovered superluminal source in M82, might be associated with an AGN or a new radio source class (e.g. intermediate-mass black holes or micro-blazars).
△ Less
Submitted 17 November, 2010;
originally announced November 2010.
-
G30.79~FIR~10: A gravitationally bound infalling high-Mass star forming clump
Authors:
P. C. Cortes,
R. Parra,
J. R. Cortes,
E. Hardy
Abstract:
Infalling motions were detected and modeled toward this source. A mean infall velocity of 0.5 km/s with an infall mass rate of 5x10^-3 Solar masses per year was obtained. Also, a previously estimated value for the magnetic field strength in the plane of the sky was refined to be 855 micro Gauss which we used to calculate a mass-to-magnetic flux ratio of 1.9, or super-critical. The virial mass from…
▽ More
Infalling motions were detected and modeled toward this source. A mean infall velocity of 0.5 km/s with an infall mass rate of 5x10^-3 Solar masses per year was obtained. Also, a previously estimated value for the magnetic field strength in the plane of the sky was refined to be 855 micro Gauss which we used to calculate a mass-to-magnetic flux ratio of 1.9, or super-critical. The virial mass from turbulent motions was also calculated finding Mvir=563 solar masses, which gives a ratio of Msubmm/Mvir=5.9. Both values strongly suggest that this clump must be in a state of gravitational collapse.
△ Less
Submitted 25 June, 2010;
originally announced June 2010.
-
COLA III. Radio Detection of AGN in Compact Moderate Luminosity Infra-Red Galaxies
Authors:
R. Parra,
J. E. Conway,
S. Aalto,
P. N. Appleton,
R. P. Norris,
Y. M. Pihlstroem,
L. J. Kewley
Abstract:
We present results from 4.8 GHz VLA and Global-VLBI observations of the northern half of the moderate FIR luminosity (median L_IR = 10^11.01 L_Sol) COLA sample of star-forming galaxies. VLBI sources are detected in a high fraction (20/90) of the galaxies observed. The radio luminosities of these cores (~10^21 W/Hz) are too large to be explained by radio supernovae or supernova remnants and we argu…
▽ More
We present results from 4.8 GHz VLA and Global-VLBI observations of the northern half of the moderate FIR luminosity (median L_IR = 10^11.01 L_Sol) COLA sample of star-forming galaxies. VLBI sources are detected in a high fraction (20/90) of the galaxies observed. The radio luminosities of these cores (~10^21 W/Hz) are too large to be explained by radio supernovae or supernova remnants and we argue that they are instead powered by AGN. These sub-parsec scale radio cores are preferentially detected toward galaxies whose VLA maps show bright 100-500 parsec scale nuclear radio components. Since these latter structures tightly follow the FIR to radio-continuum correlation for star-formation we conclude that the AGN powered VLBI sources are associated with compact nuclear starburst environments. The implications for possible starburst-AGN connections are discussed. The detected VLBI sources have a relatively narrow range of radio luminosity consistent with models in which intense compact Eddington-limited starbursts regulate the gas supply onto a central super-massive black hole. The high incidence of AGN radio cores in compact starbursts suggests little or no delay between the starburst phase and the onset of AGN activity.
△ Less
Submitted 14 June, 2010;
originally announced June 2010.
-
Envelopes of commutative rings
Authors:
Rafael Parra,
Manuel Saorin
Abstract:
Given a significative class $F$ of commutative rings, we study the precise conditions under which a commutative ring $R$ has an $F$-envelope. A full answer is obtained when $F$ is the class of fields, semisimple commutative rings or integral domains. When $F$ is the class of Noetherian rings, we give a full answer when the Krull dimension of $R$ is zero and when the envelope is required to be ep…
▽ More
Given a significative class $F$ of commutative rings, we study the precise conditions under which a commutative ring $R$ has an $F$-envelope. A full answer is obtained when $F$ is the class of fields, semisimple commutative rings or integral domains. When $F$ is the class of Noetherian rings, we give a full answer when the Krull dimension of $R$ is zero and when the envelope is required to be epimorphic. The general problem is reduced to identifying the class of non-Noetherian rings having a monomorphic Noetherian envelope, which we conjecture is the empty class.
△ Less
Submitted 23 June, 2009;
originally announced June 2009.
-
The radio spectra of the compact sources in Arp 220: A mixed population of supernovae and supernova remnants
Authors:
Rodrigo Parra,
John E. Conway,
Philip J. Diamond,
Hannah Thrall,
Colin J. Lonsdale,
Carol J. Lonsdale,
Harding E. Smith
Abstract:
We report the first detection at multiple radio wavelengths (13, 6 and 3.6 cm) of the compact sources within both nuclei of the Ultra Luminous Infra-Red Galaxy Arp 220. We present the radio spectra of the 18 detected sources. In just over half of the sources we find that these spectra and other properties are consistent with the standard model of powerful Type IIn supernovae interacting with the…
▽ More
We report the first detection at multiple radio wavelengths (13, 6 and 3.6 cm) of the compact sources within both nuclei of the Ultra Luminous Infra-Red Galaxy Arp 220. We present the radio spectra of the 18 detected sources. In just over half of the sources we find that these spectra and other properties are consistent with the standard model of powerful Type IIn supernovae interacting with their pre-explosion stellar wind. The rate of appearance of new radio sources identified with these supernova events suggests that an unusually large fraction of core collapse supernovae in Arp 220 are highly luminous; possibly implying a radically different stellar initial mass function or stellar evolution compared to galactic disks. Another possible explanation invokes very short (~3 x 10^5 year) intense (~10^3 M_Sol year^-1) star formation episodes with a duty cycle of ~10%. A second group of our detected sources, consisting of the brightest and longest monitored sources at 18 cm do not easily fit the radio supernova model. These sources show a range of spectral indexes from -0.2 to -1.9. We propose that these are young supernova remnants which have just begun interacting with a surrounding ISM with a density between 10^4 and 10^5 cm^-3. One of these sources is probably resolved at 3.6 cm wavelength with a diameter 0.9 pc. In the western nucleus we estimate that the ionized component of the ISM gives rise to foreground free-free absorption with opacity at 18 cm of <0.6 along the majority of lines of sight. Other sources may be affected by absorption with opacity in the range 1 to 2. These values are consistent with previous models as fitted to the radio recombination lines and the continuum spectrum.
△ Less
Submitted 11 December, 2006;
originally announced December 2006.
-
A compact starburst ring traced by clumpy OH megamaser emission
Authors:
R. Parra,
J. E. Conway,
M. Elitzur,
Y. M. Pihlström
Abstract:
We model the OH megamaser emission from the luminous infrared galaxy IIIZw35 as arising from a narrow rotating starburst ring of radius 22 pc enclosing a mass of 7 x 10^6 Msun. We show how both the compact and apparently diffuse maser emission from this ring can arise from a single phase of unsaturated maser clouds amplifying background radio continuum. The masing clouds are estimated to have a…
▽ More
We model the OH megamaser emission from the luminous infrared galaxy IIIZw35 as arising from a narrow rotating starburst ring of radius 22 pc enclosing a mass of 7 x 10^6 Msun. We show how both the compact and apparently diffuse maser emission from this ring can arise from a single phase of unsaturated maser clouds amplifying background radio continuum. The masing clouds are estimated to have a diameter of <0.7 pc and internal velocity dispersion of about 20 km/s. We find that the clouds are neither self-gravitating nor pressure confined but are freely expanding. Their dispersal lifetimes may set the vertical thickness of the ring. For an estimated internal density of 3 x 10^3 cm^-3, cloud masses are of order 24 Msun. The observed spectral features and velocity gradients indicate that the clouds must be outflowing and escaping the nucleus. The cloud mass outflow rate is estimated to be 0.8 Msun/yr, while the star formation rate is about 19 Msun/yr. Associated ionised gas, possibly generated from dissipated clouds, provides free-free absorption along the source axis, explaining the observed East-West asymmetries. We show that the clumpiness of a maser medium can have a dramatic effect on what is observed even in a relatively low gain OH megamaser. Specifically, in IIIZw35 our clumpy maser model naturally explains the large line to continuum ratios, the large 1667MHz:1665MHz line ratios and the wide velocity dispersions seen in the compact maser spots. Other astrophysical masers showing both compact and apparently diffuse emission might be explained by similar clumpy structures.
△ Less
Submitted 19 July, 2005;
originally announced July 2005.
-
A thin ring model for the OH megamaser in IIIZw35
Authors:
R. Parra,
J. E. Conway,
M. Elitzur
Abstract:
We present a model for the OH megamaser emission in the starburst galaxy IIIZw35. The observed diffuse and compact OH maser components in this source are explained by a single phase of unsaturated clumpy gas distributed in a thin ring structure and amplifying background continuum. We emphasize the importance of clumpiness in the OH masing medium, an effect that has not been fully appreciated pre…
▽ More
We present a model for the OH megamaser emission in the starburst galaxy IIIZw35. The observed diffuse and compact OH maser components in this source are explained by a single phase of unsaturated clumpy gas distributed in a thin ring structure and amplifying background continuum. We emphasize the importance of clumpiness in the OH masing medium, an effect that has not been fully appreciated previously.
The model explains why multiple bright spots are seen only at the ring tangents while smoother emission is found elsewhere. Both the observed velocity gradients and the line to continuum ratios around the ring enquire a geometry where most of the seed photons come from a continuum emission which lies outside the OH ring. To explain both the OH and continuum brightness, free-free absorbing gas is required along the ring axis to partially absorb the far side of the ring. It is proposed that the required geometry arises from an inwardly propagating ring of starburst activity.
△ Less
Submitted 29 June, 2004;
originally announced June 2004.
-
Spectral line and continuum radiation propagation in a clumpy medium
Authors:
J. E. Conway,
M. Elitzur,
R. Parra
Abstract:
We discuss the propagation of spectral line and continuum radiation in a clumpy medium and give general expressions for the observed absorption or emission from a cloud population. We show that the affect of the medium clumpiness can usually be characterised by a single number multiplying the mean column opacity. Our result provides a simpler proof and generalisation of the result of Martin et a…
▽ More
We discuss the propagation of spectral line and continuum radiation in a clumpy medium and give general expressions for the observed absorption or emission from a cloud population. We show that the affect of the medium clumpiness can usually be characterised by a single number multiplying the mean column opacity. Our result provides a simpler proof and generalisation of the result of Martin et al (1984). The formalism provides a simple way to understand the effects of clumping on molecular line profiles and ratios, for example how clumping effects the interpretation of 13CO(1-0) to 12CO(1-0) line ratios. It can also be used as a propagation operator in physical models of clumpy media where the incident radiation effects the spectral line emissivity. We are working to extend the formalism to the propagation of masers in a clumpy medium, but in this case there are special difficulties because formal expectation values are not characteristic of observations because they are biased by rare events.
△ Less
Submitted 22 June, 2004;
originally announced June 2004.
