Showing 1–46 of 46 results for author: Malinen, J

Student behaviour and engagement with adaptive exercises on a thermodynamics course

Authors: Matti Harjula, Ville Havu, Inkeri Kontro, Kimmo Kulmala, Jarmo Malinen, Petri Salo

Abstract: A teaching experiment was carried out in a university-level thermodynamics course using adaptive and interactive e-learning material, created in the new Moodle question type Stateful extending the original e-learning platform STACK. The system collects data about the students that is used to algorithmically classify them according to their behaviour in solving problems. It is observed that the cla… ▽ More A teaching experiment was carried out in a university-level thermodynamics course using adaptive and interactive e-learning material, created in the new Moodle question type Stateful extending the original e-learning platform STACK. The system collects data about the students that is used to algorithmically classify them according to their behaviour in solving problems. It is observed that the classification of this data predicts students' success in the other parts of the course for a majority of students. Also, the classification is statistically consistent with Thermodynamic Concept Survey and Maryland Physics Expectation Survey. △ Less

Submitted 7 March, 2024; originally announced March 2024.

MSC Class: 97U50 ACM Class: K.3.1

Gas phase Elemental abundances in Molecular cloudS (GEMS). IV. Observational results and statistical trends

Authors: M. Rodríguez-Baras, A. Fuente, P. Riviére-Marichalar, D. Navarro-Almaida, P. Caselli, M. Gerin, C. Kramer, E. Roueff, V. Wakelam, G. Esplugues, S. García-Burillo, R. Le Gal, S. Spezzano, T. Alonso-Albi, R. Bachiller, S. Cazaux, B. Commercon, J. R. Goicoechea, J. C. Loison, S. P. Treviño-Morales, O. Roncero, I. Jiménez-Serra, J. Laas, A. Hacar, J. Kirk , et al. (11 additional authors not shown)

Abstract: Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S… ▽ More Gas phase Elemental abundances in Molecular CloudS (GEMS) is an IRAM 30m Large Program designed to estimate the S, C, N, and O depletions and gas ionization degree, X(e-), in a set of star-forming filaments of Taurus, Perseus and Orion. Our immediate goal is to build up a complete database of molecular abundances that can serve as an observational basis for estimating X(e-) and the C, O, N, and S depletions through chemical modeling. We observed and derived the abundances of 14 species (13CO, C18O, HCO+, H13CO+, HC18O+, HCN, H13CN, HNC, HCS+, CS, SO, 34SO, H2S, and OCS) in 244 positions, covering the AV 3 to 100 mag, n(H2) a few 10$^{3}$ to 10$^6$ cm$^{-3}$, and Tk 10 to 30 K ranges in these clouds, avoiding protostars, HII regions, and outflows. A statistical analysis is carried out to identify general trends between different species and with physical parameters. Relations between molecules reveal strong linear correlations which define three different families: (1) 13CO and C18O; (2) H13CO+, HC18O+, H13CN, and HNC; and (3) the S-bearing molecules. The abundances of the CO isotopologs increase with the gas kinetic temperature until TK 15 K. For higher temperatures, the abundance remains constant with a scatter of a factor of 3. The abundances of H13CO+, HC18O+, H13CN, and HNC are well correlated with each other, and all of them decrease with molecular hydrogen density, following the law n(H2)$^{-0.8\pm0.2}$. The abundances of S-bearing species also decrease with n(H2) at a rate of (S-bearing/H)gas n(H2)$^{-0.6\pm0.1}$. The abundances of molecules belonging to groups 2 and 3 do not present any clear trend with gas temperature. At scales of molecular clouds, the C18O abundance is the quantity that better correlates with the cloud mass. We discuss the utility of the 13CO/C18O, HCO+/H13CO+, and H13CO+/H13CN abundance ratios as chemical diagnostics of star formation in external galaxies. △ Less

Submitted 25 February, 2021; originally announced February 2021.

Comments: 28 pages, 22 figures. Accepted for publication in A&A

Journal ref: A&A 648, A120 (2021)

Gas-phase Elemental abundances in Molecular cloudS (GEMS) III. Unlocking the CS chemistry: the CS+O reaction

Authors: N. Bulut, O. Roncero, A. Aguado, J. -C. Loison, D. Navarro-Almaida, V. Wakelam, A. Fuente, E. Roueff, R. Le Gal, P. Caselli, M. Gerinm K. M. Hickson, S. Spezzano, P. Riviere-Marichalar, T. Alonso-Albi, R. Bachiller, I. Jimenez-Serra, C. Kramer, B. Tercero, M. Rodriguez-Baras, S. Garcia-Burillo, J. R. Goicoechea, S. P. Treviño-Morales, G. Esplugues, S. Cazaux, B. Commercon , et al. (14 additional authors not shown)

Abstract: CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the el… ▽ More CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 11 pages, 10 figures, accepted in Astronomy & Astrophysics

Journal ref: Astronomy & Astrophysics 646, A5 (2021)

Distributed solution of Laplacian eigenvalue problems

Authors: Antti Hannukainen, Jarmo Malinen, Antti Ojalammi

Abstract: The purpose of this article is to approximately compute the eigenvalues of the symmetric Dirichlet Laplacian within an interval $(0,Λらむだ)$. A novel domain decomposition Ritz method, partition of unity condensed pole interpolation method, is proposed. This method can be used in distributed computing environments where communication is expensive, e.g., in clusters running on cloud computing services or… ▽ More The purpose of this article is to approximately compute the eigenvalues of the symmetric Dirichlet Laplacian within an interval $(0,Λらむだ)$. A novel domain decomposition Ritz method, partition of unity condensed pole interpolation method, is proposed. This method can be used in distributed computing environments where communication is expensive, e.g., in clusters running on cloud computing services or networked workstations. The Ritz space is obtained from local subspaces consistent with a decomposition of the domain into subdomains. These local subspaces are constructed independently of each other, using data only related to the corresponding subdomain. Relative eigenvalue error is analysed. Numerical examples on a cluster of workstations validate the error analysis and the performance of the method. △ Less

Submitted 31 March, 2021; v1 submitted 18 June, 2020; originally announced June 2020.

Comments: 28 pages

MSC Class: 65F15

Gas phase Elemental abundances in Molecular cloudS (GEMS). II. On the quest for the sulphur reservoir in molecular clouds: the $H_{2}S$ case

Authors: D. Navarro-Almaida, R. Le Gal, A. Fuente, P. Rivière-Marichalar, V. Wakelam, S. Cazaux, P. Caselli, Jacob C. Laas, T. Alonso-Albi, J. C. Loison, M. Gerin, C. Kramer, E. Roueff, R. Bachiller, B. Commerçon, R. Friesen, S. García-Burillo, J. R. Goicoechea, B. M. Giuliano, I. Jiménez-Serra, J. M. Kirk, V. Lattanzi, J. Malinen, N. Marcelino, R. Martín-Domènech , et al. (8 additional authors not shown)

Abstract: Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H… ▽ More Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10. △ Less

Submitted 7 April, 2020; originally announced April 2020.

Comments: 23 pages, 16 figures

Journal ref: A&A 637, A39 (2020)

Notes on glottal flow and acoustic inertial effects

Authors: Jarmo Malinen

Abstract: This text is a compilation of some of the notes that the author has written during the development of the low-order model "DICO" [2, 8, 10, 11] for vowel phonation and the even more rudimentary glottal flow model [9] for processing high-speed glottal video data. The following subject matters are covered: (i) Incompressible, laminar, lossless flow models for idealised rectangular and wedge shape vo… ▽ More This text is a compilation of some of the notes that the author has written during the development of the low-order model "DICO" [2, 8, 10, 11] for vowel phonation and the even more rudimentary glottal flow model [9] for processing high-speed glottal video data. The following subject matters are covered: (i) Incompressible, laminar, lossless flow models for idealised rectangular and wedge shape vocal fold geometries. Equations of motion and the pressure distribution are computed in a closed form for each model using the unsteady Bernoulli's theorem; (ii) The assumption of incompressibility and energy loss (i.e., irrecoverable pressure drop) of the airflow in airways (including the glottis) is discussed using steady compressible Bernoulli theorem as the main tool; (iii) Inertia of an uniform waveguide is studied in terms of the low-frequency limit of the the (acoustic) impedance transfer function. It is observed that the inductive loading in the boundary condition sums up with the waveguide inertance in an expected way; (iv) It is shown that an acoustic waveguide, modelled by Webster's lossless equation with Dirichlet boundary condition at the far end, will produce the expected mass inertance of the fluid column as the low-frequency limit of the impedance transfer function. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Numerical modelling of coupled linear dynamical systems

Authors: Juha Kuortti, Jarmo Malinen, Tom Gustafsson

Abstract: Numerical modelling of several coupled passive linear dynamical systems (LDS) is considered. Since such component systems may arise from partial differential equations, transfer function descriptions, lumped systems, measurement data, etc., the first step is to discretise them into finite-dimensional LDSs using, e.g., the finite element method, autoregressive techniques, and interpolation. The fin… ▽ More Numerical modelling of several coupled passive linear dynamical systems (LDS) is considered. Since such component systems may arise from partial differential equations, transfer function descriptions, lumped systems, measurement data, etc., the first step is to discretise them into finite-dimensional LDSs using, e.g., the finite element method, autoregressive techniques, and interpolation. The finite-dimensional component systems may satisfy various types of energy (in)equalities due to passivity that require translation into a common form such as the scattering passive representation. Only then can the component systems be coupled in a desired feedback configuration by computing pairwise Redheffer star products of LDSs. Unfortunately, a straightforward approach may fail due to ill-posedness of feedback loops between component systems. Adversities are particularly likely if some component systems have no energy dissipation at all, and this may happen even if the fully coupled system could be described by a finite-dimensional LDS. An approach is proposed for obtaining the coupled system that is based on passivity preserving regularisation. Two practical examples are given to illuminate the challenges and the proposed methods to overcome them: the Butterworth low-pass filter and the termination of an acoustic waveguide to an irrational impedance. △ Less

Submitted 11 November, 2019; originally announced November 2019.

MSC Class: 93A30; 93B17; 93B40; 93B50

Multi-scale analysis of the Monoceros OB 1 star-forming region: II. Colliding filaments in the Monoceros OB1 molecular cloud

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the… ▽ More We started a multi-scale analysis of G202.3+2.5, an intertwined filamentary region of Monoceros OB1. In Paper I, we examined the distributions of dense cores and protostars and found enhanced star formation (SF) activity in the junction region of the filaments. In this second paper, we aim to unveil the connections between the core and filament evolutions, and between the filament dynamics and the global evolution of the cloud. We characterise the gas dynamics and energy balance using Herschel and WISE observations and molecular tracers observed with the IRAM 30m and TRAO 14m telescopes. The velocity field of the cloud is examined and velocity-coherent structures are put in perspective with the cloud environment. Two main velocity components (VCs) are revealed, well separated in the north and merged around the location of intense N2H+ emission where Paper I found the peak of SF activity. The relative position of the two VCs along the sightline, and the velocity gradient in N2H+ emission imply that the VCs have been undergoing collision for ~10^5 yrs. The dense gas where N2H+ is detected is interpreted as the compressed region between the two filaments, which corresponds to a high mass inflow rate of ~1e-3 Msun/yr and possibly leads to an increase in its SF efficiency. We identify a protostar in the junction region that possibly powers two crossed intermittent outflows. We show that the HII region around the nearby cluster NCG 2264 is still expanding and its role in the collision is examined. However, we cannot rule out the idea that the collision arises mostly from the global collapse of the cloud. The (sub-)filament-scale observables examined in this paper reveal a collision between G202.3+2.5 sub-structures and its probable role in feeding the cores in the junction region. One must now characterise the cloud morphology, its fragmentation, and magnetic field, all at high resolution. △ Less

Submitted 13 September, 2019; originally announced September 2019.

Comments: A&A, in press; IRAM and TRAO data cubes will be available at CDS

Journal ref: A&A 631, A3 (2019)

Multi-scale analysis of the Monoceros OB 1 star-forming region: I. The dense core population

Authors: Julien Montillaud, Mika Juvela, Charlotte Vastel, J. H. He, Tie Liu, Isabelle Ristorcelli, David J. Eden, Sung-ju Kang, Kee-Tae Kim, Patrick M. Koch, Chang Won Lee, Mark G. Rawlings, Mika Saajasto, Patricio Sanhueza, Archana Soam, Sarolta Zahorecz, Dana Alina, Rebeka Bögner, David Cornu, Yasuo Doi, Johanna Malinen, Douglas Marshall, E. R. Micelotta, V. M. Pelkonen, L. V. Tóth , et al. (2 additional authors not shown)

Abstract: Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures.… ▽ More Current theories and models attempt to explain star formation globally, from core scales to giant molecular cloud scales. A multi-scale observational characterisation of an entire molecular complex is necessary to constrain them. We investigate star formation in G202.3+2.5, a ~10x3 pc sub-region of the Monoceros OB1 cloud with a complex morphology harbouring interconnected filamentary structures. We aim to connect the evolution of cores and filaments in G202.3+2.5 with the global evolution of the cloud and to identify the engines of the cloud dynamics. In this first paper, the star formation activity is evaluated by surveying the distributions of dense cores and protostars, and their evolutionary state, as characterised using both infrared observations from the Herschel and WISE telescopes and molecular line observations with the IRAM 30-m telescope. We find ongoing star formation in the whole cloud, with a local peak in star formation activity around the centre of G202.3+2.5 where a chain of massive cores (10-50 Msun) forms a massive ridge (>150 Msun). All evolutionary stages from starless cores to Class II protostars are found in G202.3+2.5, including a possibly starless, large column density (8x10^{22} cm-2), and massive (52 Msun) core. All the core-scale observables examined in this paper point to an enhanced star formation activity centred on the junction between the three main branches of the ramified structure of G202.3+2.5. This suggests that the increased star-formation activity results from the convergence of these branches. To further investigate the origin of this enhancement, it is now necessary to extend the analysis to larger scales, in order to examine the relationship between cores, filaments and their environment. We address these points through the analysis of the dynamics of G202.3+2.5 in a joint paper. △ Less

Submitted 9 September, 2019; originally announced September 2019.

Journal ref: A&A 631, L1 (2019)

SCOPE: SCUBA-2 Continuum Observations of Pre-protostellar Evolution - Survey Description and Compact Source Catalogue

Authors: D. J. Eden, Tie Liu, Kee-Tae Kim, S. -Y. Liu, K. Tatematsu, J. Di Francesco, K. Wang, Y. Wu, M. A. Thompson, G. A. Fuller, Di Li, I. Ristorcelli, Sung-ju Kang, N. Hirano, D. Johnstone, Y. Lin, J. H. He, P. M. Koch, Patricio Sanhueza, S. -L. Qin, Q. Zhang, P. F. Goldsmith, N. J. Evans II, J. Yuan, C. -P. Zhang , et al. (136 additional authors not shown)

Abstract: We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14… ▽ More We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other. △ Less

Submitted 26 February, 2019; originally announced February 2019.

Comments: 16 pages, 13 figures, 1 table. Accepted for publication in MNRAS

Dust polarisation studies on MHD simulations of molecular clouds: methods comparison for the relative orientations analysis

Authors: Elisabetta R. Micelotta, Mika Juvela, Paolo Padoan, Isabelle Ristorcelli, Dana Alina, Johanna Malinen

Abstract: Context. Polarized dust emission can be used as a probe for magnetic fields at different scales. In particular, the analysis of the relative orientation between the density structures and the magnetic field projected on the plane of the sky can provide information on the role of magnetic fields in shaping the structure of molecular clouds. Aims. The relative orientation has been investigated using… ▽ More Context. Polarized dust emission can be used as a probe for magnetic fields at different scales. In particular, the analysis of the relative orientation between the density structures and the magnetic field projected on the plane of the sky can provide information on the role of magnetic fields in shaping the structure of molecular clouds. Aims. The relative orientation has been investigated using different methods. Our goal is to explicitly compare two of them: the Rolling Hough Transform (RHT) and the gradient technique. Methods. We have applied the RHT method and the gradient technique to two regions, identified in surface brightness maps at 353 GHz generated via MHD simulations post-processed using radiative transfer modelling. For both methods we have derived the relative orientation between the magnetic field and the density structures, to which we have applied two different statistics, the histogram of relative orientation (HRO) statistic and the projected Rayleigh statistic (PRS). Results. The samples of pixels selected by each method are substantially different. When the methods are applied to the same pixel selection, the results in terms of relative orientation as a function of column density are consistent between each other. When each method is applied to its own pixel selection, the differences are apparent. In Region 1, the RHT method roughly reproduces the observed trend of the relative orientation becoming more perpendicular for increasing column density, while the gradient method, applied at the same resolution as RHT, gives the opposite trend, with the relative orientation moving towards a more parallel alignment. In Region 2, the situation is reversed. The inconsistent results, which are due to the different pixel selections operated by the methods and to the intrinsic differences between them, provide complementary valuable information. Abridged. △ Less

Submitted 23 October, 2018; originally announced October 2018.

Comments: 15 pages, 9 figures, submitted to A&A, comments welcomed

Journal ref: A&A 647, A121 (2021)

Gas phase Elemental abundances in Molecular cloudS (GEMS) I. The prototypical dark cloud TMC 1

Authors: A. Fuente, D. Navarro, P. Caselli, M. Gerin, C. Kramer, E. Roueff, T. Alonso-Albi, R. Bachiller, S. Cazaux, B. Commercon, R. Friesen, S. García-Burillo, B. M. Giuliano, J. R. Goicoechea, P. Gratier, A. Hacar, I. Jiménez-Serra, J. Kirk, V. Lattanzi, J. C. Loison, P. J. Malinen, N. Marcelino, R. Martín-Doménech, G. Muñoz-Caro, J. Pineda , et al. (8 additional authors not shown)

Abstract: GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the frac… ▽ More GEMS is an IRAM 30m Large Program whose aim is determining the elemental depletions and the ionization fraction in a set of prototypical star-forming regions. This paper presents the first results from the prototypical dark cloud TMC 1. Extensive millimeter observations have been carried out with the IRAM 30m telescope (3mm and 2mm) and the 40m Yebes telescope (1.3cm and 7mm) to determine the fractional abundances of CO, HCO+, HCN, CS, SO, HCS+, and N2H+ in three cuts which intersect the dense filament at the well-known positions TMC 1-CP, TMC 1-NH3, and TMC 1-C, covering a visual extinction range from Av~3 to ~20mag. Two phases with differentiated chemistry can be distinguished: i) the translucent envelope with molecular hydrogen densities of (1-5)x10$^3$ cm$^{-3}$; and ii) the dense phase, located at Av>10mag, with molecular hydrogen densities >10$^4$ cm$^{-3}$. Observations and modeling show that the gas phase abundances of C and O progressively decrease along the C+/C/CO transition zone where C/H~8x10$^{-5}$ and C/O~0.8-1, until the beginning of the dense phase at Av~10 mag. This is consistent with the grain temperatures being below the CO evaporation temperature in this region. In the case of sulfur, a strong depletion should occur before the translucent phase where we estimate a S/H~(0.4 - 2.2) x10$^{-6}$, an abundance ~7-40 times lower than the solar value. A second strong depletion must be present during the formation of the thick icy mantles to achieve the values of S/H measured in the dense cold cores (S/H~8x10$^{-8}$). Based on our chemical modeling, we constrain the value of $ζぜーた_{\rm H_2}$ to ~(0.5 - 1.8) x10$^{-16}$ s$^{-1}$ in the translucent cloud. △ Less

Submitted 7 March, 2019; v1 submitted 13 September, 2018; originally announced September 2018.

Comments: 26 pages, 20 figures, A&A in press

Dust spectrum and polarisation at 850 um in the massive IRDC G035.39-00.33

Authors: Mika Juvela, Vincent Guillet, Tie Liu, Isabelle Ristorcelli, Veli-Matti Pelkonen, Dana Alina, Leonardo Bronfman, David J. Eden, Kee Tae Kim, Patrick M. Koch, Woojin Kwon, Chang Won Lee, Johanna Malinen, Elisabetta Micelotta, Julien Montillaud, Mark G. Rawlings, Patricio Sanhueza, Archana Soam, Alessio Traficante, Nathalie Ysard, Chuan-Peng Zhang

Abstract: The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation… ▽ More The dust sub-millimetre polarisation of star-forming clouds carries information on dust and the role of magnetic fields in cloud evolution. With observations of a dense filamentary cloud G035.39-00.33, we aim to characterise the dust emission properties and the variations of the polarisation fraction. JCMT SCUBA-2/POL-2 data at 850um are combined with Planck 850um (353GHz) data to map polarisation fractions. With previous SCUBA-2 observations (450um and 850um) and Herschel data, the column densities are determined via modified blackbody fits and via radiative transfer modelling. Models are constructed to examine how the polarisation angles and fractions depend on potential magnetic field geometries and grain alignment. POL-2 data show clear changes in the magnetic field orientation. The filament has a peak column density of N(H2)~7 10^22 cm-2, a minimum dust temperature of T~12 K, and a mass of some 4300Msun for the area N(H2)> 5 10^21 cm-2. The estimated average value of the dust opacity spectral index is beta ~ 1.9. The ratio of sub-millimetre and J band optical depths is tau(250 um)/tau(J) ~ 2.5 10^-3, more than four times the typical values for diffuse medium. The polarisation fraction decreases as a function of column density to p ~ 1% in the central filament. Because of noise, the observed decrease of p(N) is significant only at N(H2)>2 10^22 cm-2. The observations suggest that the grain alignment is not constant. Although the data can be explained with a complete loss of alignment at densities above ~ 10^4 cm-3 or using the predictions of radiative torques alignment, the uncertainty of the field geometry and the spatial filtering of the SCUBA-2 data prevent strong conclusions. G035.39-00.33 shows strong signs of dust evolution and the low polarisation fraction is suggestive of a loss of polarised emission from its densest parts. △ Less

Submitted 6 September, 2018; v1 submitted 4 September, 2018; originally announced September 2018.

Comments: accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 620, A26 (2018)

Efficient solution of symmetric eigenvalue problems from families of coupled systems

Authors: Antti Hannukainen, Jarmo Malinen, Antti Ojalammi

Abstract: Efficient solution of the lowest eigenmodes is studied for a family of related eigenvalue problems with common $2\times 2$ block structure. It is assumed that the upper diagonal block varies between different versions while the lower diagonal block and the range of the coupling blocks remains unchanged. Such block structure naturally arises when studying the effect of a subsystem to the eigenmodes… ▽ More Efficient solution of the lowest eigenmodes is studied for a family of related eigenvalue problems with common $2\times 2$ block structure. It is assumed that the upper diagonal block varies between different versions while the lower diagonal block and the range of the coupling blocks remains unchanged. Such block structure naturally arises when studying the effect of a subsystem to the eigenmodes of the full system. The proposed method is based on interpolation of the resolvent function after some of its singularities have been removed by a spectral projection. Singular value decomposition can be used to further reduce the dimension of the computational problem. Error analysis of the method indicates exponential convergence with respect to the number of interpolation points. Theoretical results are illustrated by two numerical examples related to finite element discretisation of the Laplace operator. △ Less

Submitted 18 June, 2020; v1 submitted 19 June, 2018; originally announced June 2018.

MSC Class: 65F15

A holistic perspective on the dynamics of G035.39-00.33: the interplay between gas and magnetic fields

Authors: Tie Liu, Pak Shing Li, Mika Juvela, Kee-Tae Kim, Neal J. Evans II, James Di Francesco, Sheng-Yuan Liu, Jinghua Yuan, Ken'ichi Tatematsu, Qizhou Zhang, Derek Ward-Thompson, Gary Fuller, Paul F. Goldsmith, P. M. Koch, Patricio Sanhueza, I. Ristorcelli, Sung-ju Kang, Huei-Ru Chen, N. Hirano, Yuefang Wu, Vlas Sokolov, Chang Won Lee, Glenn J. White, Ke Wang, David Eden , et al. (32 additional authors not shown)

Abstract: Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at… ▽ More Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $μみゅー$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament. △ Less

Submitted 18 June, 2018; v1 submitted 26 March, 2018; originally announced March 2018.

Comments: Published on ApJ, 27 pages

Journal ref: Liu, T., Li, P. S., Juvela M., et al., 2018, ApJ, 859 ,151

Near-infrared scattering as a dust diagnostic

Authors: Mika Saajasto, Mika Juvela, Johanna Malinen

Abstract: We examine the possibility of using near-infrared scattering to constrain the local radiation field and the dust properties, for example, the size distribution of the grains, and maximum grain size. We use radiative transfer modelling to examine the constraints provided by J, H, and K bands in combination with mid-infrared surface brightness at 3.6 $μみゅー$m. We use a spherical one-dimensional and elli… ▽ More We examine the possibility of using near-infrared scattering to constrain the local radiation field and the dust properties, for example, the size distribution of the grains, and maximum grain size. We use radiative transfer modelling to examine the constraints provided by J, H, and K bands in combination with mid-infrared surface brightness at 3.6 $μみゅー$m. We use a spherical one-dimensional and elliptical three-dimensional cloud models to study the observable effects of different grain size distributions with varying absorption and scattering properties. As an example, we analyse observations of a molecular cloud in Taurus, TMC-1N. The observed surface brightness ratios between the bands change when the dust properties are changed. However, even a small change of 10% in the surface brightness of one channel changes the estimated powerlaw exponent of the size distribution $γがんま$ by up to ~30% and the estimated strength of the radiation field K_ISRF by up to ~60%. The maximum grain size A_max and $γがんま$ are always strongly anti-correlated. For example, overestimating the surface brightness by 10% changes the estimated radiation field strength by ~20% and the exponent of the size distribution by ~15%. The analysis of our synthetic observations indicates that the relative uncertainty of the parameter distributions are on average A_max, $γがんま$ ~25%, and the deviation between the estimated and correct values $Δでるた$ Q <15%. We show that scattered infrared light can be used to derive meaningful limits for the dust parameters. However, errors in the surface brightness data can result in considerable uncertainties in the derived parameters. △ Less

Submitted 21 February, 2018; originally announced February 2018.

Comments: 19 pages, 28 figures. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 614, A95 (2018)

Galactic cold cores IX. Column density structures and radiative transfer modelling

Authors: M. Juvela, J. Malinen, J. Montillaud, V. -M. Pelkonen, I. Ristorcelli, L. V. Toth

Abstract: The Galactic Cold Cores (GCC) project has made Herschel observations of interstellar clouds where Planck detected compact sources of cold dust emission. Our aim is to characterise the structure of the clumps and their parent clouds. We also examine the accuracy to which the structure of dense clumps can be determined from submillimetre data. We use standard statistical methods to characterise the… ▽ More The Galactic Cold Cores (GCC) project has made Herschel observations of interstellar clouds where Planck detected compact sources of cold dust emission. Our aim is to characterise the structure of the clumps and their parent clouds. We also examine the accuracy to which the structure of dense clumps can be determined from submillimetre data. We use standard statistical methods to characterise the GCC fields. Clumps are extracted using column density thresholding and we construct for each field a three-dimensional radiative transfer (RT) model. These are used to estimate the relative radiation field intensities, clump stability, and the uncertainty of column density estimates. We examine the radial column density profiles of the clumps. In the GCC fields, the structure noise follows the relations previously established at larger scales. The fractal dimension has no significant dependence on column density and the values D = 1.25 +- 0.07 are only slightly lower than in typical molecular clouds. The column density PDFs exhibit large variations, e.g. in the case of externally compressed clouds. At scales r>0.1 pc, the radial column density distributions of the clouds follow an average relation of N~r^{-1}. In spite of a great variety of clump morphology, clumps tend to follow a similar N~r^{-1} relation below r~0.1 pc. RT calculations indicate only factor of 2.5 variation in the local radiation field intensity. The fraction of gravitationally bound clumps increases significantly in regions with A_V > 5 mag but most bound objects appear to be pressure-confined. The GCC host clouds have statistical properties similar to general molecular clouds. The gravitational stability, peak column density, and clump orientation are connected to the cloud background while most other statistics (e.g. D and radial profiles) are insensitive to the environment. △ Less

Submitted 15 January, 2018; v1 submitted 8 January, 2018; originally announced January 2018.

Comments: accepted to A&A, 21 pages (37 pages with appendices, this version including language corrections)

Journal ref: A&A 614, A83 (2018)

Herschel and SCUBA-2 observations of dust emission in a sample of Planck cold clumps

Authors: Mika Juvela, Jinhua He, Katherine Pattle, Tie Liu, George Bendo, David J. Eden, Orsolya Feher, Michel Fich, Gary Fuller, Naomi Hirano, Kee-Tae Kim, Di Li, Sheng-Yuan Liu, Johanna Malinen, Douglas J. Marshall, Deborah Paradis, Harriet Parsons, Veli-Matti Pelkonen, Mark G. Rawlings, Isabelle Ristorcelli, Manash R. Samal, Ken'ichi Tatematsu, Mark Thompson, Alessio Traficante, Ke Wang , et al. (4 additional authors not shown)

Abstract: Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs… ▽ More Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths. △ Less

Submitted 27 December, 2017; v1 submitted 26 November, 2017; originally announced November 2017.

Comments: 37 pages, accepted to A&A, final version with language corrections

Journal ref: A&A 612, A71 (2018)

Correlation of gas dynamics and dust in the evolved filament G82.65-02.00

Authors: M. Saajasto, M. Juvela, K. Dobashi, T. Shimoikura, I. Ristorcelli, J. Montillaud, D. J. Marshall, J. Malinen, V. -M. Pelkonen, O. Fehér, A. Rivera-Ingraham, L. V. Toth, L. Montier, J-Ph. Bernard, T. Onishi

Abstract: The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first The combination of line and continuum observations can provide vital insights to the… ▽ More The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first The combination of line and continuum observations can provide vital insights to the formation and fragmentation of filaments and the initial conditions for star formation. We have carried out line observations to map the kinematics of an evolved, actively star forming filament G82.65-2.00. The filament was first identified from the Planck data as a region of particularly cold dust emission and was mapped at 100-500 $μみゅー$m as a part of Herschel key program Galactic Cold Cores. The observations reveal several velocity components in the field, with strongest line emission concentrated to velocity range $\sim [3,5]$ km $\rm s^{-1}$. The column density of molecular hydrogen along the filament varies from 1.0 to 2.3 $\times 10^{22}$ $\rm cm^{2}$. We have examined six cold clumps from the central part of the filament. The clumps have masses in the range $10 - 20$ $M_{\odot}$ ($\sim70M_{\odot}$ in total) and are close to or above the virial mass. Furthermore, the main filament is heavily fragmented and most of the the substructures have a mass lower than or close to the virial mass, suggesting that the filament is dispersing as a whole. Position-velocity maps of $\rm ^{12}CO$ and $\rm ^{13}CO$ lines indicate that at least one of the striations is kinematically connected to two of the clumps, potentially indicating mass accretion from the striation onto the main filament. We tentatively estimate the accretion rate to be $\dot{M}$ = $2.23 \times 10^{-6}$ $ M_{\odot} / \rm year$. △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: Accepted for publication in A&A

Journal ref: A&A 608, A21 (2017)

An acoustic glottal source for vocal tract physical models

Authors: Antti Hannukainen, Juha Kuortti, Jarmo Malinen, Antti Ojalammi

Abstract: A sound source was proposed for acoustic measurements of physical models of the human vocal tract. The physical models are produced by Fast Prototyping, based on Magnetic Resonance Imaging during prolonged vowel production. The sound source, accompanied by custom signal processing algorithms, is used for two kinds of measurements: (i) amplitude frequency response and resonant frequency measurement… ▽ More A sound source was proposed for acoustic measurements of physical models of the human vocal tract. The physical models are produced by Fast Prototyping, based on Magnetic Resonance Imaging during prolonged vowel production. The sound source, accompanied by custom signal processing algorithms, is used for two kinds of measurements: (i) amplitude frequency response and resonant frequency measurements of physical models, and (ii) signal reconstructions at the source output according to a target waveform with measurements at the mouth position of the physical model. The proposed source and the software are validated by measurements on a physical model of the vocal tract corresponding to vowel [a] of a male speaker. △ Less

Submitted 26 April, 2017; v1 submitted 5 April, 2017; originally announced April 2017.

Comments: 27 pages

MSC Class: 92C55; 93C95; 74J25

Galactic Cold Cores. VIII. Filament formation and evolution: Filament properties in context with evolutionary models

Authors: A. Rivera-Ingraham, I. Ristorcelli, M. Juvela, J. Montillaud, A. Men'shchikov, J. Malinen, V. -M. Pelkonen, A. Marston, P. G. Martin, L. Pagani, R. Paladini, D. Paradis, N. Ysard, D. Ward-Thompson, J. -P. Bernard, D. J. Marshall, L. Montier, L. V. Toth

Abstract: Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two… ▽ More Filaments are key for star formation models. As part of the study carried out by the Herschel GCC Programme, here we study the filament properties presented in GCC.VII in context with theoretical models of filament formation and evolution. A conservative sample of filaments at a distance D<500pc was extracted with the Getfilaments algorithm. Their physical structure was quantified according to two main components: the central (Gaussian) region (core component), and the power-law like region dominating the filament column density profile at larger radii (wing component). The properties and behaviour of these components relative to the total linear mass density of the filament and its environmental column density were compared with theoretical models describing the evolution of filaments under gravity-dominated conditions. The feasibility of a transition to supercritical state by accretion is dependent on the combined effect of filament intrinsic properties and environmental conditions. Reasonably self-gravitating (high Mline-core) filaments in dense environments (av\sim3mag) can become supercritical in timescales of t\sim1Myr by accreting mass at constant or decreasing width. The trend of increasing Mline-tot (Mline-core and Mline-wing), and ridge Av with background also indicates that the precursors of star-forming filaments evolve coevally with their environment. The simultaneous increase of environment and filament Av explains the association between dense environments and high Mline-core values, and argues against filaments remaining in constant single-pressure equilibrium states. The simultaneous growth of filament and background in locations with efficient mass assembly, predicted in numerical models of collapsing clouds, presents a suitable scenario for the fulfillment of the combined filament mass-environment criterium that is in quantitative agreement with Herschel observations. △ Less

Submitted 24 February, 2017; originally announced February 2017.

Comments: Accepted in A&A

Journal ref: A&A 601, A94 (2017)

Galactic Cold Cores VII: Filament Formation and Evolution - Methods & Observational Constraints

Authors: A. Rivera-Ingraham, I. Ristorcelli, M. Juvela, J. Montillaud, A. Men'shchikov, J. Malinen, V. -M. Pelkonen, A. Marston, P. G. Martin, L. Pagani, R. Paladini, D. Paradis, N. Ysard, D. Ward-Thompson, J. -P. Bernard, D. J. Marshall, L. Montier, V. Tóth

Abstract: The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role… ▽ More The association of filaments with protostellar objects has made these structures a priority target in star formation studies. The datasets of the Herschel Galactic Cold Cores Key Programme allow for a statistical study of filaments with a wide range of intrinsic and environmental characteristics. Characterisation of this sample can be used to identify key physical parameters and quantify the role of environment in the formation of supercritical filaments. Filaments were extracted from fields at D<500pc with the getfilaments algorithm and characterised according to their column density profiles and intrinsic properties. Each profile was fitted with a beam-convolved Plummer-like function and quantified based on the relative contributions from the filament 'core', represented by a Gaussian, and 'wing' component, dominated by the power-law of the Plummer-like function. These parameters were examined for populations associated with different background levels. We find that filaments increase their core (Mcore) and wing (Mwing) contributions while increasing their total linear mass density (Mtot). Both components appear to be linked to the local environment, with filaments in higher backgrounds having systematically more massive Mcore and Mwing. This dependence on the environment supports an accretion-based model for filament evolution in the local neighbourhood (D<500pc). Structures located in the highest backgrounds develop the highest central Av, Mcore, and Mwing as Mtot increases with time, favoured by the local availability of material and the enhanced gravitational potential. Our results indicate that filaments acquiring a significantly massive central region with Mcore>Mcrit/2 may become supercritical and form stars. This translates into a need for filaments to become at least moderately self-gravitating in order to undergo localised star formation or become star-forming filaments. △ Less

Submitted 30 March, 2016; originally announced March 2016.

Comments: Accepted for publication in A&A

Journal ref: A&A 591, A90 (2016)

Spectral Study of the Vocal Tract in Vowel Synthesis: A Comparison between 1D and 3D Acoustic Analysis

Authors: Negar M. Harandi, Daniel Aalto, Antti Hannukainen, Jarmo Malinen, Sidney Fels

Abstract: A state-of-the-art 1D acoustic synthesizer has been previously developed, and coupled to speaker-specific biomechanical models of oropharynx in ArtiSynth. As expected, the formant frequencies of the synthesized vowel sounds were shown to be different from those of the recorded audio. Such discrepancy was hypothesized to be due to the simplified geometry of the vocal tract model as well as the one… ▽ More A state-of-the-art 1D acoustic synthesizer has been previously developed, and coupled to speaker-specific biomechanical models of oropharynx in ArtiSynth. As expected, the formant frequencies of the synthesized vowel sounds were shown to be different from those of the recorded audio. Such discrepancy was hypothesized to be due to the simplified geometry of the vocal tract model as well as the one dimensional implementation of Navier-Stokes equations. In this paper, we calculate Helmholtz resonances of our vocal tract geometries using 3D finite element method (FEM), and compare them with the formant frequencies obtained from the 1D method and audio. We hope such comparison helps with clarifying the limitations of our current models and/or speech synthesizer. △ Less

Submitted 17 December, 2015; originally announced December 2015.

Matching dust emission structures and magnetic field in high-latitude cloud L1642: comparing Herschel and Planck maps

Authors: J. Malinen, L. Montier, J. Montillaud, M. Juvela, I. Ristorcelli, S. E. Clark, O. Berné, J. -Ph. Bernard, V. -M. Pelkonen, D. C. Collins

Abstract: The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642… ▽ More The nearby cloud L1642 is one of only two known very high latitude (|b| > 30 deg) clouds actively forming stars. It is a rare example of star formation in isolated conditions, and can reveal important details of star formation in general, e.g., of the effect of magnetic fields. We compare Herschel dust emission structures and magnetic field orientation revealed by Planck polarization maps in L1642. The high-resolution ($\sim20"$) Herschel data reveal a complex structure including a dense, compressed central clump, and low density striations. The Planck polarization data (at 10$'$ resolution) reveal an ordered magnetic field pervading the cloud and aligned with the surrounding striations. There is a complex interplay between the cloud structure and large scale magnetic field. This suggests that the magnetic field is closely linked to the formation and evolution of the cloud. CO rotational emission confirms that the striations are connected with the main clumps and likely to contain material either falling into or flowing out of the clumps. There is a clear transition from aligned to perpendicular structures approximately at a column density of $N_{\rm{H}} = 1.6 \times 10^{21}\, {\rm cm}^{-2}$. Comparing the Herschel maps with the Planck polarization maps shows the close connection between the magnetic field and cloud structure even in the finest details of the cloud. △ Less

Submitted 3 May, 2016; v1 submitted 11 December, 2015; originally announced December 2015.

Comments: accepted to MNRAS; 14 pages; changes: moderate text and figure editing

Galactic cold cores VI. Dust opacity spectral index

Authors: M. Juvela, K. Demyk, Y. Doi, A. Hughes, C. Lefevre, D. J. Marshall, C. Meny, J. Montillaud, L. Pagani, D. Paradis, I. Ristorcelli, J. Malinen, L. A. Montier, R. Paladini, V. -M. Pelkonen, A. Rivera-Ingraham

Abstract: The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the… ▽ More The Galactic Cold Cores project has made Herschel observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. The dust opacity spectral index beta and the dust colour temperature T are derived using Herschel and Planck data. The relation between beta and T is examined for the whole sample and inside individual fields. Based on IRAS and Planck data, the fields are characterised by a median colour temperature of 16.1 K and a median opacity spectral index of beta=1.84. We observe a clear T-beta anti-correlation. In Herschel observations, constrained at lower resolution by Planck data, the variations follow the column density structure and beta(FIR) can rise to ~2.2 in individual clumps. The Planck 217 GHz band shows a systematic excess that is consistent with a general flattening of the dust emission spectrum at millimetre wavelengths. When fitted separately below and above 700 um, the median spectral index values are beta(FIR) ~ 1.91 and beta(mm) ~ 1.66. The spectral index changes as a function of column density and wavelength. Beta variations are partly masked by temperature gradients and the changes in the intrinsic grain properties may be even greater. △ Less

Submitted 26 September, 2015; originally announced September 2015.

Comments: Accepted to A&A, 138 pages

Journal ref: A&A 584, A94 (2015)

Post-processing speech recordings during MRI

Authors: Juha Kuortti, Jarmo Malinen, Antti Ojalammi

Abstract: We discuss post-processing of speech that has been recorded during Magnetic Resonance Imaging (MRI) of the vocal tract. Such speech recordings are contaminated by high levels of acoustic noise from the MRI scanner. Also, the frequency response of the sound signal path is not flat as a result of severe restrictions on recording instrumentation due to MRI technology. The post-processing algorithm… ▽ More We discuss post-processing of speech that has been recorded during Magnetic Resonance Imaging (MRI) of the vocal tract. Such speech recordings are contaminated by high levels of acoustic noise from the MRI scanner. Also, the frequency response of the sound signal path is not flat as a result of severe restrictions on recording instrumentation due to MRI technology. The post-processing algorithm for noise reduction is based on adaptive spectral filtering. The speech material consists of samples of prolonged vowel productions that are used for validation of the post-processing algorithm. The comparison data is recorded in anechoic chamber from the same test subject. Formant analysis is carried out for the post-processed speech and the comparison data. Artificially noise-contaminated vowel samples are used for validation experiments to determine performance of the algorithm where using true data would be difficult. The properties of recording instrumentation or the post-processing algorithm do not explain the consistent frequency dependent discrepancy between formant data from experiments during MRI and in anechoic chamber. It is shown that the discrepancy is statistically significant, in particular, where it is largest at 1 kHzきろへるつ and 2 kHzきろへるつ. The reflecting surfaces of the MRI head and neck coil are suspected to change the speech acoustics which results in "external formants" at these frequencies. However, the role of test subject adaptation to noise and constrained space acoustics during an MRI examination cannot be ruled out. △ Less

Submitted 21 June, 2016; v1 submitted 17 September, 2015; originally announced September 2015.

Modal locking between vocal fold and vocal tract oscillations: Simulations in time domain

Authors: Atte Aalto, Tiina Murtola, Jarmo Malinen, Daniel Aalto, Martti Vainio

Abstract: During voiced speech, the human vocal folds interact with the vocal tract acoustics. The resulting glottal source-resonator coupling has been observed using mathematical and physical models as well as in in vivo phonation. We propose a computational time-domain model of the full speech apparatus that, in particular, contains a feedback mechanism from the vocal tract acoustics to the vocal fold osc… ▽ More During voiced speech, the human vocal folds interact with the vocal tract acoustics. The resulting glottal source-resonator coupling has been observed using mathematical and physical models as well as in in vivo phonation. We propose a computational time-domain model of the full speech apparatus that, in particular, contains a feedback mechanism from the vocal tract acoustics to the vocal fold oscillations. It is based on numerical solution of ordinary and partial differential equations defined on vocal tract geometries that have been obtained by Magnetic Resonance Imaging. The model is used to simulate rising and falling pitch glides of [a, i] in the fundamental frequency (f_o) interval [150 Hzへるつ, 320 Hzへるつ]. The interval contains the first vocal tract resonance f_R1 and the first formant F1 of [i] as well as the fractions of the first resonance f_R1/4 and fR1/3 of [a]. The simulations reveal a locking pattern of the fo-trajectory at f_R1 of [i] in falling and rising glides. The resonance fractions of [a] produce perturbations in the pressure signal at the lips but no locking. All these observations from the model behaviour are consistent and robust within a wide range of feasible model parameter values and under exclusion of secondary model components. △ Less

Submitted 15 March, 2017; v1 submitted 3 June, 2015; originally announced June 2015.

MSC Class: 76Q05; 92C05; 93C20

The physical state of selected cold clumps

Authors: A. Parikka, M. Juvela, V. -M. Pelkonen, J. Malinen, J. Harju

Abstract: Context. he study of prestellar cores is essential to understanding the initial stages of star formation. With $Herschel$ more cold clumps have been detected than ever before. For this study we have selected 21 cold clumps from 20 $Herschel$ fields observed as a follow-up on original $Planck$ detections. We have observed these clumps in $^{13}$CO (1-0), C$^{18}$O (1-0), and N$_2$H$^+$ (1-0) lines.… ▽ More Context. he study of prestellar cores is essential to understanding the initial stages of star formation. With $Herschel$ more cold clumps have been detected than ever before. For this study we have selected 21 cold clumps from 20 $Herschel$ fields observed as a follow-up on original $Planck$ detections. We have observed these clumps in $^{13}$CO (1-0), C$^{18}$O (1-0), and N$_2$H$^+$ (1-0) lines. Aims. Our aim is to find out if these cold clumps are prestellar. We have examined to what extent independent analysis of the dust and the molecular lines lead to similar conclusions about the masses of these objects. Methods. We calculate the clump masses and densities from the dust continuum and molecular line observations and compare these to each other and to the virial and Bonnor-Ebert masses calculated for each clump. Finally we examine two of the fields with radiative transfer models to estimate CO abundances. Results. When excitation temperatures could be estimated, the column densities derived from molecular line observations were comparable to those from dust continuum data. The median column density estimates are 4.2$\times 10^{21}$cm$^{-2}$ and 5.5$\times 10^{21}$cm$^{-2}$ for the line and dust emission data, respectively. The calculated abundances, column densities, volume densities, and masses all have large uncertainties and one must be careful when drawing conclusions. Abundance of $^{13}$CO was found in modeling the two clumps in the field G131.65$+$9.75 to be close to the usual value of 10$^{-6}$. The abundance ratio of $^{13}$CO and C$^{18}$O was $\sim$10. Molecular abundances could only be estimated with modeling, relying on dust column density data. Conclusions. The results indicate that most cold clumps, even those with dust color temperatures close to 11 K, are not necessarily prestellar. △ Less

Submitted 20 March, 2015; originally announced March 2015.

Comments: Accepted to A&A

Journal ref: A&A 577, A69 (2015)

Galactic cold cores V. Dust opacity

Authors: M. Juvela, I. Ristorcelli, D. J. Marshall, J. Montillaud, V. -M. Pelkonen, N. Ysard, P. McGehee, R. Paladini, L. Pagani, J. Malinen, A., Rivera-Ingraham, C. Lefevre, L. V. Toth, L. A. Montier, J. -P. Bernard, P. Martin

Abstract: The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observ… ▽ More The project Galactic Cold Cores has made Herschel observations of interstellar clouds where the Planck satellite survey has located cold and compact clumps. The sources range from starless clumps to protostellar cores. We examine 116 Herschel fields to estimate the submillimetre dust opacity and its variations. The submillimetre dust opacity was derived from Herschel data, and near-infrared observations of the reddening of background stars are converted into near-infrared optical depth. We studied the systematic errors affecting these parameters and used modelling to correct for the expected biases. The ratio of 250um and J band opacities is correlated with the cloud location and star formation activity. We find a median ratio of tau(250um)/tau(J)= (1.6+-0.2)*10^-3, which is more than three times the mean value in diffuse medium. Assuming a spectral index beta=1.8 instead of beta=2.0, the value would be lower by ~30%. No significant systematic variation is detected with Galactocentric distance or with Galactic height. The tau(250um)/tau(J) maps reveal six fields with clear increase of submillimetre opacity of up to tau(250um)/tau(J) ~ 4*10^-3. These are all nearby fields with spatially resolved clumps of high column density. We interpret the increase in the far-infrared opacity as a sign of grain growth in the densest and coldest regions of interstellar clouds. △ Less

Submitted 28 January, 2015; originally announced January 2015.

Comments: 28 pages, A&A, in press

A posteriori error estimates for Webster's equation in wave propagation

Authors: Teemu Lukkari, Jarmo Malinen

Abstract: We consider a generalised Webster's equation for describing wave propagation in curved tubular structures such as variable diameter acoustic wave guides. Webster's equation in generalised form has been rigorously derived in a previous article starting from the wave equation, and it approximates cross-sectional averages of the propagating wave. Here, the approximation error is estimated by an a pos… ▽ More We consider a generalised Webster's equation for describing wave propagation in curved tubular structures such as variable diameter acoustic wave guides. Webster's equation in generalised form has been rigorously derived in a previous article starting from the wave equation, and it approximates cross-sectional averages of the propagating wave. Here, the approximation error is estimated by an a posteriori technique. △ Less

Submitted 28 May, 2014; originally announced May 2014.

MSC Class: Primary 37L05. Secondary 35L05; 35L20; 47N70; 93C20

Mapping of interstellar clouds with infrared light scattered from dust: TMC-1N

Authors: J. Malinen, M. Juvela, V. -M. Pelkonen, M. G. Rawlings

Abstract: Mapping of near-infrared (NIR) scattered light is a recent method for the study of interstellar clouds, complementing other, more commonly used methods, like dust emission and extinction. Our goal is to study the usability of this method on larger scale, and compare the properties of a filament using NIR scattering and other methods. We also study the radiation field and differences in grain emiss… ▽ More Mapping of near-infrared (NIR) scattered light is a recent method for the study of interstellar clouds, complementing other, more commonly used methods, like dust emission and extinction. Our goal is to study the usability of this method on larger scale, and compare the properties of a filament using NIR scattering and other methods. We also study the radiation field and differences in grain emissivity between diffuse and dense areas. We have used scattered J, H, and K band surface brightness WFCAM-observations to map filament TMC-1N in Taurus, covering an area of 1dx1d corresponding to ~(2.44 pc)^2. We have converted the data into optical depth and compared the results with NIR extinction and Herschel observations of submm dust emission. We see the filament in scattered light in all three NIR bands. We note that our WFCAM observations in TMC-1N show notably lower intensity than previous results in Corona Australis using the same method. We show that 3D radiative transfer simulations predict similar scattered surface brightness levels as seen in the observations. However, changing the assumptions about the background can change the results of simulations notably. We derive emissivity by using optical depth in the J band as an independent tracer of column density. We obtain opacity sigma(250um) values 1.7-2.4x10^-25 cm^2/H, depending on assumptions of the extinction curve, which can change the results by over 40%. These values are twice as high as obtained for diffuse areas, at the lower limit of earlier results for denser areas. We show that NIR scattering can be a valuable tool in making high resolution maps. We conclude, however, that NIR scattering observations can be complicated, as the data can show relatively low-level artefacts. This suggests caution when planning and interpreting the observations. △ Less

Submitted 9 April, 2014; originally announced April 2014.

Comments: abstract shortened and figures reduced for astroph

Journal ref: A&A 558, A44 (2013)

Multiwavelength study of the high-latitude cloud L1642: chain of star formation

Authors: J. Malinen, M. Juvela, S. Zahorecz, A. Rivera-Ingraham, J. Montillaud, K. Arimatsu, J. -Ph. Bernard, Y. Doi, L. Haikala, R. Kawabe, G. Marton, P. McGehee, V. -M. Pelkonen, I. Ristorcelli, Y. Shimajiri, S. Takita, L. V. Toth, T. Tsukagoshi, N. Ysard

Abstract: L1642 is one of the two high galactic latitude (|b| > 30deg) clouds confirmed to have active star formation. We examine the properties of this cloud, especially the large-scale structure, dust properties, and compact sources in different stages of star formation. We present high-resolution far-infrared and submm observations with the Herschel and AKARI satellites and mm observations with the AzTEC… ▽ More L1642 is one of the two high galactic latitude (|b| > 30deg) clouds confirmed to have active star formation. We examine the properties of this cloud, especially the large-scale structure, dust properties, and compact sources in different stages of star formation. We present high-resolution far-infrared and submm observations with the Herschel and AKARI satellites and mm observations with the AzTEC/ASTE telescope, which we combined with archive data from near- and mid-infrared (2MASS, WISE) to mm observations (Planck). The Herschel observations, combined with other data, show a sequence of objects from a cold clump to young stellar objects at different evolutionary stages. Source B-3 (2MASS J04351455-1414468) appears to be a YSO forming inside the L1642 cloud, instead of a foreground brown dwarf, as previously classified. Herschel data reveal striation in the diffuse dust emission around L1642. The western region shows striation towards NE and has a steeper column density gradient on its southern side. The densest central region has a bow-shock like structure showing compression from the west and a filamentary tail extending towards east. The differences suggest that these may be spatially distinct structures, aligned only in projection. We derive values of the dust emission cross-section per H nucleon for different regions of the cloud. Modified black-body fits to the spectral energy distribution of Herschel and Planck data give emissivity spectral index beta values 1.8-2.0 for the different regions. The compact sources have lower beta values and show an anticorrelation between T and beta. Markov chain Monte Carlo calculations demonstrate the strong anticorrelation between beta and T errors and the importance of mm Planck data in constraining the estimates. L1642 reveals a more complex structure and sequence of star formation than previously known. △ Less

Submitted 9 April, 2014; v1 submitted 11 February, 2014; originally announced February 2014.

Comments: 22 pages, 18 figures, accepted to Astronomy & Astrophysics; abstract shortened and figures reduced for astroph

Measurement of acoustic and anatomic changes in oral and maxillofacial surgery patients

Authors: Daniel Aalto, Olli Aaltonen, Risto-Pekka Happonen, Päivi Jääsaari, Atle Kivelä, Juha Kuortti, Jean-Marc Luukinen, Jarmo Malinen, Tiina Murtola, Riitta Parkkola, Jani Saunavaara, Tero Soukka, Martti Vainio

Abstract: We describe an arrangement for simultaneous recording of speech and geometry of vocal tract in patients undergoing surgery involving this area. Experimental design is considered from an articulatory phonetic point of view. The speech and noise signals are recorded with an acoustic-electrical arrangement. The vocal tract is simultaneously imaged with MRI. A MATLAB-based system controls the timing o… ▽ More We describe an arrangement for simultaneous recording of speech and geometry of vocal tract in patients undergoing surgery involving this area. Experimental design is considered from an articulatory phonetic point of view. The speech and noise signals are recorded with an acoustic-electrical arrangement. The vocal tract is simultaneously imaged with MRI. A MATLAB-based system controls the timing of speech recording and MR image acquisition. The speech signals are cleaned from acoustic MRI noise by a non-linear signal processing algorithm. Finally, a vowel data set from pilot experiments is compared with validation data from anechoic chamber as well as with Helmholtz resonances of the vocal tract volume. △ Less

Submitted 11 September, 2013; originally announced September 2013.

Comments: 19 pages, 17 figures

MSC Class: 92C50; 92C55

Estimation of high-resolution dust column density maps. Comparison of modified black-body fits and radiative transfer modelling

Authors: M. Juvela, J. Malinen, T. Lunttila

Abstract: Sub-millimetre dust emission is often used to derive the column density N of dense interstellar clouds. The observations consist of data at several wavelengths but of variable resolution. We examine two procedures that been proposed for the estimation of high resolution N maps. Method A uses a low-resolution temperature map combined with higher resolution intensity data while Method B combines N e… ▽ More Sub-millimetre dust emission is often used to derive the column density N of dense interstellar clouds. The observations consist of data at several wavelengths but of variable resolution. We examine two procedures that been proposed for the estimation of high resolution N maps. Method A uses a low-resolution temperature map combined with higher resolution intensity data while Method B combines N estimates from different wavelength ranges. Our aim is to determine the accuracy of the methods relative to the true column densities and the estimates obtainable with radiative transfer modelling. We use magnetohydrodynamical (MHD) simulations and radiative transfer calculations to simulate sub-millimetre observations at the wavelengths of the Herschel Space Observatory. The observations are analysed with the methods and the results compared to the true values and to the results from radiative transfer modelling of observations. Both methods A and B give relatively reliable column density estimates at the resolution of 250um data while also making use of the longer wavelengths. For high signal-to-noise data, the results of Method B are better correlated with the true column density, while Method A is less sensitive to noise. When the cloud has internal heating, results of Method B are consistent with those that would be obtained with high-resolution data. Because of line-of-sight temperature variations, these underestimate the true column density and, because of a favourable cancellation of errors, Method A can sometimes give more correct values. Radiative transfer modelling, even with very simple 3D cloud models, can provide better results. However, the complexity of the models required for improvements increases rapidly with the complexity and opacity of the clouds. △ Less

Submitted 9 May, 2013; originally announced May 2013.

Comments: 14 pages, Accepted to A&A

Acoustic wave guides as infinite-dimensional dynamical systems

Authors: Atte Aalto, Teemu Lukkari, Jarmo Malinen

Abstract: We prove the unique solvability, passivity/conservativity and some regularity results of two mathematical models for acoustic wave propagation in curved, variable diameter tubular structures of finite length. The first of the models is the generalised Webster's model that includes dissipation and curvature of the 1D waveguide. The second model is the scattering passive, boundary controlled wave eq… ▽ More We prove the unique solvability, passivity/conservativity and some regularity results of two mathematical models for acoustic wave propagation in curved, variable diameter tubular structures of finite length. The first of the models is the generalised Webster's model that includes dissipation and curvature of the 1D waveguide. The second model is the scattering passive, boundary controlled wave equation on 3D waveguides. The two models are treated in an unified fashion so that the results on the wave equation reduce to the corresponding results of approximating Webster's model at the limit of vanishing waveguide intersection. △ Less

Submitted 30 May, 2013; v1 submitted 29 November, 2012; originally announced November 2012.

Comments: 29 pp with Appendix 7pp

MSC Class: Primary 35L05; secondary 35L20; 93C20; 47N70

Microspectral analysis of quasinilpotent operators

Authors: Jarmo Malinen, Olavi Nevanlinna, Jaroslav Zemánek

Abstract: We develop a microspectral theory for quasinilpotent linear operators $Q$ (i.e., those with $σしぐま(Q) = \{0}$) in a Banach space. When such $Q$ is not compact, normal, or nilpotent, the classical spectral theory gives little information, and a somewhat deeper structure can be recovered from microspectral sets in $\C$. Such sets describe, e.g., semigroup generation, resolvent properties, power boundedn… ▽ More We develop a microspectral theory for quasinilpotent linear operators $Q$ (i.e., those with $σしぐま(Q) = \{0}$) in a Banach space. When such $Q$ is not compact, normal, or nilpotent, the classical spectral theory gives little information, and a somewhat deeper structure can be recovered from microspectral sets in $\C$. Such sets describe, e.g., semigroup generation, resolvent properties, power boundedness as well as Tauberian properties associated to $zQ$ for $z \in \C$. △ Less

Submitted 20 November, 2012; originally announced November 2012.

MSC Class: 47A10; 47B06; 47B10; 47B60

Modal locking between vocal fold and vocal tract oscillations: Experiments and statistical analysis

Authors: Daniel Aalto, Jarmo Malinen, Martti Vainio

Abstract: The human vocal folds are known to interact with the vocal tract acoustics during voiced speech production; namely a nonlinear source-filter coupling has been observed both by using models and in \emph{in vivo} phonation. These phenomena are approached from two directions in this article. We first present a computational dynamical model of the speech apparatus that contains an explicit filter-sour… ▽ More The human vocal folds are known to interact with the vocal tract acoustics during voiced speech production; namely a nonlinear source-filter coupling has been observed both by using models and in \emph{in vivo} phonation. These phenomena are approached from two directions in this article. We first present a computational dynamical model of the speech apparatus that contains an explicit filter-source feedback mechanism from the vocal tract acoustics back to the vocal folds oscillations. The model was used to simulate vocal pitch glideswhere the trajectory was forced to cross the lowest vocal tract resonance, i.e., the lowest formant $F_1$. Similar patterns produced by human participants were then studied. Both the simulations and the experimental results reveal an effect when the glides cross the first formant (as may happen in \textipa{[i]}). Conversely, this effect is not observed if there is no formant within the glide range (as is the case in \textipa{[\textscripta]}). The experiments show smaller effect compared to the simulations, pointing to an active compensation mechanism. △ Less

Submitted 16 November, 2015; v1 submitted 20 November, 2012; originally announced November 2012.

Comments: PACS: Primary 4370Bk. Secondary 4370Gr, 4370Aj, 4370Mn

How far are vowel formants from computed vocal tract resonances?

Authors: Daniel Aalto, Antti Huhtala, Atle Kivelä, Jarmo Malinen, Pertti Palo, Jani Saunavaara, Martti Vainio

Abstract: We compare numerically computed resonances of the human vocal tract with formants that have been extracted from speech during vowel pronunciation. The geometry of the vocal tract has been obtained by MRI from a male subject, and the corresponding speech has been recorded simultaneously. The resonances are computed by solving the Helmholtz partial differential equation with the Finite Element Metho… ▽ More We compare numerically computed resonances of the human vocal tract with formants that have been extracted from speech during vowel pronunciation. The geometry of the vocal tract has been obtained by MRI from a male subject, and the corresponding speech has been recorded simultaneously. The resonances are computed by solving the Helmholtz partial differential equation with the Finite Element Method (FEM). Despite a rudimentary exterior space acoustics model, i.e., the Dirichlet boundary condition at the mouth opening, the computed resonance structure differs from the measured formant structure by $\approx$ 0.7 semitones for [i] and [u] having small mouth opening area, and by $\approx$ 3 semitones for vowels [a] and [ae] that have a larger mouth opening. The contribution of the possibly open velar port has not been taken into considaration at all which adds the discrepancy for [a] in the present data set. We conclude that by improving the exterior space model and properly treating the velar port opening, it is possible to computationally attain four lowest vowel formants with an error less than a semitone. The corresponding wave equation model on MRI-produced vocal tract geometries is expected to have a comparable accuracy. △ Less

Submitted 8 October, 2012; v1 submitted 29 August, 2012; originally announced August 2012.

Comments: 13 pages, 1 figure, 3 tables

MSC Class: 37N25; 92C10

Profiles of interstellar cloud filaments. Observational effects in synthetic sub-millimetre observations

Authors: M. Juvela, J. Malinen, T. Lunttila

Abstract: Sub-millimetre observations suggest that the filaments of interstellar clouds have rather uniform widths and can be described with the so-called Plummer profiles. The shapes of the filament profiles are linked to their physical state. Before drawing conclusions on the observed column density profiles, we must evaluate the observational uncertainties. We want to estimate the bias that could result… ▽ More Sub-millimetre observations suggest that the filaments of interstellar clouds have rather uniform widths and can be described with the so-called Plummer profiles. The shapes of the filament profiles are linked to their physical state. Before drawing conclusions on the observed column density profiles, we must evaluate the observational uncertainties. We want to estimate the bias that could result from radiative transfer effects or from variations of submm dust emissivity. We use cloud models obtained with magnetohydrodynamic simulations and carry out radiative transfer calculations to produce maps of sub-millimetre emission. Column densities are estimated based on the synthetic observations. For selected filaments, the estimated profiles are compared to those derived from the original column density. Possible effects from spatial variations of dust properties are examined. With instrumental noise typical of the Herschel observations, the parameters derived for nearby clouds are correct to within a few percent. The radiative transfer effects have only a minor effect on the results. If the signal-to-noise ratio is degraded by a factor of four, the errors become significant and for half of the examined filaments the values cannot be constrained. The errors increase in proportion to the cloud distance. Assuming the resolution of Herschel instruments, the model filaments are barely resolved at a distance of ~400 pc and the errors in the parameters of the Plummer function are several tens of per cent. The Plummer parameters, in particular the power-law exponent p, are sensitive to noise but can be determined with good accuracy using Herschel data. One must be cautious about possible line-of-sight confusion. In our models, a large fraction of the filaments seen in the column density maps are not continuous structures in three dimensions. △ Less

Submitted 16 July, 2012; originally announced July 2012.

Comments: 12 pages, 14 figures, accepted to A&A

Profiling filaments: comparing near-infrared extinction and submillimetre data in TMC-1

Authors: J. Malinen, M. Juvela, M. G. Rawlings, D. Ward-Thompson, P. Palmeirim, Ph. Andre

Abstract: Interstellar filaments are an important part of star formation. To understand the structure of filaments, cross-section profiles are often fitted with Plummer profiles. This profiling is often done with submm studies, such as Herschel. It would be convenient if filament properties could also be studied using groundbased NIR data. We compare the filament profiles obtained by NIR extinction and subm… ▽ More Interstellar filaments are an important part of star formation. To understand the structure of filaments, cross-section profiles are often fitted with Plummer profiles. This profiling is often done with submm studies, such as Herschel. It would be convenient if filament properties could also be studied using groundbased NIR data. We compare the filament profiles obtained by NIR extinction and submm observations to find out if reliable profiles can be derived using NIR data. We use J-, H-, and K-band data of a filament north of TMC-1 to derive an extinction map from colour excesses of background stars. We compare the Plummer profiles obtained from extinction maps with Herschel dust emission maps. We present 2 methods to estimate profiles from NIR: Plummer profile fits to median Av of stars or directly to the Av of individual stars. We compare the methods by simulations. In simulations extinction maps and the new methods give correct results to within ~10-20 for modest densities. Direct fit to data on individual stars gives more accurate results than extinction map, and can work in higher density. In profile fits to real observations, values of Plummer parameters are generally similar to within a factor of ~2. Although parameter values can vary significantly, estimates of filament mass usually remain accurate to within some tens of per cent. Our results for TMC-1 are in agreement with earlier results. High resolution NIR data give more details, but 2MASS data can be used to estimate profiles. NIR extinction can be used as an alternative to submm observations to profile filaments. Direct fits of stars can also be a valuable tool. Plummer profile parameters are not always well constrained, and caution should be taken when making fits. In the evaluation of Plummer parameters, one can use the independence of dust emission and NIR data and the difference in the shapes of the confidence regions. △ Less

Submitted 21 June, 2012; originally announced June 2012.

Comments: accepted to Astronomy & Astrophysics; abstract has been shortened for astroph

Webster's equation with curvature and dissipation

Authors: Teemu Lukkari, Jarmo Malinen

Abstract: Wave propagation in curved tubular domains is considered. A general version of Webster's equation is derived from the scattering passive wave equation. More precisely, it is shown that planar averages of a sufficiently smooth solution of the wave equation satisfy the corresponding Webster's equation when the latter includes additional control signals determined by the solution. Wave propagation in curved tubular domains is considered. A general version of Webster's equation is derived from the scattering passive wave equation. More precisely, it is shown that planar averages of a sufficiently smooth solution of the wave equation satisfy the corresponding Webster's equation when the latter includes additional control signals determined by the solution. △ Less

Submitted 28 June, 2013; v1 submitted 18 April, 2012; originally announced April 2012.

Comments: 27pp., 3pp. appendix, 3 figures

MSC Class: 35L05; 35L20; 93C20; 47N70

Modelling the dust emission from dense interstellar clouds: disentangling the effects of radiative transfer and dust properties

Authors: N. Ysard, M. Juvela, K. Demyk, V. Guillet, A. Abergel, J. -P. Bernard, J. Malinen, C. Mény, L. Montier, D. Paradis, I. Ristorcelli, L. Verstraete

Abstract: With Planck and Herschel, we now have the spectral coverage and angular resolution required to observe dense and cold molecular clouds. As these clouds are optically thick at short wavelength but optically thin at long wavelength, it is tricky to conclude anything about dust properties without a proper treatment of the radiative transfer (RT). Our aim is to disentangle the effects of RT and of dus… ▽ More With Planck and Herschel, we now have the spectral coverage and angular resolution required to observe dense and cold molecular clouds. As these clouds are optically thick at short wavelength but optically thin at long wavelength, it is tricky to conclude anything about dust properties without a proper treatment of the radiative transfer (RT). Our aim is to disentangle the effects of RT and of dust properties on the variations in the dust emission to provide observers with keys to analyse the emission arising from dense clouds. We model cylindrical clouds, illuminated by the ISRF, and carry out full RT calculations. Dust temperatures are solved using DustEM for amorphous carbons and silicates, grains coated with carbon mantles, and mixed aggregates of carbon and silicate. We allow variations of the grain optical properties with wavelength and temperature. We determine observed colour temperatures, T, and emissivity spectral indices, beta, by fitting the dust emission with modified blackbodies, to compare our models with observations. RT effects can neither explain the low T nor the increased submm emissivity measured at the centre of dense clouds, nor the observed beta-T anti-correlation. Adding noise to the modelled data, we show that it is not likely to be the unique explanation for the beta-T anti-correlation observed in starless clouds. It may be explained by intrinsic variations in the grain optical properties with temperature. As for the increased submm emissivity and the low T, they have to originate in variations in the grain optical properties, probably caused by their growth to form porous aggregates. We find it difficult to track back the nature of grains from the spectral variations in their emission. Finally, the column density is underestimated when determined with blackbody fitting because of the discrepancy between T and the true dust temperature at the cloud centre. △ Less

Submitted 25 April, 2012; v1 submitted 27 February, 2012; originally announced February 2012.

Comments: accepted for publication by A&A, abstract modified to match arXiv format, version 2 => correction of numerical values throughout the text/figures

Galactic cold cores III. General cloud properties

Authors: M. Juvela, I. Ristorcelli, L. Pagani, Y. Doi, V. -M. Pelkonen, D. J. Marshall, J. -P. Bernard, E. Falgarone, J. Malinen, G. Marton, P. McGehee, L. A. Montier, F. Motte, R. Paladini, L. V. Toth, N. Ysard, S. Zahorecz, A. Zavagno

Abstract: We examine the cloud structure around the Planck detections in 71 fields observed with the Herschel SPIRE instrument. We wish to determine the general physical characteristics of the fields and to examine the morphology of the clouds where the cold high column density clumps are found. We derive colour temperature and column density maps of the fields. We examine the infrared spectral energy distr… ▽ More We examine the cloud structure around the Planck detections in 71 fields observed with the Herschel SPIRE instrument. We wish to determine the general physical characteristics of the fields and to examine the morphology of the clouds where the cold high column density clumps are found. We derive colour temperature and column density maps of the fields. We examine the infrared spectral energy distributions of the main clumps. The clouds are categorised according to their large scale morphology. With the help of recently released WISE satellite data, we look for signs of enhanced mid-infrared scattering (coreshine), an indication of growth of the dust grains, and examine the star formation activity associated with the cold clumps. The mapped clouds have distances ranging from ~100pc to several kiloparsecs and cover a range of sizes and masses from cores of less than 10 solar masses to clouds with masses in excess of 10000 solar mass. Most fields contain some filamentary structures and in about half of the cases a filament or a few filaments dominate the morphology. In one case out of ten, the clouds show a cometary shape or have sharp boundaries indicative of compression by an external force. The width of the filaments is typically ~0.2-0.3pc. However, there is significant variation from 0.1pc to 1pc and the estimates are sensitive to the methods used and the very definition of a filament. Enhanced mid-infrared scattering, coreshine, was detected in four clouds with six additional tentative detections. The cloud LDN183 is included in our sample and remains the best example of this phenomenon. About half of the fields are associated with active star formation as indicated by the presence of mid-infrared point sources. The mid-infrared sources often coincide with structures whose sub-millimetre spectra are still dominated by the cold dust. △ Less

Submitted 8 February, 2012; originally announced February 2012.

Comments: Accepted to A&A. The pdf with high resolution figures and with the figures of the appendices is available at https://wiki.helsinki.fi/display/PlanckHerschel/Publications

Galactic cold cores II. Herschel study of the extended dust emission around the first Planck detections

Authors: M. Juvela, I. Ristorcelli, V. -M. Pelkonen, D. J. Marshall, L. A. Montier, J. -P. Bernard, R. Paladini, T. Lunttila, A. Abergel, Ph. Andre, C. Dickinson, X. Dupac, J. Malinen, P. Martin, P. McGehee, L. Pagani, N. Ysard, A. Zavagno

Abstract: Within the project Galactic Cold Cores we are carrying out Herschel photometric observations of cold interstellar clouds detected with the Planck satellite. The three fields observed as part of the Herschel science demonstration phase (SDP) provided the first glimpse into the nature of these sources. We examine the properties of the dust emission within the fields. We determine the dust sub-millim… ▽ More Within the project Galactic Cold Cores we are carrying out Herschel photometric observations of cold interstellar clouds detected with the Planck satellite. The three fields observed as part of the Herschel science demonstration phase (SDP) provided the first glimpse into the nature of these sources. We examine the properties of the dust emission within the fields. We determine the dust sub-millimetre opacity, look for signs of spatial variations in the dust spectral index, and estimate how the apparent variations of the parameters could be affected by different sources of uncertainty. We use the Herschel observations where the zero point of the surface brightness scale is set with the help of the Planck satellite data. We derive the colour temperature and column density maps of the regions and determine the dust opacity by a comparison with extinction measurements. By simultaneously fitting the colour temperature and the dust spectral index values we look for spatial variations in the apparent dust properties. With a simple radiative transfer model we estimate to what extent these can be explained by line-of-sight temperature variations, without changes in the dust grain properties. The analysis of the dust emission reveals cold and dense clouds that coincide with the Planck sources and confirm those detections. The derived dust opacity varies in the range kappa(250um) ~ 0.05-0.2 cm^2/g, higher values being observed preferentially in regions of high column density. The average dust spectral index beta is ~ 1.9-2.2. There are indications that beta increases towards the coldest regions. The spectral index decreases strongly near internal heating sources but, according to radiative transfer models, this can be explained by the line-of-sight temperature variations without a change in the dust properties. △ Less

Submitted 15 January, 2011; originally announced January 2011.

Comments: Accepted to A&A

Planck Early Results: The submillimetre properties of a sample of Galactic cold clumps

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. Balbi, A. J. Banday, R. B. Barreiro, J. G. Bartlett, E. Battaner, K. Benabed, A. Benoît, J. -P. Bernard, M. Bersanelli, R. Bhatia, J. J. Bock, A. Bonaldi, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana , et al. (179 additional authors not shown)

Abstract: (abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sour… ▽ More (abridged) We perform a detailed investigation of sources from the Cold Cores Catalogue of Planck Objects (C3PO). Our goal is to probe the reliability of the detections, validate the separation between warm and cold dust emission components, provide the first glimpse at the nature, internal morphology and physical characterictics of the Planck-detected sources. We focus on a sub-sample of ten sources from the C3PO list, selected to sample different environments, from high latitude cirrus to nearby (150pc) and remote (2kpc) molecular complexes. We present Planck surface brightness maps and derive the dust temperature, emissivity spectral index, and column densities of the fields. With the help of higher resolution Herschel and AKARI continuum observations and molecular line data, we investigate the morphology of the sources and the properties of the substructures at scales below the Planck beam size. △ Less

Submitted 11 January, 2011; originally announced January 2011.

Comments: This paper is part of a package of papers describing first results of the Planck mission, to be published together in a special issue of A&A. Corresponding author: Isabelle.Ristorcelli@cesr.fr. Full resolution version available from http://www.sciops.esa.int/index.php?project=PLANCK&page=Planck_Published_Papers

Accuracy of core mass estimates in simulated observations of dust emission

Authors: J. Malinen, M. Juvela, D. C. Collins, T. Lunttila, P. Padoan

Abstract: We study the reliability of mass estimates obtained for molecular cloud cores using sub-millimetre and infrared dust emission. We use magnetohydrodynamic simulations and radiative transfer to produce synthetic observations with spatial resolution and noise levels typical of Herschel surveys. We estimate dust colour temperatures using different pairs of intensities, calculate column densities and c… ▽ More We study the reliability of mass estimates obtained for molecular cloud cores using sub-millimetre and infrared dust emission. We use magnetohydrodynamic simulations and radiative transfer to produce synthetic observations with spatial resolution and noise levels typical of Herschel surveys. We estimate dust colour temperatures using different pairs of intensities, calculate column densities and compare the estimated masses with the true values. We compare these results to the case when all five Herschel wavelengths are available. We investigate the effects of spatial variations of dust properties and the influence of embedded heating sources. Wrong assumptions of dust opacity and its spectral index beta can cause significant systematic errors in mass estimates. These are mainly multiplicative and leave the slope of the mass spectrum intact, unless cores with very high optical depth are included. Temperature variations bias colour temperature estimates and, in quiescent cores with optical depths higher than for normal stable cores, masses can be underestimated by up to one order of magnitude. When heated by internal radiation sources the observations recover the true mass spectra. The shape, although not the position, of the mass spectrum is reliable against observational errors and biases introduced in the analysis. This changes only if the cores have optical depths much higher than expected for basic hydrostatic equilibrium conditions. Observations underestimate the value of beta whenever there are temperature variations along the line of sight. A bias can also be observed when the true beta varies with wavelength. Internal heating sources produce an inverse correlation between colour temperature and beta that may be difficult to separate from any intrinsic beta(T) relation of the dust grains. This suggests caution when interpreting the observed mass spectra and the spectral indices. △ Less

Submitted 18 April, 2011; v1 submitted 23 September, 2010; originally announced September 2010.

Comments: Revised version, 17 pages, 17 figures, submitted to A&A

Journal ref: A&A, 530 (2011), A101