• ### Stochastic modeling of multiwavelength variability of the classical BL Lac object OJ 287 on timescales ranging from decades to hours(1709.04457)

July 10, 2018 astro-ph.HE
We present the results of our power spectral density analysis for the BL Lac object OJ\,287, utilizing the {\it Fermi}-LAT survey at high-energy $\gamma$-rays, {\it Swift}-XRT in X-rays, several ground-based telescopes and the {\it Kepler} satellite in the optical, and radio telescopes at GHz frequencies. The light curves are modeled in terms of continuous-time auto-regressive moving average (CARMA) processes. Owing to the inclusion of the {\it Kepler} data, we were able to construct \emph{for the first time} the optical variability power spectrum of a blazar without any gaps across $\sim6$ dex in temporal frequencies. Our analysis reveals that the radio power spectra are of a colored-noise type on timescales ranging from tens of years down to months, with no evidence for breaks or other spectral features. The overall optical power spectrum is also consistent with a colored noise on the variability timescales ranging from 117 years down to hours, with no hints of any quasi-periodic oscillations. The X-ray power spectrum resembles the radio and optical power spectra on the analogous timescales ranging from tens of years down to months. Finally, the $\gamma$-ray power spectrum is noticeably different from the radio, optical, and X-ray power spectra of the source: we have detected a characteristic relaxation timescale in the {\it Fermi}-LAT data, corresponding to $\sim 150$\,days, such that on timescales longer than this, the power spectrum is consistent with uncorrelated (white) noise, while on shorter variability timescales there is correlated (colored) noise.
• ### Primary black hole spin in OJ287 as determined by the General Relativity centenary flare(1603.04171)

March 14, 2016 astro-ph.HE
OJ287 is a quasi-periodic quasar with roughly 12 year optical cycles. It displays prominent outbursts which are predictable in a binary black hole model. The model predicted a major optical outburst in December 2015. We found that the outburst did occur within the expected time range, peaking on 2015 December 5 at magnitude 12.9 in the optical R-band. Based on Swift/XRT satellite measurements and optical polarization data, we find that it included a major thermal component. Its timing provides an accurate estimate for the spin of the primary black hole, chi = 0.313 +- 0.01. The present outburst also confirms the established general relativistic properties of the system such as the loss of orbital energy to gravitational radiation at the 2 % accuracy level and it opens up the possibility of testing the black hole no-hair theorem with a 10 % accuracy during the present decade.
• ### New transit observations for HAT-P-30 b, HAT-P-37 b, TrES-5 b, WASP-28 b, WASP-36 b, and WASP-39 b(1603.03268)

March 10, 2016 astro-ph.EP
We present new transit light curves for planets in six extrasolar planetary systems. They were acquired with 0.4-2.2 m telescopes located in west Asia, Europe, and South America. When combined with literature data, they allowed us to redetermine system parameters in a homogeneous way. Our results for individual systems are in agreement with values reported in previous studies. We refined transit ephemerides and reduced uncertainties of orbital periods by a factor between 2 and 7. No sign of any variations in transit times was detected for the planets studied.
• ### Departure from the constant-period ephemeris for the transiting exoplanet WASP-12 b(1602.09055)

Feb. 29, 2016 astro-ph.EP
Most hot Jupiters are expected to spiral in towards their host stars due to transfering of the angular momentum of the orbital motion to the stellar spin. Their orbits can also precess due to planet-star interactions. Calculations show that both effects could be detected for the very-hot exoplanet WASP-12 b using the method of precise transit timing over a timespan of the order of 10 yr. We acquired new precise light curves for 29 transits of WASP-12 b, spannning 4 observing seasons from November 2012 to February 2016. New mid-transit times, together with literature ones, were used to refine the transit ephemeris and analyse the timing residuals. We find that the transit times of WASP-12 b do not follow a linear ephemeris with a 5 sigma confidence level. They may be approximated with a quadratic ephemeris that gives a rate of change in the orbital period of -2.56 +/- 0.40 x 10^{-2} s/yr. The tidal quality parameter of the host star was found to be equal to 2.5 x 10^5 that is comparable to theoretical predictions for Sun-like stars. We also consider a model, in which the observed timing residuals are interpreted as a result of the apsidal precession. We find, however, that this model is statistically less probable than the orbital decay.
• ### 2MASSJ22560844+5954299: the newly discovered cataclysmic star with the deepest eclipse(1510.04094)

Oct. 14, 2015 astro-ph.SR
Context: The SW Sex stars are assumed to represent a distinguished stage in CV evolution, making it especially important to study them. Aims: We discovered a new cataclysmic star and carried out prolonged and precise photometric observations, as well as medium-resolution spectral observations. Modelling these data allowed us to determine the psysical parameters and to establish its peculiarities. Results: The newly discovered vataclysmic variable 2MASSJ22560844+5954299 shows the deepest eclipse amongst the known nova-like stars. It was reproduced by totally covering a very luminous accretion disk by a red secondary component. The temperature distribution of the disk is flatter than that of steady-state disk. The target is unusual with the combination of a low mass ratio q~1.0 (considerably below the limit q=1.2 of stable mass transfer of CVs) and an M-star secondary. The intensity of the observed three emission lines, H_alpha, He 5875, and He 6678, sharply increases around phase 0.0, accompanied by a Doppler jump to the shorter wavelength. The absence of eclipses of the emission lines and their single-peaked profiles means that they originate mainly in a vertically extended hot-spot halo. The emission H_alpha line reveals S-wave wavelength shifts with semi-amplitude of around 210 km/s and phase lag of 0.03. Conclusions: The non-steady-state emission of the luminous accretion disk of 2MASSJ22560844+5954299 was attributed to the low viscosity of the disk matter caused by its unusually high temperature. The star shows all spectral properties of an SW Sex variable apart from the 0.5 central absorption.
• ### Ground-based transit observations of the HAT-P-18, HAT-P-19, HAT-P-27/WASP-40 and WASP-21 systems(1508.06215)

Aug. 25, 2015 astro-ph.SR, astro-ph.EP
As part of our ongoing effort to investigate transit timing variations (TTVs) of known exoplanets, we monitored transits of the four exoplanets HAT-P-18b, HAT-P-19b, HAT-P-27b/WASP-40b and WASP-21b. All of them are suspected to show TTVs due to the known properties of their host systems based on the respective discovery papers. During the past three years 46 transit observations were carried out, mostly using telescopes of the Young Exoplanet Transit Initiative. The analyses are used to refine the systems orbital parameters. In all cases we found no hints for significant TTVs, or changes in the system parameters inclination, fractional stellar radius and planet to star radius ratio. However, comparing our results with those available in the literature shows that we can confirm the already published values.
• ### No variations in transit times for Qatar-1 b(1503.07191)

March 24, 2015 astro-ph.EP
The transiting hot Jupiter planet Qatar-1 b was presented to exhibit variations in transit times that could be of perturbative nature. A hot Jupiter with a planetary companion on a nearby orbit would constitute an unprecedented planetary configuration, important for theories of formation and evolution of planetary systems. We performed a photometric follow-up campaign to confirm or refute transit timing variations. We extend the baseline of transit observations by acquiring 18 new transit light curves acquired with 0.6-2.0 m telescopes. These photometric time series, together with data available in the literature, were analyzed in a homogenous way to derive reliable transit parameters and their uncertainties. We show that the dataset of transit times is consistent with a linear ephemeris leaving no hint for any periodic variations with a range of 1 min. We find no compelling evidence for the existence of a close-in planetary companion to Qatar-1 b. This finding is in line with a paradigm that hot Jupiters are not components of compact multi-planetary systems. Based on dynamical simulations, we place tighter constraints on a mass of any fictitious nearby planet in the system. Furthermore, new transit light curves allowed us to redetermine system parameters with the precision better than that reported in previous studies. Our values generally agree with previous determinations.
• ### The unusual photometric variability of the PMS star GM Cep(1502.07579)

Feb. 26, 2015 astro-ph.SR
Results from UBVRI photometric observations of the pre-main sequence star GM Cep obtained in the period April 2011 - August 2014 are reported in the paper. Presented data are a continuation of our photometric monitoring of the star started in 2008. GM Cep is located in the field of the young open cluster Trumpler 37 and over the past years it has been an object of intense photometric and spectral studies. The star shows a strong photometric variability interpreted as a possible outburst from EXor type in previous studies. Our photometric data for a period of over six years show a large amplitude variability (Delta V ~ 2.3 mag.) and several deep minimums in brightness are observed. The analysis of the collected multicolor photometric data shows the typical of UX Ori variables a color reversal during the minimums in brightness. The observed decreases in brightness have a different shape, and evidences of periodicity are not detected. At the same time, high amplitude rapid variations in brightness typical for the classical T Tauri stars also present on the light curve of GM Cep. The spectrum of GM Cep shows the typical of classical T Tauri stars wide H/alpha emission line and absorption lines of some metals. We calculate the outer radius of the H/alpha emitting region as 10.4 +/-0.5 Rsun and the accretion rate as 1.8 x 10 E-7 Msun/yr.
• ### The HU Aqr planetary system hypothesis revisited(1412.5899)

Dec. 18, 2014 astro-ph.EP
We study the mid-egress eclipse timing data gathered for the cataclysmic binary HU Aquarii during the years 1993-2014. The (O-C) residuals were previously attributed to a single ~7 Jupiter mass companion in ~5 au orbit or to a stable 2-planet system with an unconstrained outermost orbit. We present 22 new observations gathered between June, 2011 and July, 2014 with four instruments around the world. They reveal a systematic deviation of ~60 - 120 seconds from the older ephemeris. We re-analyse the whole set of the timing data available. Our results provide an erratum to the previous HU Aqr planetary models, indicating that the hypothesis for a third and fourth body in this system is uncertain. The dynamical stability criterion and a particular geometry of orbits rule out coplanar 2-planet configurations. A putative HU Aqr planetary system may be more complex, e.g., highly non-coplanar. Indeed, we found examples of 3-planet configurations with the middle planet in a retrograde orbit, which are stable for at least 1Gyr, and consistent with the observations. The (O-C) may be also driven by oscillations of the gravitational quadrupole moment of the secondary, as predicted by the Lanza et al. modification of the Applegate mechanism. Further systematic, long-term monitoring of HU Aqr is required to interpret the (O-C) residuals.
• ### The 2014 Eclipse of EE Cep. Announcement for a Third International Observational Campaign(1407.1541)

July 6, 2014 astro-ph.SR
EE Cep is a unique system in which a Be star is eclipsed by a dark dusty disk, making this star similar to the famous epsilon Aur in many respects. The depth and the duration of the EE Cep eclipses change to a large extent. The last two eclipses were observed in the framework of extensive international campaigns. The joint analysis of these campaigns data and historical photometry, enabled us to propose a model of this system, which implies a disk precession with a period approximately 11-12 times larger than the orbital period. This model predicts that the forthcoming eclipse should be among the deepest observed, reaching about 2 mag. The next eclipse approaches - the photometric minimum should occur around August 23, 2014. Here we would like to announce a new, third international campaign with purpose to verify the disk precession model and to put more constraints on the physical parameters of this system.
• ### Transit Timing Analysis in the HAT-P-32 system(1404.6354)

April 25, 2014 astro-ph.SR, astro-ph.EP
We present the results of 45 transit observations obtained for the transiting exoplanet HAT-P-32b. The transits have been observed using several telescopes mainly throughout the YETI network. In 25 cases, complete transit light curves with a timing precision better than $1.4\:$min have been obtained. These light curves have been used to refine the system properties, namely inclination $i$, planet-to-star radius ratio $R_\textrm{p}/R_\textrm{s}$, and the ratio between the semimajor axis and the stellar radius $a/R_\textrm{s}$. First analyses by Hartman et al. (2011) suggest the existence of a second planet in the system, thus we tried to find an additional body using the transit timing variation (TTV) technique. Taking also literature data points into account, we can explain all mid-transit times by refining the linear ephemeris by 21ms. Thus we can exclude TTV amplitudes of more than $\sim1.5$min.
• ### Photometric reverberation mapping of Markarian 279(1403.4753)

March 19, 2014 astro-ph.CO, astro-ph.GA
By using standard broad-band VRI photometry we were able to discriminate the variations of the broad hydrogen alpha line from the continuum variations for the active galaxy Mkn 279. Cross-correlating both light curves enabled us to determine the time lag of the broad line variations behind the continuum and thus to determine the BLR size (about 8 light days). Our preliminary results are rather consistent with the spectroscopic reverberation mapping results (about 12/17 days). This study is a part of an ambitious program to perform photometric reverberation mapping and determine BLR sizes (respectively - the central black hole masses) for more that 100 nearby AGN.
• ### Constraints on a second planet in the WASP-3 system(1309.6733)

Sept. 26, 2013 astro-ph.EP
There have been previous hints that the transiting planet WASP-3 b is accompanied by a second planet in a nearby orbit, based on small deviations from strict periodicity of the observed transits. Here we present 17 precise radial velocity measurements and 32 transit light curves that were acquired between 2009 and 2011. These data were used to refine the parameters of the host star and transiting planet. This has resulted in reduced uncertainties for the radii and masses of the star and planet. The radial-velocity data and the transit times show no evidence for an additional planet in the system. Therefore, we have determined the upper limit on the mass of any hypothetical second planet, as a function of its orbital period.
• ### Seismic investigation of the {\gamma} Dor star KIC 6462033: The first results of Kepler and ground-based follow up observations(1309.3129)

Sept. 12, 2013 astro-ph.SR
We present the first preliminary results on the analysis of ground-based time series of the {\gamma} Dor star KIC 6462033 (TYC 3144-646-1, V = 10.83, P = 0.69686 d) as well as Kepler photometry in order to study pulsational behaviour in this star.{\gamma} Dor variables, which exhibit g-mode pulsations, are promising asteroseismic targets to understand their rich complexity of pulsational characteristics in detail. In order to achieve this goal, intensive and numerous multicolour and high resolution spectroscopic observations are also required, to complete space-based data aimed at the determination of their physical parameters.
• Compilation of CTA contributions to the proceedings of the 33rd International Cosmic Ray Conference (ICRC2013), which took place in 2-9 July, 2013, in Rio de Janeiro, Brazil
• ### Multi-site photometric campaign on the high amplitude {\delta} Scuti star KIC 6382916(1304.6993)

April 25, 2013 astro-ph.SR
We present results of a multi-site photometric campaign on the high-amplitude $\delta$\,Scuti star KIC\,6382916 in the {\it Kepler} field. The star was observed over a 85-d interval at five different sites in North America and Europe during 2011. {\it Kepler} photometry and ground-based multicolour light curves of KIC\,6382916 are used to investigate the pulsational content and to identify the principal modes. High-dispersion spectroscopy was also obtained in order to derive the stellar parameters and projected rotational velocity. From an analysis of the {\it Kepler} time series, three independent frequencies and a few hundred combination frequencies are found. The light curve is dominated by two modes with frequencies $f_{1}$= 4.9107 and $f_{2}$= 6.4314\,d$^{-1}$. The third mode with $f_{3}$= 8.0350\,d$^{-1}$ has a much lower amplitude. We attempt mode identification by examining the amplitude ratios and phase differences in different wavebands from multicolour photometry and comparing them to calculations for different spherical harmonic degree, $l$. We find that the theoretical models for $f_1$ and $f_2$ are in a best agreement with the observations and lead to value of l = 1 modes, but the mode identification of $f_3$ is uncertain due to its low amplitude. Non-adiabatic pulsation models show that frequencies below 6\,d$^{-1}$ are stable, which means that the low frequency of $f_1$ cannot be reproduced. This is further confirmation that current models predict a narrower pulsation frequency range than actually observed.
• ### New approach for modeling of transiting exoplanets for arbitrary limb-darkening law(1303.2573)

March 11, 2013 astro-ph.EP, astro-ph.IM
We present a new solution of the direct problem of planet transits based on transformation of double integrals to single ones. On the basis of our direct problem solution we created the code TAC-maker for rapid and interactive calculation of synthetic planet transits by numerical computations of the integrals. The validation of our approach was made by comparison with the results of the wide-spread Mandel & Agol (2002) method for the cases of linear, quadratic and squared root limb-darkening laws and various combinations of model parameters. For the first time our approach allows the use of arbitrary limb-darkening law of the host star. This advantage together with the practically arbitrary precision of the calculations make the code a valuable tool that faces the challenges of the continuously increasing photometric precision of the ground-based and space observations.
• ### The SW Sex-type star 2MASS J01074282+4845188: an unusual bright accretion disk with non-steady emission and a hot white dwarf(1303.1675)

March 7, 2013 astro-ph.SR
We present new photometric and spectral observations of the newly discovered nova-like eclipsing star 2MASS J01074282+4845188. To obtain a light curve solution we used model of a nova-like star whose emission sources are a white dwarf surrounded by an accretion disk, a secondary star filling its Roche lobe, a hot spot and a hot line. 2MASS J01074282+4845188 shows the deepest permanent eclipse among the known nova-like stars. It is reproduced by covering the very bright accretion disk by the secondary component. The luminosity of the disk is much bigger than that of the rest light sources. The determined high temperature of the disk is typical for that observed during the outbursts of CVs. The primary of 2MASS J01074282+4845188 is one of the hottest white dwarfs in CVs. The temperature of 5090 K of its secondary is also quite high and more appropriate for a long-period SW Sex star. It might be explained by the intense heating from the hot white dwarf and the hot accretion disk of the target. The high mass accretion rate $\dot{M} = 8\times 10^{-9} M_{\sun}$ yr$^{-1}$, the broad and single-peaked H$\alpha$ emission profile, and the presence of an S-wave are sure signs for the SW Sex classification of 2MASS J01074282+4845188. The obtained flat temperature distribution along the disk radius as well as the deviation of the energy distribution from the black-body law are evidence of the non-steady emission of the disk. The close values of the disk temperature and the parameter $\alpha_{g}$ of 2MASS J01074282+4845188 and those of the cataclysmic stars at eruptions might be considered as an additional argument for the permanent active state of nova-like stars.
• ### Multi-site campaign for transit timing variations of WASP-12 b: possible detection of a long-period signal of planetary origin(1301.5976)

Jan. 25, 2013 astro-ph.EP
The transiting planet WASP-12 b was identified as a potential target for transit timing studies because a departure from a linear ephemeris was reported in the literature. Such deviations could be caused by an additional planet in the system. We attempt to confirm the existence of claimed variations in transit timing and interpret its origin. We organised a multi-site campaign to observe transits by WASP-12 b in three observing seasons, using 0.5-2.6-metre telescopes. We obtained 61 transit light curves, many of them with sub-millimagnitude precision. The simultaneous analysis of the best-quality datasets allowed us to obtain refined system parameters, which agree with values reported in previous studies. The residuals versus a linear ephemeris reveal a possible periodic signal that may be approximated by a sinusoid with an amplitude of 0.00068+/-0.00013 d and period of 500+/-20 orbital periods of WASP-12 b. The joint analysis of timing data and published radial velocity measurements results in a two-planet model which better explains observations than single-planet scenarios. We hypothesize that WASP-12 b might be not the only planet in the system and there might be the additional 0.1 M_Jup body on a 3.6-d eccentric orbit. A dynamical analysis indicates that the proposed two-planet system is stable over long timescales.
• ### Mode identification in the high-amplitude {\delta} Scuti star V2367 Cyg(1210.7147)

Oct. 26, 2012 astro-ph.SR
We report on a multi-site photometric campaign on the high-amplitude $\delta$ Scuti star V2367 Cyg in order to determine the pulsation modes. We also used high-dispersion spectroscopy to estimate the stellar parameters and projected rotational velocity. Time series multicolour photometry was obtained during a 98-d interval from five different sites. These data were used together with model atmospheres and non-adiabatic pulsation models to identify the spherical harmonic degree of the three independent frequencies of highest amplitude as well as the first two harmonics of the dominant mode. This was accomplished by matching the observed relative light amplitudes and phases in different wavebands with those computed by the models. In general, our results support the assumed mode identifications in a previous analysis of Kepler data.
• ### International observational campaigns of the last two eclipses in EE Cephei: 2003 and 2008/9(1205.0028)

Sept. 26, 2012 astro-ph.SR
Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete set of observational data collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous with a "bump" that is very clear in all the UBVRI light curves. Variations in the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and strong absorption is superimposed on the H_alpha emission. These observations confirm that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the "bump" observed during the last two eclipses in terms of the different times of obscuration of the hot polar regions of the Be star by the disk; and (iii) our suggested period of the disk precession (~11-12 P_orb) and predicted depth of about 2 mag the forthcoming eclipse in 2014.
• ### The DWARF project: Eclipsing binaries - precise clocks to discover exoplanets(1206.6709)

June 28, 2012 astro-ph.SR
We present a new observational campaign, DWARF, aimed at detection of circumbinary extrasolar planets using the timing of the minima of low-mass eclipsing binaries. The observations will be performed within an extensive network of relatively small to medium-size telescopes with apertures of ~20-200 cm. The starting sample of the objects to be monitored contains (i) low-mass eclipsing binaries with M and K components, (ii) short-period binaries with sdB or sdO component, and (iii) post-common-envelope systems containing a WD, which enable to determine minima with high precision. Since the amplitude of the timing signal increases with the orbital period of an invisible third component, the timescale of project is long, at least 5-10 years. The paper gives simple formulas to estimate suitability of individual eclipsing binaries for the circumbinary planet detection. Intrinsic variability of the binaries (photospheric spots, flares, pulsation etc.) limiting the accuracy of the minima timing is also discussed. The manuscript also describes the best observing strategy and methods to detect cyclic timing variability in the minima times indicating presence of circumbinary planets. First test observation of the selected targets are presented.
• ### Multi-site, multi-year monitoring of the oscillating Algol-type eclipsing binary CT Her(1105.5064)

May 26, 2011 astro-ph.SR
We present the results of a multi-site photometric campaign carried out in 2004-2008 for the Algol-type eclipsing binary system CT Her, the primary component of which shows Delta Scuti-type oscillations. Our data consist of differential light curves collected in the filters B and V which have been analysed using the method of Wilson-Devinney (PHOEBE). After identification of an adequate binary model and removal of the best-matching light curve solution, we performed a Fourier analysis of the residual B and V light curves to investigate the pulsational behaviour. We confirm the presence of rapid pulsations with a main period of 27.2 min. Up to eight significant frequencies with semi-amplitudes in the range 3 to 1 mmag were detected, all of which surprisingly lie in the frequency range 43.5-53.5 c\d. This result is independent from the choice of the primary's effective temperature (8200 or 8700 K) since the light curve models for the binary are very similar in both cases. This is yet another case of a complex frequency spectrum observed for an accreting Delta Scuti-type star (after Y Cam). In addition, we demonstrate that the amplitudes of several of these pulsation frequencies show evidence of variability on time scales as short as 1-2 years, perhaps even less. Moreover, our analysis takes into account some recently acquired spectra, from which we obtained the corresponding radial velocities for the years 2007-2009. Investigation of the O-C diagram shows that further monitoring of the epochs of eclipse minima of CT Her will cast a new light on the evolution of its orbital period.
• ### UBVRI observations of the flickering of RS Ophiuchi at Quiescence(1001.0852)

Jan. 6, 2010 astro-ph.SR
We report observations of the flickering variability of the recurrent nova RS Oph at quiescence on the basis of simultaneous observations in 5 bands (UBVRI). RS Oph has flickering source with (U-B)_0=-0.62 \pm 0.07, (B-V)_0=0.15 \pm 0.10, (V-R)_0=0.25 \pm 0.05. We find for the flickering source a temperature T_fl = 9500 \pm 500 K, and luminosity L_fl = 50 - 150 L_sun (using a distance of d=1.6kpc). We also find that on a (U-B) vs (B-V) diagram the flickering of the symbiotic stars differs from that of the cataclysmic variables. The possible source of the flickering is discussed. The data are available upon request from the authors and on the web www.astro.bas.bg/~rz/RSOph.UBVRI.2010.MNRAS.tar.gz.
• ### Capillary filling with wall corrugations] Capillary filling in microchannels with wall corrugations: A comparative study of the Concus-Finn criterion by continuum, kinetic and atomistic approaches(0906.3635)

Oct. 12, 2009 physics.flu-dyn
We study the impact of wall corrugations in microchannels on the process of capillary filling by means of three broadly used methods - Computational Fluid Dynamics (CFD), Lattice-Boltzmann Equations (LBE) and Molecular Dynamics (MD). The numerical results of these approaches are compared and tested against the Concus-Finn (CF) criterion, which predicts pinning of the contact line at rectangular ridges perpendicular to flow for contact angles theta > 45. While for theta = 30, theta = 40 (no flow) and theta = 60 (flow) all methods are found to produce data consistent with the CF criterion, at theta = 50 the numerical experiments provide different results. Whilst pinning of the liquid front is observed both in the LB and CFD simulations, MD simulations show that molecular fluctuations allow front propagation even above the critical value predicted by the deterministic CF criterion, thereby introducing a sensitivity to the obstacle heigth.