Dense Blobs Research Articles

Dropouts of Fully Stripped Ions in the Solar Wind: A Diagnostic for Wave Heating versus Reconnection

The SWICS instrument on board the ACE satellite has detected frequent intervals in the slow solar wind and interplanetary coronal mass ejections in which C6+ and other fully stripped ions are strongly depleted, though the ionization states of elements such as Si and Fe indicate that those ions should be present. It has been suggested that these “outlier” or “dropout” events can be explained by the resonant cyclotron heating process, because these ions all have the same cyclotron frequency as He2+. We investigate the region in the corona where these outlier events form. It must be above the ionization freeze-in height and the transition to collisionless plasma conditions, but low enough that the wind still feels the effects of solar gravity. We suggest that the dropout events correspond to relatively dense blobs of gas in which the heating is reduced because local variations in the Alfvén speed change the reflection of Alfvén waves and the turbulent cascade. As a result, the wave power at the cyclotron frequency of the fully stripped ions is absorbed by He2+ and may not be able to heat the other fully stripped ions enough to overcome solar gravity. If this picture is borne out, it may help to discriminate between resonant cyclotron heating and stochastic heating models of the solar wind.

Exploring γ-Ray Flares in the Long-term Light Curves of CTA 102 at GeV Energies

Blazar CTA 102 experienced an intense multiwavelength activity phase from 2015 to 2018; in particular, an unprecedented outburst was observed from 2016 October to 2017 February. In this work, we extract a 7 day binned γ-ray light curve from 2008 August to 2018 March in the energy range 0.1–300 GeV and identify three main outbursts. We study in detail the short-timescale variability of these three outbursts via an exponential function with parameterized rise and decay timescales. The obtained shortest rise and decay timescales are 0.70 ± 0.05 hr and 0.79 ± 0.27 hr, respectively. Based on these variability timescales, the physical parameters of the flaring region (e.g., the minimum Doppler factor and the emission region size) are constrained. The short-timescale flares exhibit a symmetric temporal profile within the error bars, implying that the rise and decay timescales are dominated by the light-crossing timescale or by disturbances caused by dense plasma blobs passing through the standing shock front in the jet region. We also find that the best-fitting form of the γ-ray spectra during the flare period is a power law with an exponential cutoff. The derived jet parameters from the spectral behavior and the temporal characteristics of the individual flares suggest that the γ-ray emission region is located upstream of the radio core. The extreme γ-ray flare of CTA 102 is likely to have been caused by magnetic reconnection.

Emergent Collective Locomotion in an Active Polymer Model of Entangled Worm Blobs

Numerous worm and arthropod species form physically-connected aggregations in which interactions among individuals give rise to emergent macroscale dynamics and functionalities that enhance collective survival. In particular, some aquatic worms such as the California blackworm (Lumbriculus variegatus) entangle their bodies into dense blobs to shield themselves against external stressors and preserve moisture in dry conditions. Motivated by recent experiments revealing emergent locomotion in blackworm blobs, we investigate the collective worm dynamics by modeling each worm as a self-propelled Brownian polymer. Though our model is two-dimensional, compared to real three-dimensional worm blobs, we demonstrate how a simulated blob can collectively traverse temperature gradientsviathe coupling between the active motion and the environment. By performing a systematic parameter sweep over the strength of attractive forces between worms, and the magnitude of their directed self-propulsion, we obtain a rich phase diagram which reveals that effective collective locomotion emerges as a result of finely balancing a tradeoff between these two parameters. Our model brings the physics of active filaments into a new meso- and macroscale context and invites further theoretical investigation into the collective behavior of long, slender, semi-flexible organisms.

BALQSO spectra explained by shock disruption of galactic clouds

ABSTRACT Blue-shifted broad absorption lines (BALs) detected in quasar’s spectra are indicative of active galactic nuclei (AGN) outflows. We show, using 2D hydrodynamical simulations, that disruption of interstellar clouds by a fast AGN wind can lead to formation of cold, dense high-speed blobs that give rise to broad absorption features in the transmission spectrum of the AGN continuum source. For a wind velocity of 0.1c and sufficiently high cloud density (nc ≳ 104 cm−3; depending on size), absorption troughs with velocities up to about 3000 km s−1can be produced. For slower winds and/or lower cloud density, the anticipated velocity of the absorbing clouds should be smaller.

In Situ Measurements of the Variable Slow Solar Wind near Sector Boundaries

Abstract The release of density structures at the tip of the coronal helmet streamers, likely as a consequence of magnetic reconnection, contributes to the mass flux of the slow solar wind (SSW). In situ measurements in the vicinity of the heliospheric plasma sheet of the magnetic field, protons, and suprathermal electrons reveal details of the processes at play during the formation of density structures near the Sun. In a previous article, we exploited remote-sensing observations to derive a 3D picture of the dynamic evolution of a streamer. We found evidence of the recurrent and continual release of dense blobs from the tip of the streamers. In the present paper, we interpret in situ measurements of the SSW during solar maximum. Through both case and statistical analysis, we show that in situ signatures (magnetic field magnitude, smoothness and rotation, proton density, and suprathermal electrons, in the first place) are consistent with the helmet streamers producing, in alternation, high-density regions (mostly disconnected) separated by magnetic flux ropes (mostly connected to the Sun). This sequence of emission of dense blobs and flux ropes also seems repeated at smaller scales inside each of the high-density regions. These properties are further confirmed with in situ measurements much closer to the Sun using Helios observations. We conclude on a model for the formation of dense blobs and flux ropes that explains both the in situ measurements and the remote-sensing observations presented in our previous studies.

Intermittent electron density and temperature fluctuations and associated fluxes in the Alcator C-Mod scrape-off layer

The Alcator C-Mod mirror Langmuir probe system has been used to sample data time series of fluctuating plasma parameters in the outboard mid-plane far scrape-off layer. We present a statistical analysis of one second long time series of electron density, temperature, radial electric drift velocity and the corresponding particle and electron heat fluxes. These are sampled during stationary plasma conditions in an ohmically heated, lower single null diverted discharge. The electron density and temperature are strongly correlated and feature fluctuation statistics similar to the ion saturation current. Both electron density and temperature time series are dominated by intermittent, large-amplitude burst with an exponential distribution of both burst amplitudes and waiting times between them. The characteristic time scale of the large-amplitude bursts is approximately 15 . Large-amplitude velocity fluctuations feature a slightly faster characteristic time scale and appear at a faster rate than electron density and temperature fluctuations. Describing these time series as a superposition of uncorrelated exponential pulses, we find that probability distribution functions, power spectral densities as well as auto-correlation functions of the data time series agree well with predictions from the stochastic model. The electron particle and heat fluxes present large-amplitude fluctuations. For this low-density plasma, the radial electron heat flux is dominated by convection, that is, correlations of fluctuations in the electron density and radial velocity. Hot and dense blobs contribute only a minute fraction of the total fluctuation driven heat flux.

Excitation and evolution of vertically polarised transverse loop oscillations by coronal rain

Context. Coronal rain is composed of cool dense blobs that form in solar coronal loops and are a manifestation of catastrophic cooling linked to thermal instability. The nature and excitation of oscillations associated with coronal rain is not well understood. Aims. We consider observations of coronal rain in a bid to elucidate the excitation mechanism and evolution of wave characteristics. Methods. We analyse IRIS and Hinode/SOT observations of an oscillating coronal rain event on 17 Aug 2014 and determine the wave characteristics as a function of time using tried and tested time-space analysis techniques. Results. We exploit the seismological capability of the oscillation to deduce the relative rain mass from the oscillation amplitude. This is consistent with the evolution of the oscillation period showing the loop loosing a third of its mass due to falling coronal rain in a 10-15 minute time period. Conclusions. We present first evidence of the excitation of vertically polarised transverse loop oscillations triggered by a catastrophic cooling at the loop top and consistent with two thirds of the loop mass comprising of cool rain mass.

Formation of Galactic Prominence in the Galactic Central Region

Abstract We carried out 2.5-dimensional resistive MHD simulations to study the formation mechanism of molecular loops observed by Fukui et al. in the Galactic central region. Since it is hard to form molecular loops by lifting up dense molecular gas, we study the formation mechanism of molecular gas in rising magnetic arcades. This model is based on the in situ formation model of solar prominences, in which prominences are formed by cooling instability in helical magnetic flux ropes formed by imposing converging and shearing motion at footpoints of the magnetic arch anchored to the solar surface. We extended this model to Galactic center scale (a few hundreds of parsecs). Numerical results indicate that magnetic reconnection taking place in the current sheet that formed inside the rising magnetic arcade creates dense blobs confined by the rising helical magnetic flux ropes. Thermal instability taking place in the flux ropes forms dense molecular filaments floating at high Galactic latitude. The mass of the filament increases with time and can exceed 10 5 M ⊙ .

Magnetic shuffling of coronal downdrafts

Channelled fragmented downflows are ubiquitous in magnetized atmospheres, and have been recently addressed from an observation after a solar eruption. We study the possible back-effect of the magnetic field on the propagation of confined flows. We compare two 3D MHD simulations of dense supersonic plasma blobs downfalling along a coronal magnetic flux tube. In one, the blobs move strictly along the field lines; in the other, the initial velocity of the blobs is not perfectly aligned to the magnetic field and the field is weaker. The aligned blobs remain compact while flowing along the tube, with the generated shocks. The misaligned blobs are disrupted and merged by the chaotic shuffling of the field lines, and structured into thinner filaments; Alfven wave fronts are generated together with shocks ahead of the dense moving front. Downflowing plasma fragments can be chaotically and efficiently mixed if their motion is misaligned to field lines, with broad implications, e.g., disk accretion in protostars, coronal eruptions and rain.

Kinematics of coronal rain in a transversely oscillating loop: Ponderomotive force and rain-excited oscillations

Context. Coronal rain are cool dense blobs that form in solar coronal loops and are a manifestation of catastrophic cooling linked to thermal instability. Once formed, rain falls towards the solar surface at sub-ballistic speeds, which is not well-understood. Pressure forces seem to be the prime candidate to explain this. In many observations rain is accompanied by transverse oscillations and the interaction between the two needs to be explored. Aims. Therefore, an alternative kinematic model for coronal rain kinematics in transversely oscillating loops is developed to understand the physical nature of the observed sub-ballistic falling motion of rain. It explicitly explores the role of the ponderomotive force arising from the transverse oscillation on the rain motion as well as the capacity of rain to excite wave motion. Methods. An analytical model is presented that describes a rain blob guided by the coronal magnetic field supporting a onedimensional shear Alfven wave as a point mass on an oscillating string. The model includes gravity and the ponderomotive force from the oscillation acting on the mass, as well as the inertia of the mass acting on the oscillation. Results. The kinematics of rain in the limit of negligible rain mass are explored and falling and trapped regimes are found, depending on wave amplitude. In the trapped regime for the fundamental mode, the rain blob bounces back and forth around the loop top at a long period inversely proportional to the oscillation amplitude. The model is compared with several observational rain studies, including one in-depth comparison with an observation that shows rain with up-and down bobbing motion. The role of rain inertia in exciting transverse oscillations is explored in inclined loops. Conclusions. It is found that the model requires displacement amplitudes of the transverse oscillation that are typically an order of magnitude larger than observed to explain the measured sub-ballistic motion of the rain. Therefore, it is concluded that the ponderomotive force is not the primary reason for understanding sub-ballistic motion, but it plays a role in cases of large loop oscillations. The appearance of rain causes the excitation of small-amplitude transverse oscillations that may explain observed events and provide a seismological tool to measure rain mass.

Impacts of fragmented accretion streams onto classical T Tauri stars: UV and X-ray emission lines

Context. The accretion process in Classical T Tauri Stars (CTTSs) can be studied through the analysis of some UV and X-ray emission lines which trace hot gas flows and act as diagnostics of the post-shock downfalling plasma. In the UV band, where higher spectral resolution is available, these lines are characterized by rather complex profiles whose origin is still not clear. Aims. We investigate the origin of UV and X-ray emission at impact regions of density structured (fragmented) accretion streams.We study if and how the stream fragmentation and the resulting structure of the post-shock region determine the observed profiles of UV and X-ray emission lines. Methods. We model the impact of an accretion stream consisting of a series of dense blobs onto the chromosphere of a CTTS through 2D MHD simulations. We explore different levels of stream fragmentation and accretion rates. From the model results, we synthesize C IV (1550 {\AA}) and OVIII (18.97 {\AA}) line profiles. Results. The impacts of accreting blobs onto the stellar chromosphere produce reverse shocks propagating through the blobs and shocked upflows. These upflows, in turn, hit and shock the subsequent downfalling fragments. As a result, several plasma components differing for the downfalling velocity, density, and temperature are present altoghether. The profiles of C IV doublet are characterized by two main components: one narrow and redshifted to speed $\approx$ 50 km s$^{-1}$ and the other broader and consisting of subcomponents with redshift to speed in the range 200 $\approx$ 400 km s$^{-1}$. The profiles of OVIII lines appear more symmetric than C IV and are redshifted to speed $\approx$ 150 km s$^{-1}$. Conclusions. Our model predicts profiles of C IV line remarkably similar to those observed and explains their origin in a natural way as due to stream fragmentation.

How the internal structures of the imprinted and the random hydrogels change upon washing?

The effect of washing with acidic and basic solutions on the internal structure (fluctuations in the density and the size of the polymer clusters through the gel) of the imprinted and the random PNIPA (polyN-isopropylacrylamide) hydrogels has been studied by using steady-state fluorescence technique, swelling kinetics and SEM images. When the gel is imprinted, fluctuations in the density and the size of the polymer clusters decrease through the gel; i.e., the gel becomes more homogeneous. Washing with acid causes no considerable change in the internal structure of the both gels, but washing with base causes a drastic change in the structure by breaking the bonds between the polymer strands. The release rate of target molecules for the imprinted gel becomes considerably higher than for the random one regardless of the washing condition since the target molecules are trapped in the dense blobs of the random gel. The change in the internal structure of the gels upon washing is modeled by a spring approximation. The effect of washing with acidic and basic solutions on the morphology of the imprinted and the random PNIPA hydrogels was studied. When the gel is imprinted, it becomes more homogeneous. Washing with acid causes no considerable change in the morphology of the both gels, but washing with base causes a drastic change by breaking the bonds between the polymer strands.

Formation and evolution of coronal rain observed by SDO/AIA on February 22, 2012

The formation and dynamics of coronal rain are currently not fully understood. Coronal rain is the fall of cool and dense blobs formed by thermal instability in the solar corona towards the solar surface with acceleration smaller than gravitational free fall. We aim to study the observational evidence of the formation of coronal rain and to trace the detailed dynamics of individual blobs. We used time series of the 171 \AA\, and 304 \AA\, spectral lines obtained by the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) above active region AR 11420 on February 22, 2012. Observations show that a coronal loop disappeared in the 171 \AA\ channel and appeared in the 304 \AA\ line$\text{}\text{}$ more than one hour later, which indicates a rapid cooling of the coronal loop from 1 MK to 0.05 MK. An energy estimation shows that the radiation is higher than the heat input, which indicates so-called catastrophic cooling. The cooling was accompanied by the formation of coronal rain in the form of falling cold plasma. We studied two different sequences of falling blobs. The first sequence includes three different blobs. The mean velocities of the blobs were estimated to be 50 km s$^{-1}$, 60 km s$^{-1}$ and 40 km s$^{-1}$. A polynomial fit shows the different values of the acceleration for different blobs, which are lower than free-fall in the solar corona. The first and second blob move along the same path, but with and without acceleration, respectively. We performed simple numerical simulations for two consecutive blobs, which show that the second blob moves in a medium that is modified by the passage of the first blob. Therefore, the second blob has a relatively high speed and no acceleration, as is shown by observations. The second sequence includes two different blobs with mean velocities of 100 km s$^{-1}$ and 90 km s$^{-1}$, respectively.

A vortical dawn flank boundary layer for near‐radial IMF: Wind observations on 24 October 2001

AbstractWe present an example of a boundary layer tailward of the dawn terminator which is entirely populated by rolled‐up flow vortices. Observations were made by Wind on 24 October 2001 as the spacecraft moved across the region at X ∼−13 RE. Interplanetary conditions were steady with a near‐radial interplanetary magnetic field (IMF). Approximately 15 vortices were observed over the 1.5 h duration of Wind's crossing, each lasting ∼5 min. The rolling up is inferred from the presence of a hot tenuous plasma being accelerated to speeds higher than in the adjoining magnetosheath, a circumstance which has been shown to be a reliable signature of this in single‐spacecraft observations. A blob of cold dense plasma was entrained in each vortex, at whose leading edge abrupt polarity changes of field and velocity components at current sheets were regularly observed. In the frame of the average boundary layer velocity, the dense blobs were moving predominantly sunward and their scale size along X was ∼7.4 RE. Inquiring into the generation mechanism of the vortices, we analyze the stability of the boundary layer to sheared flows using compressible magnetohydrodynamic Kelvin‐Helmholtz theory with continuous profiles for the physical quantities. We input parameters from (i) the exact theory of magnetosheath flow under aligned solar wind field and flow vectors near the terminator and (ii) the Wind data. It is shown that the configuration is indeed Kelvin‐Helmholtz (KH) unstable. This is the first reported example of KH‐unstable waves at the magnetopause under a radial IMF.

DYNAMICS OF CORONAL RAIN AND DESCENDING PLASMA BLOBS IN SOLAR PROMINENCES. I. FULLY IONIZED CASE

Observations of active regions and limb prominences often show cold, dense blobs descending with an acceleration smaller than that of free fall. The dynamics of these condensations falling in the solar corona is investigated in this paper using a simple fully ionised plasma model. We find that the presence of a heavy condensation gives rise to a dynamical rearrangement of the coronal pressure that results in the formation of a large pressure gradient that opposes gravity. Eventually this pressure gradient becomes so large that the blob acceleration vanishes or even points upwards. Then, the blob descent is characterised by an initial acceleration phase followed by an essentially constant velocity phase. These two stages can be identified in published time-distance diagrams of coronal rain events. Both the duration of the first stage and the velocity attained by the blob increase for larger values of the ratio of blob to coronal density, for larger blob mass, and for smaller coronal temperature. Dense blobs are characterised by a detectable density growth (up to 60% in our calculations) and by a steepening of the density in their lower part, that could lead to the formation of a shock. They also emit sound waves that could be detected as small intensity changes with periods of the order of 100 s and lasting between a few and about ten periods. Finally, the curvature of the falling path is only relevant when a very dense blob falls along inclined magnetic field lines.

Molecular cloud evolution – V. Cloud destruction by stellar feedback

We present a numerical study of the evolution of molecular clouds, from their formation by converging flows in the warm interstellar medium, to their destruction by the ionizing feedback of the massive stars they form. We improve with respect to our previous simulations by including a different stellar-particle formation algorithm, which allows them to have masses corresponding to single stars rather than to small clusters, and with a mass distribution following a near-Salpeter stellar initial mass function. We also employ a simplified radiative-transfer algorithm that allows the stellar particles to feedback on the medium at a rate that depends on their mass and the local density. Our results are as follows: (a) contrary to the results from our previous study, where all stellar particles injected energy at a rate corresponding to a star of ∼10 M⊙, the dense gas is now completely evacuated from 10 pc regions around the stars within 10–20 Myr, suggesting that this feat is accomplished essentially by the most massive stars. (b) At the scale of the whole numerical simulations, the dense gas mass is reduced by up to an order of magnitude, although star formation (SF) never shuts off completely, indicating that the feedback terminates SF locally, but new SF events continue to occur elsewhere in the clouds. (c) The SF efficiency (SFE) is maintained globally at the ∼10 per cent level, although locally, the cloud with largest degree of focusing of its accretion flow reaches SFE ∼30 per cent. (d) The virial parameter of the clouds approaches unity before the stellar feedback begins to dominate the dynamics, becoming much larger once feedback dominates, suggesting that clouds become unbound as a consequence of the stellar feedback, rather than unboundness being the cause of a low SFE. (e) The erosion of the filaments that feed the star-forming clumps produces chains of isolated dense blobs reminiscent of those observed in the vicinity of the dark globule B68.

On-Disk Coronal Rain

Small and elongated, cool and dense blob-like structures are being reported with high resolution telescopes in physically different regions throughout the solar atmosphere. Their detection and the understanding of their formation, morphology, and thermodynamical characteristics can provide important information on their hosting environment, especially concerning the magnetic field, whose understanding constitutes a major problem in solar physics. An example of such blobs is coronal rain, a phenomenon of thermal non-equilibrium observed in active region loops, which consists of cool and dense chromospheric blobs falling along loop-like paths from coronal heights. So far, only off-limb coronal rain has been observed, and few reports on the phenomenon exist. In the present work, several data sets of on-disk Hα observations with the CRisp Imaging SpectroPolarimeter (CRISP) at the Swedish 1-m Solar Telescope (SST) are analyzed. A special family of on-disk blobs is selected for each data set, and a statistical analysis is carried out on their dynamics, morphology, and temperature. All characteristics present distributions which are very similar to reported coronal rain statistics. We discuss possible interpretations considering other similar blob-like structures reported so far and show that a coronal rain interpretation is the most likely one. The chromospheric nature of the blobs and the projection effects (which eliminate all direct possibilities of height estimation) on one side, and their small sizes, fast dynamics, and especially their faint character (offering low contrast with the background intensity) on the other side, are found as the main causes for the absence until now of the detection of this on-disk coronal rain counterpart.

THE CIRCUMSTELLAR ENVIRONMENT OF R CORONAE BOREALIS: WHITE DWARF MERGER OR FINAL-HELIUM-SHELL FLASH?

In 2007, R Coronae Borealis (R CrB) went into an historically deep and long decline. In this state, the dust acts like a natural coronagraph at visible wavelengths, allowing faint nebulosity around the star to be seen. Imaging has been obtained from 0.5 to 500 micron with Gemini/GMOS, HST/WFPC2, Spitzer/MIPS, and Herschel/SPIRE. Several of the structures around R CrB are cometary globules caused by wind from the star streaming past dense blobs. The estimated dust mass of the knots is consistent with their being responsible for the R CrB declines if they form along the line of sight to the star. In addition, there is a large diffuse shell extending up to 4 pc away from the star containing cool 25 K dust that is detected all the way out to 500 micron. The SED of R CrB can be well fit by a 150 AU disk surrounded by a very large diffuse envelope which corresponds to the size of the observed nebulosity. The total masses of the disk and envelope are 10^-4 and 2 M(Sun), respectively, assuming a gas-to-dust ratio of 100. The evidence pointing toward a white-dwarf merger or a final-helium-shell flash origin for R CrB is contradictory. The shell and the cometary knots are consistent with a fossil planetary nebula. Along with the fact that R CrB shows significant Lithium in its atmosphere, this supports the final-helium-shell flash. However, the relatively high inferred mass of R CrB and its high fluorine abundance support a white-dwarf merger.

Entropy mode at a magnetic null point as a possible tool for indirect observation of nanoflares in the solar corona

Aims. We aim to explore the dynamics of the entropy mode perturbation excited by an energy release in the vicinity of a magnetic null point that is embedded in a gravitationally stratified solar corona. Methods. We solve two-dimensional, time-dependent magnetohydrodynamic equations numerically to find spatial and temporal signatures of the entropy mode that is triggered impulsively by a spatially localized pulse of the gas pressure. Results. We find that the properties of the entropy mode are determined by the sign of the initial pressure pulse. The initial increase in the gas pressure creates, together with the magnetoacoustic-gravity waves, a stationary void of the rarefied plasma at the launching place, associated with the entropy mode. In contrast, an initial decrease in the gas pressure, which corresponds to a rapid (or catastrophic) cooling, forms a blob of the dense plasma at the launching place. Conclusions. The cool, dense blobs at magnetic null points may be observed in transition region and chromospheric spectral lines at coronal heights off the solar limb and may be associated with the places of nanoflare occurrence. Therefore, extensions of entropy mode studies may produce a diagnostic tool for indirect observations of nanoflares. The dense cool blobs may be affected by the gravity or carried by downflows, hence may initiate a coronal rain.

3D Magneto-hydrodynamic simulations of Parker instability with cosmic rays

This study investigates Parker instability in an interstellar medium (ISM) near the Galactic plane using three-dimensional magneto-hydrodynamic simulations. Parker instability arises from the presence of a magnetic field in a plasma, wherein the magnetic buoyant pressure expels the gas and cause the gas to move along the field lines. The process is thought to induce the formation of giant molecular clouds in the Galaxy. In this study, the effects of cosmic-ray (CR) diffusion are examined. The ISM at equilibrium is assumed to comprise a plasma fluid and a CR fluid at various temperatures, with a uniform magnetic field passing through it in the azimuthal direction of the Galactic disk. After a small perturbation, the unstable gas aggregates at the footpoint of the magnetic fields and forms dense blobs. The growth rate of the instability increases with the strength of the CR diffusion. The formation of dense clouds is enhanced by the effect of cosmic rays (CRs), whereas the shape of the clouds depends sensitively on the initial conditions of perturbation.