the dynamics of micron and nanoscale cantilevers are important to a wide variety of technologies . for example , the invention of the atomic force microscope ( afm )  @xcite , which relies upon the dynamics of a cantilever a few hundred microns in length , has revolutionized surface science paving the way for direct measurements of intermolecular forces and topographical mappings with atomic precision for a broad array of materials including semiconductors , polymers , carbon nanotubes , and biological cells  @xcite ( see  @xcite for current reviews ) . in conventional dynamic atomic force microscopy
the cantilever is used to measure the force interactions between the cantilever tip and sample .
cantilevers smaller than conventional afm have also been used to unfold single protein molecules with improved force and time resolution  @xcite .
it has also been proposed to exploit the inherent thermal motion of small cantilevers to make dynamic measurements of single molecules  @xcite .
a passive undriven cantilever placed in a viscous fluid will exhibit stochastic oscillations caused by the thermal bombardment of fluid molecules by brownian motion .
in fact , measuring the thermal spectra of the cantilever is a commonly used afm calibration technique  @xcite .    in all of these applications
the ultimate force sensitivity of a particular measurement is limited by the inherent thermal noise of the experimental system .
there are at least two ways to improve upon this limitation : ( i ) make the cantilevers smaller  @xcite , ( ii ) use correlated measurements of multiple cantilevers  @xcite .    uniformly decreasing the dimensions of a cantilever results in the favorable combination of decreasing the cantilever s equivalent spring constant while increasing its resonant frequency yielding improved force sensitivity and time resolution . by measuring the cross - correlations between two cantilevers in fluid the independent fluctuations of the two
do not contribute , leaving only the smaller correlated fluctuations due to coupling through the fluid .
this type of approach has been used to measure femtonewton forces on millisecond time scales between two micron scale beads placed in separate optical traps ( an improvement of a hundredfold from prior measurements )  @xcite .
the ability to significantly increase the force resolution by making correlated measurements has yet to be exploited for micron and nanoscale cantilevers .
the magnitude of the fluid coupled noise will depend upon the spacing and exact geometries of the cantilevers . combining ( i ) and
( ii ) and measuring the correlations of multiple nanoscale cantilevers offer the potential for experimental measurements with unprecedented force and time resolution  @xcite .    as experimental measurement continues to push toward the stochastic limit it is important that we build a physical understanding of the stochastic dynamics of micron and nanoscale cantilevers for the precise conditions of experiment including complex cantilever geometries  @xcite , the effects of nearby walls  @xcite , and the fluid - coupled dynamics of multiple cantilevers in an array configuration  @xcite .
although micron and nanoscale cantilevers exhibit stochastic motion due to the thermal motion of matter the elastic structures are still large compared to individual fluid molecules and the equations of continuum mechanics remain valid . in what follows we
are concerned with situations where the knudsen number ( the ratio of the mean free path of the fluid particles to the width of a cantilever ) remains sufficiently small so that this statement remains true .
this means that the fluid is described by the usual navier - stokes equations with no - slip and stress continuity boundary conditions at the solid surfaces .
the navier - stokes equations governing the motion of an incompressible fluid , and written in nondimensional form , are @xmath0 where @xmath1 is the fluid velocity , @xmath2 is the pressure , and @xmath3 is time .
there are two inertial terms on the left hand side of eq .
( [ eq : ns ] ) multiplied by the nondimensional parameters @xmath4 and @xmath5 .
the strouhal number @xmath6 plays the role of a frequency based reynolds number expressing the ratio between local inertia forces and viscous forces where @xmath7 and @xmath8 are characteristic length and time scales , respectively .
the velocity based reynolds number @xmath9 expresses the ratio between convective inertial forces and viscous forces .
micron and nanoscale cantilevers are characterized by high oscillation frequencies and small oscillation amplitudes . in this case
@xmath10 so that the nonlinear convective inertial term @xmath11 is negligible and the equations become linear .
however , the strouhal number , @xmath12 ( using the half - width @xmath13 as the appropriate length scale ) is often not negligible . as a result
, the local inertia term must be kept in eq .
( [ eq : ns ] ) making the analysis more difficult .
in addition , experimentally motivated cantilevers are often of complex geometry , near surfaces , or in an array configuration with multiple cantilevers in close proximity .
these difficulties have led to the development of a thermodynamic approach to calculate the stochastic dynamics of micron and submicron scale cantilevers for the precise conditions of experiment ( discussed below )  @xcite . in the following it
is assumed that @xmath5 is negligible , and we leave off the subscript denoting @xmath14 by @xmath15 .
the fluid equations are coupled with the governing equations of elasticity , @xmath16 where @xmath17 is the stress tensor , @xmath18 is the cantilever density , and @xmath19 is the cantilever deflection .
equations  ( [ eq : ns]),([eq : mass ] ) , and ( [ eq : beam ] ) represent the governing deterministic continuum equations .
if one considers lengths scales on the order of a few atomic lengths these equations would include stochastic terms  @xcite .
this represents a difficult fluid - solid interaction problem governing the stochastic cantilever dynamics . a brute force molecular dynamics approach which would resolve the stochastic motion of all of the fluid and solid molecules is computationally prohibitive
however , the system is always near equilibrium permitting a much more accessible thermodynamic based solution strategy .
the paper is organized as follows .
a very general thermodynamic approach , valid for complex geometries as well as for arrays of closely spaced cantilevers , is discussed in section  [ section : thermodynamic_approach ] .
analytical theory based upon the thermodynamic approach , valid for long slender cantilevers , is discussed in section  [ section : analytical_theory ] . in section  [ section : micron_scale_beam ] a micron scale cantilever , similar to an atomic force microscope , is explored and in section  [ section : nanoscale_cantileves ] the stochastic dynamics of an array of experimentally motivated nanoscale cantilevers are quantified using the thermodynamic approach with deterministic finite element numerical simulations to generate results for the precise conditions of experiment
. a general discussion of the approach can be found in ref .
@xcite .    .
in the thermodynamic approach of section  [ section : thermodynamic_approach ] a step force @xmath20 is removed from the cantilever on the left and the deterministic cantilever deflections @xmath21 and @xmath22 are calculated.,width=288 ]
the thermodynamic approach discussed here is based upon the fluctuation - dissipation theorem which states that for equilibrium systems the manner in which the system returns from a linear macroscopic perturbation is related to time correlations of equilibrium microscopic fluctuations  @xcite .
underlying this statement is the fact that equilibrium fluctuations and the dissipation of the system responding to a macroscopic perturbation are governed by the same physics . for the case of miniature elastic objects in a fluid
the dissipation is dominated by the fluid viscosity and the fluctuations by the motion due to the bombardment by the fluid molecules . in what follows ,
we assume that all of the dissipation comes from the fluid and that elastic dissipation in the cantilever is negligible .
this assumption , however , is not required and other sources of dissipation could be included if desired . the essence of our approach is to calculate the dissipation _ deterministically _ and to use the fluctuation - dissipation theorem to determine the cantilever s dynamics . the approach is exact and the only assumptions that have been made are classical dynamics and linear perturbations from equilibrium .
the deterministic calculation of the dissipation can come from analytical theory , simplified models , or from detailed numerical simulations .
the major benefit of this approach is that the deterministic calculations are straightforward , not computationally prohibitive , and methods of calculation are sophisticated and readily available .
we will introduce the use of the thermodynamic approach for the case of two opposing cantilevers as shown in fig .
[ fig : adjacent ] .
consider one dynamical variable to be the stochastic displacement of the cantilever on the left @xmath23 where @xmath23 is a function of the microscopic phase space variables consisting of @xmath24 coordinates , @xmath25 , and conjugate momenta , @xmath26 , of the system and @xmath27 is the number of particles in the system .
the statistical treatment that follows determines the ensemble average of the cantilever deflections @xmath23 over all experimental possibilities . for the undisturbed equilibrium system
the ensemble average , and correlations , are denoted as @xmath28 .
we now take the system to a prescribed macroscopic excursion from equilibrium and observe how the system returns to equilibrium .
the macroscopic cantilever deflection as it returns to equilibrium from a prepared initial condition will be denoted as @xmath21 . for clarity of presentation ,
lower case variables are reserved for stochastic quantities and upper case variables are for deterministic quantities .
the connection between the stochastic and deterministic quantities is particularly straightforward when the excursion from equilibrium is achieved through the application of a step force some time in the distant past that is removed at time @xmath29 , i.e. the force given by , @xmath30{cc}f_1 & \text{for } t \le 0\\ 0 & \text{for } t > 0 . \end{array } \right .
\label{eq : step_force}\ ] ] the applied force should be conjugate to the variable for which the fluctuations are to be calculated . for the cantilever system under consideration here , to determine the fluctuations of the tip displacements , the force is applied to the tip of the left cantilever . the step force couples with the cantilever deflection @xmath23 and the full hamiltonian of the system @xmath31
is given by , @xmath32{cc}h_0 + \delta h & \text{for } t \le 0\\
h_0 & \text{for } t > 0 \end{array } \right .
\label{eq : hamiltonian}\ ] ] where @xmath33 and @xmath34 is the unperturbed hamiltonian . for small perturbations , @xmath20 ,
the quantity @xmath35 will also be small , and the results are greatly simplified . the nonequilibrium ensemble average @xmath21 is given by , @xmath36 where @xmath37 , @xmath38 is boltzmann s constant , and @xmath8 is the absolute temperature . in eq .
( [ eq : noneq_ave ] ) the notation @xmath23 is used to represent the value of @xmath39 evaluated at the phase space coordinates that have evolved from the values @xmath40 and @xmath41 at time @xmath29  @xcite . in the limit of linear perturbations , @xmath42 ,
this simplifies to , @xmath43 where an equilibrium ensemble average is given by , @xmath44 if we now assume that @xmath39 and @xmath45 have the equilibrium average @xmath46 subtracted and recall that @xmath47 , we have our desired result , @xmath48 which can be rearranged to yield , @xmath49 the analysis is similar for the cross correlations of the deflections for two cantilevers in an array , @xmath50 where @xmath51 is the displacement of tip 2 arising from the step force @xmath20 applied to tip 1 .    it is interesting to highlight that the cantilever motions are indeed correlated as indicated by eq .
( [ eq : cross ] ) .
the correlated fluctuations appear contrary to the naive idea that random molecular impacts upon the individual cantilevers should only lead to uncorrelated motion .
additionally , it should be emphasized that the correlated motion of two ( or more ) cantilevers is no more difficult a calculation than determining the autocorrelation of a single cantilever .    on the left hand side of eqs .
( [ eq : auto ] ) and  ( [ eq : cross ] ) are stochastic quantities and on the right hand side are deterministic quantities .
this permits a significant reduction in effort in the calculations of the auto- and cross - correlations of the cantilever deflections .
the macroscopic deflections @xmath21 and @xmath22 can be found from deterministic theory , model equations , or numerical simulations and the stochastic quantities are found simply through application of eqs .
( [ eq : auto ] ) and  ( [ eq : cross ] ) .
the deterministic problem then reduces to the fluid - solid interaction problem given by the removal of a step force on a cantilever . for the step force eq .
( [ eq : step_force ] ) , for @xmath52 , @xmath21 has a finite deflection . after the removal of the step force the cantilever returns to its equilibrium position of @xmath53 . the precise manner in which the cantilever returns to equilibrium
is dominated by the dissipation in the fluid .
an adjacent cantilever will exhibit time dependent deflections given by @xmath22 . in this case
both the initial and final equilibrium deflections are @xmath54 .
the motion of the second cantilever is a result of the fluid motion caused by the first cantilever as it returns to equilibrium .
the spectral properties of the auto- and cross - correlations are determined by taking the cosine fourier transform of eqs .
( [ eq : auto ] ) and  ( [ eq : cross ] ) with appropriate factors ( see eq .
( [ eq : spectral_density_definition ] ) ) .
this yields the noise spectra , @xmath55 and @xmath56 , given by , @xmath57 where @xmath58 is the angular frequency .
the noise spectra are precisely the measured quantity in experiment . in the above expression
we have defined the spectral density of the random process @xmath59 to be , @xmath60 e^{i \omega t } dt \right|^2 , \label{eq : spectral_density_definition}\ ] ] where @xmath61 is the average of @xmath59 over time @xmath8 .    in summary , the explicit steps necessary to determine the stochastic dynamics of a cantilever array are :    1 .
perform the _ deterministic _ calculation . given
some arrangement of cantilevers in a fluid choose one cantilever and apply a step force some time in the distant past .
this will eventually result in a finite steady deflection of this cantilever , whereas the other cantilevers will have zero deflection .
2 .   remove the step force and calculate the deflections of the cantilevers as a function of time as they return to equilibrium , for two cantilevers this is @xmath21 and @xmath22 ( see for example fig .
[ fig : beam_fit ] for the case of one cantilever ) .
this deterministic calculation can be done analytically or numerically depending upon the complexity of the particular situation .
3 .   calculate the _
stochastic _ quantities . using @xmath21 and
@xmath22 calculate the auto- and cross - correlations of the fluctuations from eqs .
( [ eq : auto ] ) and  ( [ eq : cross ] ) .
in essence , the cantilever deflections are merely rescaled to produce the stochastic response ( this is illustrated by the two ordinate axes in fig .
[ fig : beam_fit ] ) .
4 .   use the auto- and cross - correlations to calculate the noise spectra given by eqs .
( [ eq : noise1 ] ) and  ( [ eq : noise2 ] ) .
the equipartition theorem permits the calculation of the mean square mode displacements of micron and nanoscale cantilevers in vacuum  @xcite .
sader  @xcite further developed this approach to include the damping forces of a surrounding viscous fluid and solved the equations of elasticity in the thin beam limit coupled together with the navier - stokes equations .
he used the approximation that the effect of the fluid on each element of cantilever is the same for an element of a long and slender cantilever moving with the same speed . in the limit of an infinitely long cylinder the flow field around the cantilever
can be assumed to be the flow field around an infinite two - dimensional cylinder to a good approximation  @xcite .
in addition , because the forces on each element of cantilever due to the fluid depend only on the velocity of that element , and not otherwise on its position in the cantilever , the mode structure of the damped cantilever is unchanged from the undamped limit .
we follow this approach here and use the infinite cylinder approximation for the flow field .
however we do not assume that the fluid noise is white as was done in  @xcite but use the correct spectral density of the brownian force given by @xmath62 in eq .
( [ eq : fdt ] ) .
the frequency dependence of the brownian noise @xmath62 is determined by the fluid damping @xmath63 ( illustrated in fig .
[ fig : mass_damping ] and discussed further below ) .    to connect with the previous work we note that the susceptibility @xmath64 connecting the tip displacement to an applied force is given by the cantilever response to a unit impulse of force , which is the derivative of the step force response given by eq .
( [ eq : auto ] ) so that @xmath65 the fourier transform of @xmath64 is the response to a deterministic sinusoidal force and the imaginary part @xmath66 , where im\ { } indicates the imaginary component is , @xmath67 where we have adopted the fourier transform convention given by , @xmath68 and @xmath69 . integrating eq .
( [ eq : dissipation ] ) by parts and using the definition of the noise spectrum eq .
( [ eq : noise1 ] ) yields , @xmath70 the susceptibility @xmath64 can be determined by solving for the deterministic cantilever response to a force impulse .    in the following we present results only for the fundamental mode ; however higher modes can be included if desired .
the equation governing the cantilever motion is then , @xmath71 where the amplitude of the mode displacement is characterized by the tip displacement @xmath72 , @xmath73 is the effective mass of the cantilever in vacuum , @xmath74 is the force acting on the cantilever due to the fluid , and @xmath75 is the force impulse ( @xmath76 is the dirac delta ) .
the effective mass of the cantilever @xmath73 is chosen to yield the same kinetic energy as in the cantilever mode , and is related to the cantilever mass by , @xmath77 where @xmath78 is the actual cantilever mass and for the fundamental mode of oscillation for a beam @xmath79 . taking the fourier transform of this equation yields , @xmath80
the force from the fluid can be written in the form , @xmath81 where @xmath82 here @xmath83 is the effective mass of a fluid cylinder of radius @xmath13 where @xmath84 is the fluid density , @xmath85 is the cantilever width , and @xmath7 is the cantilever length . again
the prefactor @xmath86 is to take into account the variation of the fluid force along the cantilever due to the varying velocity given by the mode structure .
the fluid loading and damping are captured by the hydrodynamic function @xmath87 given by , @xmath88 where @xmath89 and @xmath90 are bessel functions  @xcite . in this definition
the frequency dependence on the right hand side appears through the frequency dependent strouhal number @xmath15 .
the cantilever is loaded by the fluid which can be characterized by an effective mass , @xmath91 , larger than @xmath73 that takes into account the fluid mass that is also being moved .
the fluid also damps the motion of the cantilever which can be expressed as an effective damping @xmath92 .
relations for @xmath91 and @xmath92 can be found by expanding @xmath87 into its real and imaginary parts in eq .
( [ eq : eom_fourier_deterministic ] ) and rearranging such that , @xmath93 to give , @xmath94 and , @xmath95 where @xmath96 and @xmath97 are the real and imaginary parts of @xmath98 , respectively .
@xmath99 is the mass loading parameter , which is the ratio of the mass of a cylinder of fluid with radius @xmath13 to the actual mass of the cantilever , and is given by , @xmath100 where @xmath18 is the density of the cantilever and @xmath101 is the cantilever thickness .
it is evident from eqs .
( [ eq : fluid_loaded_mass ] ) and  ( [ eq : damping ] ) that both the fluid loaded mass of the cantilever and the fluidic damping are functions of frequency .
the ratio of the mass of the fluid loaded cantilever to the effective mass of the cantilever in vacuum , @xmath73 , as a function of frequency is given by @xmath102 .
the factor @xmath99 is a constant for any particular choice of cantilever and fluid combination .
the frequency dependence of the added mass is then given by @xmath103 and is illustrated by the solid line in fig .
[ fig : mass_damping ] and uses the left ordinate axis .
the added mass increases rapidly as the frequency is decreased . over the range of four orders of magnitude in frequency the added mass
is seen to change by a factor of approximately 25 .
the frequency dependence of the fluid damping is given by @xmath104 and is shown by the dashed line in fig .
[ fig : mass_damping ] using the right ordinate axis .
there is a weaker frequency dependence for the fluid damping when compared to the mass loading .
the fluid damping decreases as the frequency decreases and over 4 orders of magnitude in frequency the damping changes by a factor of approximately 7 .     illustrates the frequency dependence of the fluid loaded mass .
( dashed line ) @xmath105 illustrates the frequency dependence of the fluid damping.,width=288 ]    solving for the cantilever response @xmath106 and taking the imaginary part yields , @xmath107 using this result in eq .
( [ eq : gxf_susceptibility ] ) and rearranging gives the desired result for the spectral density of the stochastic fluctuations in cantilever displacement , @xmath108 \nonumber},\end{aligned}\ ] ] where @xmath109 is a nondimensional reduced frequency , and @xmath110 is the strouhal number evaluated at the resonant frequency in vacuum given by , @xmath111 we emphasize that eq .
( [ eq : gxf_noise ] ) is _ not _ the same as in previous work  @xcite which failed to include the frequency dependent of the brownian force .    to understand the difference from the previous work we now connect eq .
( [ eq : fdt_displacements ] ) with the calculation in terms of a fluctuating force with spectral density @xmath62
. the fluctuating displacement can be described in terms of the response to this force through the susceptibility @xmath112 inserting eq .
( [ eq : gxf_susceptibility ] ) into this expression yields , @xmath113 it is useful at this point to introduce the impedence @xmath114 where @xmath115 is the velocity of the cantilever tip @xmath116 , it is then straight forward to show , @xmath117 inserting this into eq .
( [ eq : gf_susceptibility ] ) and rearranging yields the desired expression of the fluctuation dissipation theorem , @xmath118 where @xmath119 is the resistance or dissipation and re\ { } indicates the real part . for the case of oscillating cantilevers in fluid @xmath92 is the effective fluid damping .
it is evident from the frequency dependence of the damping that the fluctuating force described by eq .
( [ eq : fdt ] ) is not white noise as considered previously  @xcite ( see fig .  [
fig : mass_damping ] for the variation of @xmath92 with frequency ) .    the spectral density of fluctuations in cantilever displacement @xmath55 can also be determined by directly solving the governing stochastic equation of motion , @xmath120 where @xmath121 is the random force due to brownian motion .
taking the fourier transform of this equation yields , @xmath122 solving for the magnitude of the cantilever response and using eq .
( [ eq : fdt ] ) for the spectral density of the brownian force again yields eq .  ( [ eq : gxf_noise ] ) .
note that integrating this result automatically leads to the equipartition result , @xmath123 we would like to point out that it is not possible to start with the equipartition result given by eq .
( [ eq : equipartition ] ) and to then calculate the spectral properties of the fluctuations , one must first start with the fluctuation - dissipation result given by eq .
( [ eq : fdt ] ) as done here .
the noise spectrum of the cantilever in fluid may be used to calibrate an afm . a simple way to do
this is to extract the effective spring constant from the the frequency @xmath124 giving the maximum of the noise intensity . in previous work
@xcite the fluctuations of the cantilever have been estimated without the frequency dependence of the numerator in eq .
( [ eq : gxf_noise ] ) , leading to some inaccuracy in the estimate of the spring constant .
once @xmath124 is known , the quality factor @xmath125 can be estimated from , @xmath126 where eqs .
( [ eq : fluid_loaded_mass ] ) and  ( [ eq : damping ] ) have been used for the mass and damping .
this expression for @xmath125 is only approximate again because the frequency response of the cantilever is not precisely like that of simple harmonic oscillator , since both @xmath91 and @xmath92 are frequency dependent .    the error incurred by neglecting the frequency dependence in the numerator in eq .
( [ eq : gxf_noise ] ) in fitting the spectrum of a brownian driven cantilever is illustrated in fig .
[ fig : mass_loading_noise_compare ] as a function of @xmath110 and @xmath99 . in fig .
[ fig : mass_loading_noise_compare ] @xmath127 is the frequency at the maximum of the noise spectrum when the frequency dependence has not been included ( specifically , in eq .
( [ eq : gxf_driven ] ) @xmath62 has been assumed constant ) and @xmath124 is the correct value of the frequency using the fluctuation - dissipation theorem as discussed in section  [ section : analytical_theory ] .
the results indicate that as the strouhal number increases the error in this assumption becomes smaller . in particular , when the frequency dependence is neglected the approximate theory always underpredicts @xmath124 .
this can be understood in light of fig .
[ fig : mass_damping ] where it is clear that the magnitude of the damping decreases as the frequency of oscillation is reduced .
when this decrease is not included the predicted value of @xmath124 will be unnecessarily reduced .
although the fluid induced damping slowly decreases as the oscillation frequency is reduced it is important to note that the fluid loaded mass increases rapidly and it is this interplay which results in the small @xmath125 associated with oscillating micron and nanoscale cantilevers .    .
@xmath127 is the predicted frequency at the maximum of the noise spectrum when the frequency dependence of the brownian force is neglected as done in previous work  @xcite and @xmath124 is the correct value given by the results of section  [ section : analytical_theory ] .
results are presented as a function of mass loading @xmath99 and the strouhal number @xmath110 based upon the cantilever s resonant frequency in vacuum .
the curves are for @xmath128 where @xmath129 are labelled and the remaining curves are in sequential order.,width=312 ]    figure  [ fig : mass_loading_noise ] summarizes the theoretical results and allows for easy determination of the stochastic dynamics of a single cantilever of arbitrary geometry placed in an arbitrary viscous fluid . given a particular cantilever and fluid combination the procedure is :    1 .   determine the cantilever spring constant @xmath130 , resonant frequency in vacuum , @xmath131 , and characteristic length scale .
these can be experimental measurements or theoretical calculations .
the characteristic length is a measure of the characteristic half - width of the cantilever .
2 .   determine @xmath110 from eq .
( [ eq : r0 ] ) and @xmath99 from eq .
( [ eq : t0 ] ) .
3 .   use @xmath110 and @xmath99 to determine @xmath132 from fig .
[ fig : mass_loading_noise](a ) and @xmath125 from fig .
[ fig : mass_loading_noise](b ) .     of oscillation for a cantilever undergoing stochastic oscillations in a viscous fluid as a function of mass loading @xmath99 and the strouhal number @xmath110 based upon the cantilever s resonant frequency in vacuum .
panel  ( b ) , the quality factor @xmath125 as a function of @xmath110 and @xmath99 as given by eq .
( [ eq : qapprox ] ) . in each panel
the largest value of @xmath110 is labelled and each successive curve represents @xmath133 , i.e. in panel  ( a ) the last curve is for @xmath134 and in panel  ( b ) the last curve is for @xmath135.,title="fig:",width=312 ]   of oscillation for a cantilever undergoing stochastic oscillations in a viscous fluid as a function of mass loading @xmath99 and the strouhal number @xmath110 based upon the cantilever s resonant frequency in vacuum .
panel  ( b ) , the quality factor @xmath125 as a function of @xmath110 and @xmath99 as given by eq .
( [ eq : qapprox ] ) . in each panel
the largest value of @xmath110 is labelled and each successive curve represents @xmath133 , i.e. in panel  ( a ) the last curve is for @xmath134 and in panel  ( b ) the last curve is for @xmath135.,title="fig:",width=312 ]    in the case that the cantilever is driven externally the response can be found in a similar manner .
if the cantilever is driven by a force @xmath136 the equation of motion becomes , @xmath137 if the driving force is @xmath138 where @xmath139 is a constant and @xmath140 is the driving frequency the amplitude of the response is given by , @xmath141 } \nonumber,\end{aligned}\ ] ] where @xmath142 is the reduced driving frequency .
we emphasize that eqs .
( [ eq : eom_driven ] ) and  ( [ eq : gxf_driven ] ) are for externally driven cantilevers and neglect brownian noise .
( [ eq : gxf_driven ] ) and  ( [ eq : gxf_noise ] ) are similar in that they share a common denominator , whereas the additional frequency dependence in the numerator of eq .
( [ eq : gxf_noise ] ) is from the frequency dependence of the brownian noise .
the stochastic dynamics of a micron scale cantilever , similar to what is commonly used in atomic force microscopy , are now quantified ( see fig .  [
fig : beam ] ) . the beam is chosen from chon and sader  @xcite where both theoretical and experimental values are present for comparison .
the beam properties are summarized in table  [ table : beam_geometry ] ,    , width @xmath85 , and height @xmath101.,width=240 ]    [ h ]    [ c]l@l@l@ l@l@l@l @xmath7 & @xmath85 & @xmath101 & @xmath143 & @xmath110 & @xmath99 +   + 197@xmath144 m & 29@xmath144 m & 2@xmath144 m & 71.56khz & 109.7 & 4.89 +    using classical beam theory  @xcite the equivalent spring constant @xmath130 is predicted to be , @xmath145 which yields a value of @xmath146n / m and the resonant frequency of the cantilever in vacuum is @xmath147 which yields @xmath148khz where @xmath149 . using eq .
( [ eq : r0 ] ) the strouhal number in vacuum is @xmath150 and using eq .
( [ eq : t0 ] ) the mass loading factor is @xmath151 .
[ h ]    l@l@l@l@l@l@l & @xmath125 & @xmath152 & @xmath153 & @xmath92 ( kg / s ) & @xmath154 +   + ( 1 ) & 3.24 & 8.16 & 0.34 & @xmath155 & 37.3 + ( 2 ) & 2.93 & 8.05 & 0.35 & @xmath156 & 38.7 +    we now use the thermodynamic approach to calculate the stochastic dynamics of the cantilever .
as previously stated this can be accomplished in a straightforward manner by determining the deterministic response of the cantilever to the removal of a step force applied to the tip .
we do this using finite element numerical simulations of the full three dimensional , time - dependent , fluid - solid interaction problem ( algorithm discussed elsewhere  ) .
the simulation is initiated with the removal of a step force applied to the tip of the cantilever and the deterministic dynamics of the beam @xmath21 are shown by the solid line in fig .
[ fig : beam_fit ] using the right ordinate axis .    in order to calculate commonly used diagnostics , such as @xmath124 and @xmath125 ,
the deterministic cantilever deflection @xmath21 was fit to the deflections of an underdamped simple harmonic oscillator ( valid for @xmath157 ) given by the equation of motion , @xmath158 where @xmath159 and @xmath160 .
the solution is , @xmath161 where , @xmath162 and @xmath163 .
using a nonlinear least squares curve fit algorithm  @xcite the numerical results are fit to the model to yield the values of @xmath125 , @xmath124 , @xmath92 , and @xmath91 shown in table  [ table : beam_fluid_sader ] .
the curve fit is nearly indistinguishable from the numerical simulation results for @xmath21 and is shown by the dashed line in fig .
[ fig : beam_fit ] .    using @xmath21 in eq .
( [ eq : auto ] ) yields the autocorrelation of the equilibrium deflections , @xmath164 , which are shown in fig .
[ fig : beam_fit ] using the left ordinate axis .
the noise spectrum is given by eq .
( [ eq : noise1 ] ) and is shown in comparison with predictions of the infinite - beam approximation given by eq .
( [ eq : gxf_noise ] ) in fig .
[ fig : beam_noise ] .
the numerical simulations contain all of the oscillation modes as indicated by the presence of the second mode in the figure . in the infinite beam approximation
the different modes behave independently and the calculation could be extended to include higher modes if desired . in table
[ table : beam_fluid_sader ] the analytical results of section  [ section : analytical_theory ] are compared with results using the full thermodynamic approach and deterministic finite element numerical simulations . in general , from the table and the figures , we see that the differences between the simple model and the full calculations are small for this shape cantilever . the experimentally measured value of the resonant frequency in water is @xmath165  @xcite which agrees well with both the analytical results and the numerical predictions using the thermodynamic approach .
however , table  [ table : beam_fluid_sader ] suggests that the analytical results based upon the oscillating infinite cylinder model under predict the amount of damping @xmath92 .
the noise spectra show differences that are small , but may be significant for quantitative applications such as calibration .    in all of the deterministic finite element numerical simulations
performed we were careful to ensure that the bounding no - slip surfaces of the computational domain did not affect the results .
comparisons with numerical simulations performed with larger domains did not significantly alter the results .
for a micron scale cantilever in water .
( solid line ) deterministic finite element numerical simulations using the thermodynamic approach .
( dashed line ) approximate analytical theory for the fundamental mode only.,width=288 ]
as the cantilever dimensions become smaller the effective cantilever spring constant decreases while the resonant frequency increases .
this favorable combination potentially provides access to the biologically important parameter regime characterized by 10 s of piconewtons with microsecond scale time resolution , a range that is difficult to reach using other methods .
this has led to the development of nanoscale cantilevers  @xcite .    in what follows
we quantify the stochastic dynamics of two adjacent nanoscale cantilevers immersed in water .
a schematic of the nanoscale cantilever under consideration here is shown in fig .  [
fig : paddle_config ]  @xcite .
this is an experimentally motivated cantilever whose paddle shaped geometry decreases the effective spring constant while localizing the strain for piezoresistive measurement .
the configuration of the cantilever array is shown in fig .
[ fig : adjacent ] with two opposing cantilevers separated by a distance @xmath166 .
this configuration was chosen because of its experimental accessibility as well as its potential use in single molecule measurements with the tethering of a target biomolecule between the two cantilevers . in the work presented here
we build a baseline understanding of the cantilever dynamics in the absence of target biomolecules .
further consideration of the dynamics caused by a tethered biomolecule is beyond the scope of the present paper .
the dimensions of the nanoscale cantilevers investigated here are summarized in table  [ table : paddle - geometry ] .
as discussed previously , the physical properties describing the dynamics of the cantilever in vacuum are important parameters in determining the stochastic dynamics of the cantilevers .
these properties are summarized in table  [ table : paddle_vacuum ] where the values have been determined from finite element simulations of the elastic cantilever structure in the absence of the surrounding fluid .
it is important to emphasize that although these calculations can sometimes be performed analytically , this is not the case for many geometries of interest .
however , even for cantilevers with very complex geometries , the vacuum based finite element calculations are straight forward .
[ h ]    [ c]l@l@l@l@l@l @xmath7 & @xmath85 & @xmath101 & @xmath167 & @xmath168 +   + 3@xmath144 m & 100 nm & 30 nm & @xmath169 m & 33 nm +    [ h ]    [ c]l@l@l@l@l@l@l @xmath130 & @xmath131 & @xmath110 & @xmath99 +   + 8.7 mn / m & @xmath170 rads / s & 0.11 & 1.12 +    , width @xmath85 , and height @xmath101.,width=264 ]    the stochastic dynamics of the cantilever array are determined using the thermodynamic approach discussed in section  [ section : thermodynamic_approach ] with deterministic finite element numerical simulations .
a series of numerical simulations have been performed to determine the variation in the cantilever dynamics as a function of the separation distance between the two cantilevers in the array .
in particular we have explored @xmath171 where @xmath166 is the cantilever separation and @xmath172 nm is the cantilever thickness .
the numerical simulations are initiated by the removal of a step force applied to the tip of one of the cantilevers ( in fig .
[ fig : adjacent ] this is the cantilever on the left ) while the other cantilever is initially undeflected and at mechanical equilibrium .
the deterministic dynamics of the nanoscale cantilever for which the step force has been removed , @xmath21 , is shown by the solid line in fig .
[ fig : one_paddle_auto1 ] using the right ordinate axis .
the deflection @xmath45(t ) was not significantly affected by the presence of the second cantilever for the separations considered here . as with the micron scale cantilever
, we again fit the result for @xmath21 from the finite element numerical simulation to the solution of the simple harmonic oscillator equation given by eq .
( [ eq : sho_fit ] ) . for the nanoscale cantilever
the dynamics are overdamped ( @xmath173 and the deterministic deflection predicted by the simple harmonic oscillator approximation is , @xmath174 where , @xmath175 the agreement with the curve fit is very good and @xmath124 and @xmath125 are shown in table  [ table : paddle_fluid ] . also shown in table  [ table : paddle_fluid ] are the theoretical predictions based upon the infinite cylinder approximation . for the analytical calculations of the damping @xmath92 from eq .
( [ eq : damping ] ) the fundamental mode of the nanoscale cantilever has been modelled as that of a hinge where all of the bending occurs in the short legs near the base .
this is guided by numerical calculations showing the strain localized in the short legs . in this case
we find the mass factor is @xmath176 .
( note that @xmath177 is not needed to determine the other quantities in table  [ table : paddle_fluid ] ) .
it is clear from table  [ table : paddle_fluid ] that the infinite cylinder approximation is no longer valid .
the analytical prediction for the added mass is an order of magnitude too large .
this suggests that three - dimensional flow effects become significant for the shape of cantilever under consideration here .
for example there would be flow around the tip of the cantilever , as well as the flow through the open region near the base of the cantilever .
in addition , the predicted frequency shift is an order of magnitude too small .
[ h ]    l@l@l@l@l@l@l & @xmath125 & @xmath152 & @xmath153 & @xmath92 ( kg / s ) & @xmath154 +   + ( 1 ) & 0.265 & 106.3 & 0.032 & @xmath178 & 0.0036 + ( 2 ) & 0.327 & 18.16 & 0.236 & @xmath179 & 0.026 +     ( left ordinate axis ) for a single nanoscale cantilever in fluid : ( solid line ) results from finite element numerical simulations using the thermodynamic approach ; ( dashed line ) results based upon the approximate analytical theory of section  [ section : analytical_theory].,width=288 ]    inserting @xmath21 into eq .
( [ eq : auto ] ) yields the autocorrelation of the equilibrium fluctuations in cantilever displacement , @xmath164 .
these are also shown by the solid line in fig .
[ fig : one_paddle_auto1 ] using the ordinate axis on the left .
the noise spectrum , @xmath55 , is found by inserting the autocorrelation in eq .
( [ eq : noise1 ] ) and the result is shown by the solid line in fig .
[ fig : one_paddle_noise_comparison ] . fig .
[ fig : one_paddle_noise_comparison ] also presents a comparison of the noise spectrum @xmath55 given from the theoretical predictions of eq .
( [ eq : gxf_noise ] ) and the results from the numerical simulations using the thermodynamic approach .
the noise spectra are different for small frequencies , in particular the model result has a peak at finite frequency not seen in the full calculation .
the inverse cosine transform of the noise spectrum from the infinite - cylinder model yields the deflection of the cantilever as a function of time shown by the dashed line in fig .
[ fig : one_paddle_auto1 ] using the right ordinate axis . the model prediction for the cantilever
deflection exhibits negative values that are not seen in the simulations suggesting that the response of the overdamped cantilever is not precisely that of a simple harmonic oscillator .     for a single nanoscale cantilever immersed in water .
( solid line ) results from deterministic finite element simulations using the thermodynamic approach described in section  [ section : thermodynamic_approach ] .
( dashed line ) results from approximate analytical theory described in section  [ section : analytical_theory].,width=288 ]    we now turn to the deterministic deflection of the second cantilever , @xmath22 , after removing the force on the first cantilever , evaluated using the full fluid - elasticity simulations .
this is shown in fig .
[ fig : two_paddles_cross_all ] using the right ordinate axis . for close separations
@xmath180 the second cantilever exhibits negative deflections for all time .
however , as the cantilever separation increases the second cantilever initially exhibits a positive deflection followed by a negative deflection .
the flow field around an oscillating object will have both potential and non - potential components .
for an incompressible fluid the potential component is instantaneous whereas the non - potential component diffuses with a diffusion coefficient given by the kinematic viscosity @xmath181 .
although the deterministic flow field around simple oscillating objects is well known the fluid coupled motion of multiple elastic objects is not well understood .
the fluid dynamics resulting from the motion of two adjacent cantilevers will be discussed in a forthcoming article  @xcite .
using eq .
( [ eq : cross ] ) the stochastic correlations between the two cantilevers can be determined from @xmath22 .
the cross - correlations of the equilibrium displacement fluctuations @xmath182 are shown in fig .
[ fig : two_paddles_cross_all ] using the ordinate axis on the left . for small cantilever separations
@xmath180 the cross - correlations are anticorrelated for all time .
this is in agreement with experimental measurements of the fluid correlations of two closely spaced micron scale beads in water  @xcite .
however , as the cantilever separation is increased the cross - correlations changes : for short times the cross - correlations are positively correlated whereas for larger times they are anti - correlated .
the noise spectra @xmath56 are found from @xmath182 by using eq .
( [ eq : noise2 ] ) . figure  [ fig : two_paddles_noise_all ] shows the variation in the noise spectra as a function of cantilever separation .
the noise spectra contain both positive and negative values and for each cantilever separation there is a frequency at which correlated noise vanishes , @xmath183 .
this null point could be exploited by a measurement scheme to minimize the correlated noise .
for @xmath172 nm found using deterministic finite element numerical simulations with the thermodynamic approach discussed in section  [ section : thermodynamic_approach ] .
the curve for @xmath184 is labelled and the others follow in sequential order .
the right ordinate axis illustrates the deterministic cantilever deflection as a function of time @xmath22 scaled by the deflection of the adjacent cantilever , @xmath185).,width=288 ]     for two nanoscale cantilevers separated by distance @xmath186 for @xmath172 nm found using deterministic finite element numerical simulations with the thermodynamic approach discussed in section  [ section : thermodynamic_approach ] .
the curve for @xmath184 is labelled and the others follow in sequential order.,width=288 ]    .
panel  ( a ) , the maximum magnitude of the force induced by fluid - coupled correlations in cantilever displacement @xmath187 .
the data is empirically fit with an exponential given by @xmath188 ( in units of pn ) .
the magnitude of the force felt by a single cantilever due to brownian noise , @xmath189 pn , is nearly 6 times larger than the correlated noise at close separations ( @xmath180 ) .
the right ordinate axis scales @xmath187 by @xmath190 .
panel  ( b ) , characteristic time scales , the circles represent the time at which the maximum value of the cross correlation occurs , the squares represent the time at which the fluid induced correlations vanish .
the ordinate axis on the right hand side illustrates the time scales when scaled by the period of oscillation of the cantilever in vacuum @xmath191 .
the solid lines represent a linear curve fit to the data.,title="fig:",width=264 ] .
panel  ( a ) , the maximum magnitude of the force induced by fluid - coupled correlations in cantilever displacement @xmath187 .
the data is empirically fit with an exponential given by @xmath188 ( in units of pn ) .
the magnitude of the force felt by a single cantilever due to brownian noise , @xmath189 pn , is nearly 6 times larger than the correlated noise at close separations ( @xmath180 ) .
the right ordinate axis scales @xmath187 by @xmath190 .
panel  ( b ) , characteristic time scales , the circles represent the time at which the maximum value of the cross correlation occurs , the squares represent the time at which the fluid induced correlations vanish . the ordinate axis on the right hand side illustrates the time scales when scaled by the period of oscillation of the cantilever in vacuum @xmath191 .
the solid lines represent a linear curve fit to the data.,title="fig:",width=264 ]    from these results it is possible to characterize the force sensitivity and time resolution of a correlation measurement technique using an array of closely spaced nanoscale cantilevers .
an estimate of the force sensitivity can be found using the auto and cross correlation functions as , @xmath192 @xmath190 represents the approximate magnitude of the stochastic brownian force acting on a single cantilever and @xmath187 the force induced by fluid correlations between two cantilevers . for the case of a single nanoscale cantilever
the maximum value occurs at @xmath29 and , using the equipartition theorem this yields a value of @xmath193pn .    in a cross - correlation measurement between two cantilevers the brownian noise felt by the two individual cantilevers is uncorrelated and does not contribute .
this leaves only the correlations due to the viscous fluidic coupling and @xmath194 represents the approximate magnitude of the force due to this hydrodynamic coupling .
for example , from fig .
[ fig : two_paddles_cross_all ] for the case with the closest separation @xmath184 , the magnitude of the maximum value of the cross - correlation is @xmath195nm@xmath196 which yields a force sensitivity @xmath197pn . therefore @xmath198 at @xmath184 , indicating a 6-fold reduction in thermal noise over a single cantilever measurement .    the variation in @xmath187 as a function of cantilever separation is shown in fig .
[ fig : force_and_time](a ) .
the data is fit with an exponential curve , given by @xmath199 ( in units of pn ) , shown in fig .
[ fig : force_and_time](a ) as the solid line .
for very small cantilever separations @xmath180 the noise reduction possible through the use of a correlated measurement technique is approximately 6-fold .
as expected , the correlated noise decreases as the cantilever separation is increased : for a separation of @xmath200 there is a 12-fold noise reduction . the ordinate axis on the right hand side of fig .
[ fig : force_and_time](a ) shows @xmath201 to show this performance improvement directly .
the characteristic time scales of the fluid correlated motion as a function of cantilever separation are illustrated in fig .
[ fig : force_and_time](b ) . in fig .
[ fig : force_and_time](b ) circles represent the time at which the maximum value in the cross correlation occur and squares represent the time at which the fluid induced correlations vanish ( i.e. , the zero - crossing in fig .
[ fig : two_paddles_cross_all ] ) . as expected from the finite rate at which momentum diffuses away from an oscillating cantilever , both of these time scales increase with cantilever separation . using the ordinate axis on the right in fig .
[ fig : force_and_time](b ) the ratio of these time scales to the period of a single oscillation for the cantilever in vacuum @xmath202 are given .
the stochastic dynamics of micron and nanoscale cantilevers can be quantified for the precise conditions of experiment using straightforward deterministic calculations coupled with the fluctuation - dissipation theorem .
this can be done for complex cantilever geometries including wall effects , and the stochastic correlated behavior of closely spaced cantilevers for which theoretical predictions currently are not available .
cantilever geometries that are long and thin are well described by approximate analytical theory based upon the oscillation of an infinite cylinder in fluid .
we have used the fluctuation - dissipation method to correct the previous results of this theory . for the short and wide nanocantilevers
currently being proposed the infinite - cylinder approximation is no longer appropriate .
the thermodynamic approach using the fluctuation - dissipation theorem provides an important means of developing a physical understanding of the stochastic dynamics of closely spaced micron and nanoscale objects that will be important as micron and nanotechnology progress and a quantitative design tool for experiment .
this research has been partially supported by darpa / mto simbiosys under grant f49620 - 02 - 1 - 0085 and an aspires grant from virginia tech .
this work has been carried out in collaboration with the caltech bionems effort ( m. l. roukes , pi ) and we gratefully acknowledge extensive interactions with this team .


two observations of with the rosat pspc resulted in the detection of x - ray emission from 25 members of this old open cluster ( belloni et al .
1993 , 1998 ) . the x - ray emission of many of these sources is readily understood .
for example , the x - ray emission originates in deep , hot atmospheric layers in a hot white dwarf ; is due to mass transfer in a cataclysmic variable ; and is caused by magnetic activity in two contact binaries and several rs cvn - type binaries . however , belloni et al .
( 1998 ) point out several x - ray sources in for which the x - ray emission is unexplained .
all but one of these objects are located away from the isochrone formed by the main sequence and the ( sub)giant branch of ( fig .
[ cmd ] ) .    in this paper
we investigate the nature of the x - ray emission of these stars through low- and high - resolution optical spectra .
in particular , we investigate whether the emission could be coronal as a consequence of magnetic activity , by looking for emission cores in the h&k lines .
tidal interaction in a close binary orbit is thought to enhance magnetic activity at the stellar surface by spinning up the stars in the binary .
therefore , we also derive projected rotational velocities with the crosscorrelation method .
finally , we study the h@xmath0 profile as a possible indicator of activity or mass transfer .
the observations and the data reduction are described in sect .  2 , and the analysis of the spectra in sect .
comparison with chromospherically active binaries is made in sect .
a discussion of our results is given in sect .
5 . in the remainder of the introduction
we give brief sketches of the stars studied in this paper ; details on many of them are given by mathieu et al .
the stars are indicated with their number in sanders ( 1977 ) , and are listed in table  [ tab1 ] .
l@l@l@l@l@l@l i d & @xmath1 & @xmath2-@xmath1 & sp.type & @xmath3 & @xmath4 & ctrate + & & & & ( d ) & & ( s@xmath5 ) +   +   + s1040 & 11.52 & 0.87 & g4iii & 42.83 & 0.027(28 ) & 0.0050(6 ) + s1063 & 13.79 & 1.05 & g8iv & 18.39 & 0.217(14 ) & 0.0047(6 ) + s1072 & 11.32 & 0.61 & g3iii  iv & 1495 . & 0.32(7 ) & 0.0013(3 ) + s1082 & 11.25 & 0.42 & f5iv & & & 0.0046(6 ) + s1113 & 13.77 & 1.01 & & 2.823 & 0.031(14 ) & 0.0047(6 ) + s1237 & 10.78 & 0.94 & g8iii & 697.8 & 0.105(15 ) & 0.0010(3 ) + s1242 & 12.72 & 0.68 & & 31.78 & 0.664(18 ) & 0.0007(2 ) +    s1063 and s1113 are two binaries located below the subgiant branch in the colour - magnitude diagram of .
their orbital periods , 18.4 and 2.82 days respectively , are too long for them to be contact binaries ; also they are too far above the main sequence to be binaries of main - sequence stars . in principle , a ( sub)giant can become underluminous when it transfers mass to its companion , as energy is taken from the stellar luminosity to restore hydrostatic equilibrium ( e.g.  kippenhahn &  weigert 1967 ) .
however , mass transfer through roche lobe overflow very rapidly leads to circularization of the binary orbit , whereas s1063 has an eccentricity @xmath6 .
the orbit of s1113 is circular , so mass transfer could be occurring in that system . for the moment ,
the nature of these binaries is not understood . in both , pasquini &  belloni ( 1998 )
observed emission cores in the h&k lines .
s1063 is reported to be photometrically variable with @xmath7 0.10 mag ( rajamohan et al .  1988 ; kaluzny & radczynska 1991 ) , but no period is found . for s1113 , photometric variability with a period of 0.313 days and a total amplitude of 0.6 mag
was claimed by kurochkin ( 1960 ) , but this has not been confirmed by kaluzny & radczynska ( 1991 ) , who find variability with only 0.05 mag .
s1063 is the only star in our sample that shows significantly variable x - ray emission ( between 0.0081 and 0.0047 cts s@xmath5 ; belloni et al .
1998 ) .
s1072 and s1237 are binaries with orbital periods of 1495 and 698 days , and with eccentricities @xmath8 and 0.105 , respectively .
the colour and magnitude of s1072 can not be explained with the pairing of a giant and a blue straggler , since this is not compatible with its @xmath9 photometry ( nissen et al .
1987 ; mathieu &  latham 1986 ) , nor with superposition of three subgiants , since this is excluded by the radial velocity correlations ( mathieu et al .
the absence of the 6708   lithium feature in the spectrum of s1072 indicates that the surface material has undergone mixing ( hobbs &  mathieu 1991 ; pritchett & glaspey 1991 ) .
s1237 could be a binary of a giant and a star at the top of the evolved main sequence ( janes & smith 1984 ) ; high - resolution spectroscopy should be able to detect the main - sequence star in that case ( mathieu et al .  1990 ) .
the wide orbits and significant eccentricities appear to exclude both mass transfer and tidal interaction as explanations for the x - ray emission .
s1242 has the largest eccentricity of the binaries in our sample , at @xmath10 in an orbit of 31.8@xmath11days .
its position on the subgiant branch is explained if a subgiant of 1.25@xmath12 has a secondary with @xmath13 ( mathieu et al .
k line emission is reported by pasquini &  belloni ( 1998 ) .
photometric variability with a period of 4.88@xmath11days and amplitude of 0.0025 mag has been found by gilliland et al .
we note that this photometric period corresponds to corotation with the orbit at periastron , which suggests that the x - ray emission may be due to tidal interaction taking place at periastron .
the binary would then be an interesting example of a system in transition from an eccentric to a circular orbit .
indeed , according to the diagnostic diagram of verbunt &  phinney ( 1995 ) a giant of 1.25@xmath12 with a current radius of @xmath14 ( as derived from the location of s1242 in the colour - magnitude diagram ) can not have circularized an orbit of 31.8days .
s1040 is a binary consisting of a giant and a white dwarf .
the progenitor of the white dwarf circularized the orbit during a phase of mass transfer ( verbunt &  phinney 1995 ) ; as a result the mass of the white dwarf is very low ( landsman et al .
the white dwarf is probably too cool , at 16160@xmath11k , to be the x - ray emitter .
indications for magnetic activity are h&k ( pasquini &  belloni 1998 ) and ( @xmath15 2800  , landsman et al .
1997 ) emission lines .
if the x - rays are due to coronal emission of the giant , this must be the consequence of the past evolution of the binary , since the giant is too small for significant tidal interaction to be taking place in the current orbit .
s1082 is a blue straggler .
photometric variability of 0.08 mag within a few hours was observed by simoda ( 1991 ) .
goranskii et al .
( 1992 ) found eclipses with a total amplitude of 0.12 mag and a binary period of 1.07 days ; however , the radial velocities of the star do not show this period , and vary by about 2 km s@xmath5 , far too little for a 1day eclipsing binary ( mathieu et al .
landsman et al .
( 1998 ) detect a significant excess at 1520   with the ultraviolet imaging telescope , and ascribe this to a hot , subluminous secondary .
such a secondary was suggested already by mathys ( 1991 ) on the basis of a broad component in the d and absorption lines .
optical spectra were obtained on february 28/29 , 1996 with the 4.2 m william herschel telescope on la palma , under good weather conditions ( seeing @xmath16 until 4@xmath1730 ut , @xmath18 thereafter ) .
in addition to the x - ray sources in we observed two ordinary member giants of , s1288 and s1402 , for comparison .
furthermore one flux standard and three velocity standards were observed .
the blue high - resolution spectra of s1113 were obtained on april 7/8 , 1998 with the same telescope through a service observation ( seeing 12@xmath19 ) .
a log of the observations is given in table  [ tab11 ] .
llrlrlr i d & & & + & & & + & ut & & ut & & ut & + & & & & & & +   +   +   + s1040 & 02:43&60 & 22:03&600 & 00:46&300 + s1063 & 02:30&180 & 20:51&1200 & 00:11&1200 + & & & 21:15&1200 & 00:33&600 + & & & 21:37&1200 + s1072 & 02:49&60 & 22:24&360 & 01:01&240 + s1082 & 02:39&60 & 20:39&300 & 00:01&360 + & & & 23:24&600 & 01:41&360 + & & & & & 04:11&360 + s1113 & 03:01&180 & 22:37@xmath20 & 900 & 03:36&1200 + & & & 22:58@xmath20 & 900 & 03:58&600 + s1237 & 02:46&60 & 22:16&300 & 00:54&180 + s1242 & 02:55&120 & 22:33&1500 & 01:08&900 +   +   + s1288 & 02:58&60 & 23:12&600 & 01:33&300 + s1402 & 02:52&60 & 23:01&450 & 01:26&240 +   +   + & 02:21&80 & 23:42&600 & 01:57&600 + & 03:13&80 + & 04:43&480 +   +   + & & & 06:04 & 90 & 05:13&45 + & & & 05:55&100 & 05:42&50 + & & & & & 05:34&45 +    all spectra have been reduced using the image reduction and analysis facility ( iraf ) .
low - resolution spectra were taken with the isis double - beam spectrograph ( carter et al .
the blue arm of isis was used with the 300 lines per mm grating and tek - ccd , resulting in a wavelength coverage of 3831 to 5404   and a dispersion of 1.54   per pixel at 4000  .
the red arm , combined with the 316 lines per mm grating and eev - ccd , covered a wavelength region of 5619 to 7135   with a dispersion of 1.40   per pixel at 6500  .
the format of the frames is 1124 @xmath21 200 pixels which includes the under- and overscan regions . for
the object exposures the slit width was set to @xmath22 .
flatfields were made with a tungsten lamp while cuar and cune lamp exposures were taken for the purpose of wavelength calibration .
for the isis - spectra , basic reduction steps have been done within the iraf ccdred - package .
these steps include removing the bias signal making use of the under- and overscan regions and zero frames , trimming the frames to remove the under- and overscan , and flatfielding to correct for small pixel - to - pixel gain variations .
the remaining reduction has been done with iraf specred - package tasks . with the optimal extraction algorithm ( horne 1986 )
the two dimensional images are reduced to one dimensional spectra .
next , the spectra are calibrated in wavelength with the arc frames .
a dispersion solution is found by fitting third ( blue ) and fourth ( red ) order polynomials to the positions on the ccd of the arclamp lines .
the fluxes of the spectra are calibrated with the absolute fluxes of hz44 , tabulated at 50   intervals ( massey et al .
1988 ) , and adopting the standard atmospheric extinction curve for la palma as given by king ( 1985 ) .
the estimated accuracy of the flux calibration is @xmath23% .
high - resolution echelle spectra were taken with  the utrecht echelle spectrograph ( ues , unger et al .
observations were done with a @xmath24 slitwidth .    for the 1996 observations ,
the ues was used in combination with a 1024 @xmath21 1024 pixels tek - ccd , and the 31.6 lines per mm grating ( e31 ) , which resulted in a broad wavelength coverage , but small separation of the echelle orders on the ccd . in this setup ,
the ues resolving power is 49000 per resolution element ( two pixels ) , corresponding to a dispersion of 3 km s@xmath5 per pixel or 0.06   per pixel at 6000  .
the frames were centered on @xmath25   and @xmath26   in order to get a blue ( 3820 to 4920  ) and red ( 4890 to 7940  ) echelle spectrum .
the number of orders recorded on the ccd is 34 in the blue and 45 in the red , each covering @xmath7 45 to 80 increasing for longer wavelengths . towards the red ,
gaps occur between the wavelength coverage of adjacent orders .
exposures of a quartz lamp were taken to make the flatfield corrections .
thar exposures served as wavelength calibration frames .    for the 1998 observations of s1113 ,
a 2048 @xmath21 2048 pixels site - ccd was used .
two spectra were taken with the 79.0 lines per mm ( e79 ) grating ( @xmath27  ) .
the difference between the e79 and the e31 gratings is that e79-spectra have a larger separation of the echelle - orders on the detector , which can improve the determination of the sky - background .
the spectral resolution of the gratings is the same ( the central dispersion in these observations is @xmath28   per pixel ) . the specified wavelength - coverage for this combination of grating and detector is 3546 to 6163  but only the central orders were bright enough to extract spectra ( 3724 to 5998  ) . flatfield and thar exposures
were made for calibration purposes .
the reduction of the ues spectra has been performed using the routines available within the iraf echelle - package .
first , the frames are debiased and the under- and overscan regions removed . after locating the orders on the ccd for both the quartz lamp and the object exposures , we flatfielded the frames .
spectra are extracted with optimal extraction .
the small order separation makes sky subtraction difficult ; however , our targets are bright , and the resulting error is negligible . in the step of wavelength calibration ,
the dispersion solution is derived by fitting third and fourth order polynomials leaving rms - residuals of 0.004   ( red ) and 0.002   ( blue , 0.003   for the 1998-spectra ) . to find absolute fluxes for the k ( @xmath29  ) & h ( @xmath30  ) emission lines ( sect .
3.1 ) , the fluxes of the relevant blue orders of an object have been calibrated with the calibrated isis spectrum of the same object .
continuum normalization of the orders in the red spectra , required for the rotational velocity analysis , is done by fitting third to fifth order polynomials to the wavelength - calibrated spectra .
we study two indicators of magnetic activity .
the direct indicator is emission in the cores of the h&k lines .
another indicator is the rotational speed : rapid ( differential ) rotation and convective motions are thought to generate magnetic fields through a dynamo .      to estimate the amount of flux emitted in the h&k line cores , @xmath31 , we add the fluxes above the h&k absorption profiles as follows .
an upper and a lower limit of the level of the absorption pseudo - continuum is estimated by eye and is marked by a straight line . for s1113 this is illustrated in fig .
[ subca ] .
we obtain a lower and upper limit of the emitted flux by adding the fluxes in each wavelength - bin above these levels .
the value given in table  [ fluxes ] is the average of these two results , the uncertainty is half their difference .
use of higher order fits ( following fernndez - figueroa et al .
1994 ) to the absorption profile gives similar results .
if an emission line is not clearly visible , we obtain an upper limit by estimating the minimal detectable emission flux at the h&k line centers within a 1   wide region ( typical width of the emission lines ) .
six of our sources show h&k line emission ( fig .  [ cahk ] ) .
the profiles of s1113 appear to be double - lined , suggesting that we see activity of both stars ( fig .
[ subca ] ) .
the fluxes given in table  [ fluxes ] are the total fluxes , i.e.  no attempt was made to deblend the emission lines .
no emission is visible in the spectrum of s1082 .
lrrl i d & & & + & & & +   +   + s1040 & @xmath32 & @xmath33 & 3.0(1.0 ) + s1063 & @xmath34 & @xmath35 & 3.9(0.8 ) + s1072 & @xmath36 & @xmath37 & 8.1(1.1 ) + s1082 & @xmath38 & @xmath39 & 5.1(0.7 ) + s1113 & @xmath40 & @xmath41 & 45(6 ) + & & & 12(1 ) + s1237 & @xmath42 & @xmath43 & @xmath44 + s1242 & @xmath45 & @xmath46 & @xmath47 +        in order to derive the projected rotational velocity @xmath48 of our targets , we apply the crosscorrelation technique ( e.g. tonry & davis 1979 ) .
this method computes the correlation between the object spectrum and an appropriately chosen template spectrum as function of relative shift .
the position of the maximum of the crosscorrelation function ( ccf ) provides the value of the radial - velocity difference between object and template .
the width of the peak is indicative for the width of the spectral lines and can therefore be used as a measure of the rotational velocity of the stars .    for rotational velocities not too large
the line profiles may be approximated with gaussians , allowing analytical treatment of the crosscorrelation method .
assuming that the binary spectrum is a shifted , scaled and broadened version of the template spectrum , the broadening can be related to the width of the ccf peak as follows ( gunn et al .
with @xmath49 the dispersion in the template s and @xmath50 the dispersion in the target s spectral lines , @xmath51 the dispersion of the ccf peak and @xmath52 the dispersion of the gaussian that describes the broadening of the object s spectrum with respect to the template , one can write : @xmath53 eq .
[ sigma ] applies to both components in the binary spectrum and their corresponding ccf peaks .
the crosscorrelations are performed with the iraf task fxcor that uses fourier transforms of the spectra to compute the ccf . before performing the crosscorrelation ,
the continuum is subtracted from the normalized spectra .
filtering in the fourier domain is applied to avoid undesirable contributions originating from noise or intrinsically broad lines ( see wyatt 1985 ) .
the templates are chosen from the radial - velocity standards such that their spectral types resemble those of the targets .
the value of @xmath49 for these stars is determined for each order separately by autocorrelation of the template spectrum adopting the same filter used for the crosscorrelations . in this case
@xmath52 is zero and therefore @xmath49 is found directly from the width of the ccf peak : @xmath54 .
the template spectra are correlated with our target stars order by order , where we limit ourselves to orders in the red spectra that do not suffer from strong telluric lines .
most ccf peaks can be fitted well with a gaussian .
as the final value for @xmath48 we give the broadening @xmath52 averaged over the different orders , and for the uncertainty we take the rms of the spread around the average of @xmath52 ( table  [ fluxes ] ) . equating @xmath48 to @xmath52
implicitly assumes that @xmath49 is the width of the lines not related to rotation .
an upper limit to @xmath48 is found from the other extreme in which we assume that the total width of the spectral line @xmath50 follows from rotation .
this creates uncertainties of the order of @xmath52 for s1242 . in the case of s1237
we find that @xmath55 ( @xmath56 ) , i.e. the lines in s1237 are narrower than in the template .
for these two stars , we give an upper limit of @xmath57 in table  [ fluxes ] .
s1113 is the only binary observed whose ccf shows two peaks .
the ccf peaks of both stars overlap in the 1998 observation .
therefore we do not use these spectra in the following analysis .
for the 1996 spectra , the ccf shows two peaks , one of which is broad indicating the presence of a fast rotating star .
both peaks are clearly separated with a center - to - center velocity separation of @xmath7 110 km s@xmath5 .
the lines in the spectrum are smeared out and less pronounced , resulting in a noisy ccf . to improve this
, we combine four sequential orders before crosscorrelating ( being constrained by the maximum number of points fxcor can handle ) . from eq .
17 in gunn et al .
we derive the relative light contribution of both components to the spectrum from the height and dispersion of the crosscorrelation peaks , assuming that the binary stars have the same spectrum as the template ( @xmath58 in eq .
according to this the rapidly rotating star contributes @xmath59%  of the light . note that luminosity ratios derived from cross - correlations are uncertain and
should be confirmed photometrically .
the line profile of the fast rotating star in s1113 is not compatible anymore with a gaussian , therefore we adopt another method to determine its @xmath48 , described in piters et al .
this method uses the property that the fourier - bessel transform of a spectral line that is purely rotationally broadened has a maximum at the position of the projected rotational velocity . in practice , this position is a function of the limits over which the fourier transform is performed . the local maxima of this velocity - versus - cutoff - frequency ( vcf ) function approach @xmath48 within half a percent , the result growing more accurate for maxima at higher frequencies .
this error is negligible when compared to errors arising from noise , other line broadening mechanisms , etc .
( see piters et al .
1996 ) . in our determination of @xmath48 of the primary of s1113
, we have used the first local maximum in the vcf - plot of the transformation of four isolated lines at @xmath15 6265.14 , 6400.15 , 6408.03 and 6411.54  .
the @xmath48 in table  [ fluxes ] is an average of the resulting values ; the @xmath48 of the secondary is found from crosscorrelation .
the application of the fourier - bessel transformation method is limited on the low velocity side by the spectral resolution : the fourier transform can not be performed beyond the nyquist frequency which for slow rotators lies at a frequency that is lower than the cutoff - frequency at which the first maximum occurs in the vcf - plot . for our spectra
this means that the method can not be used for @xmath60 5.2 km s@xmath5 . indeed , for every star for which the crosscorrelation method gives a @xmath48 smaller than this value , the vcf - plot does not reach the first local maximum , except for s1082 for which we find @xmath48 = 9.5(1.6 ) . for s1072 ,
the fourier - bessel transform gives a @xmath48 of 12.7(1.0 ) km s@xmath5 .
spectral lines were selected from those used in groot et al .
( 1996 ) and from the additional lines used for s1113 .
the results of our search for emission cores in the h&k lines are displayed in fig .
[ cahk ] and in table  [ fluxes ] .
the emission lines are strong in s1063 , s1113 and s1040 ; and still detectable in s1242 which indicates chromospheric activity in these stars . in s1072 and s1237 the emission cores are marginal , and on s1082 we can only determine an upper limit .
the ( projected ) rotational velocities of all of our stars are relatively small @xmath61 km s@xmath5 , with the exception of s1113 .    in sect .
4.1 we investigate whether the relations between x - ray emission , strength of the emission cores in the h&k lines , and the rotational velocities of the unusual x - ray emitters in are similar to the relations found for well - known magnetically active stars , the rs cvn binaries . in sect .
4.2 we briefly discuss the behaviour of the h@xmath0 line and spectral lines other than h&k that are indicators of chromospheric activity .
individual systems are discussed in sect .
4.3 .      to investigate whether the x - rays of the stars studied in this paper are related to magnetic activity
, we compare their optical activity indicators and x - ray fluxes with those of a sample of rs cvn binaries . in particular , we select rs cvn binaries for which fluxes of the emission cores in the h&k lines have been determined from high - resolution spectra by fernndez - figueroa et al .
( 1994 ) . to obtain x - ray countrates for these binaries
, we searched the rosat data archive for pspc observations of them .
we then analyzed all these observations , and determined the countrates , in the same bandpass as used in the analysis of , using the standard procedure described in zimmermann et al .
all pointings that we have analyzed actually led to a positive detection of the rs cvn system : even when not the target of the observation , the rs cvn system is usually the brightest object in the field of view .
the results of our analysis are listed in table  [ tabros ] .
lllrl@r name & d & jd & & ctrate & + & ( pc ) & & & ( s@xmath5 ) & +   +   +   + & 107@xmath6212 & 8420.575 & 16156 & 0.322(5 ) & 36 + & & 8775.123 & 8345 & 0.126(5 ) & 36 + & 21.7@xmath620.5 & 8864.531 & 4267 & 6.072(15 ) & 0 + & & 9003.872 & 5776 & 6.228(11 ) & 0 + & & 8864.595 & 3651 & 6.84(2 ) & 24 + & & 9004.871 & 3151 & 5.41(2 ) & 24 + & 37.7@xmath621.9 & 9049.117 & 14531 & 1.206(4 ) & 0 + & 16.4@xmath620.2 & 9247.809 & 9895 & 1.009(7 ) & 0 + @xmath63 & 21.4@xmath620.5 & 9085.714 & 2140 & 1.10(2 ) & 0 + & 49.9@xmath622.1 & 9148.099 & 1209 & 1.23(3 ) & 0 + & 25.3@xmath624.9 & 8604.712 & 5249 & 0.888(16 ) & 36 +   +   + & 29.0@xmath627 & 8648.446 & 3098 & 6.29(2 ) & 0 + & 168@xmath6251 & 8426.486 & 21428 & 0.146(3 ) & 33 + & 108@xmath6212 & 8810.279 & 2526 & 0.336(17 ) & 44 + & & 8991.040 & 6214 & 0.372(11 ) & 44 + & & 8796.874 & 5076 & 0.305(12 ) & 51 + & & 8966.506 & 2904 & 0.366(18 ) & 51 + & 44.5@xmath621.2 & 8431.422 & 71803 & 2.3152(14 ) & 44 + & & 9158.435 & 37658 & 2.403(3 ) & 44 + @xmath64 & 202@xmath6257 & 9030.359 & 10327 & 0.059(2 ) & 0 + @xmath64 & 193@xmath6239 & 8789.392 & 8668 & 0.198(5 ) & 0 + @xmath65 & 42.0@xmath621.0 & 8620.767 & 13460 & 1.717(4 ) & 0 + & & 9136.873 & 4892 & 2.902(13 ) & 0 +   +   + & 192@xmath6234 & 8322.477 & 3516 & 0.574(13 ) & 3 + @xmath66 & 37.5@xmath621.1 & 8346.927 & 4745 & 5.269(14 ) & 0 + & & 8904.569 & 1438 & 3.43(4 ) & 0 + & & 8735.159 & 7940 & 4.753(8 ) & 0 + & 138@xmath6210 & 9272.484 & 4156 & 1.32(2 ) & 41 + & 106@xmath627.6 & 8328.919 & 14690 & 0.624(6 ) & 40 + @xmath67 & 103.@xmath628.5 & 8863.981 & 10576 & 1.554(8 ) & 31 + @xmath67 & 323@xmath6253 & 9088.631 & 24889 & 0.112(2 ) & 0 + & 96.8@xmath627.1 & 8589.200 & 6335 & 1.692(15 ) & 44 + & & 8769.630 & 8150 & 1.422(12 ) & 44 + & & 8973.263 & 22143 & 1.909(4 ) & 44 + @xmath68 & 25.8@xmath620.5 & 8448.155 & 31165 & 4.077(2 ) & 0 +    to compare systems at different distances , we multiply the rosat countrate and the flux of the emission cores for each system with the square of the distance listed in table  [ tabros ] ; for we adopt a distance of 850pc ( twarog & anthony - twarog 1989 ) .
no corrections are made for interstellar absorption .
the choice of the 0.42.4kev bandpass minimizes the effects of interstellar absorption , which are severe at energies @xmath69kev .
as it is unknown which component of the binary emits the x - rays , we plot the total x - ray and fluxes , adding the contributions of both components where these are given separately by fernndez - figueroa et al .
( 1994 ) .
the resulting absolute countrates and fluxes are shown in fig .
the systems with h&k emission clearly visible in fig .
[ cahk ] , viz .
s1063 , s1113 , s1040 and s1242 lie on the relation between x - ray and h&k emission defined by the rs cvn systems , in agreement with the hypothesis that the x - ray flux of these objects is related to the magnetic activity .
it is also seen that the upper limits or marginally detected emission cores in s1082 , s1072 and s1237 are high enough that we can not exclude the hypothesis that the x - ray emission in these systems is related to magnetic activity .
the rotational velocity is another indicator of magnetic activity .
we investigate the relation between rotational velocity and ca emission by selecting those stars from the sample of fernndez - figueroa for which a value of @xmath48 is given in the catalogue of chromospherically active binary stars ( strassmeier et al . 1993 ) . in fig .
[ vsinica ] the h&k emission of these stars is compared with their @xmath48 . in this figure we do discriminate between the separate contributions of both stars to @xmath31 , with the exception of s1113 for which we combine the total flux @xmath31 with the @xmath48 of the primary .
the stars are found within the range occupied by chromospherically active stars .
we note that the correlation between the observed h&k flux and @xmath48 is not tight .
in particular , high and low h&k emission flux is found at low values of @xmath48 .
some of the scatter may be due to the use of @xmath48 instead of the stellar rotation period .
parameters depending on the spectral type ( e.g. properties of the convective region ) have been used to reduce the scatter in the activity - rotation relation ; whereas this is successful for main - sequence stars with @xmath70 , it fails for other main - sequence stars and for giants ( see discussion in stpie 1994 ) .
for example , the three giants ( k0 iii ) , ( k0 iii ) and ( k0 - 2 iii ) have @xmath48 values of 10 , 10 and 8 km s@xmath5 , respectively , but differ in @xmath71 by three orders of magnitude ( see fig .
[ vsinica ] ) .
the h@xmath0 lines ( @xmath72 6562.76  ) of only two stars , s1063 and s1113 , show clear evidence of emission , as shown in fig .
[ subha ] .
this is described in more detail in their individual subsection in sect .
4.3 . for the other stars ,
we have used the few ( sub)giants in the library of ues - spectra of montes &  martn ( 1998 ) to investigate the behaviour of h@xmath0 .
we have chosen library spectra of stars that match the spectra of the stars as closely as possible ( see fig .  [ halpha ] ) .
s1072 and s1237 show no evidence for filling in of the h@xmath0 profile compared to a g0iv  v and a k0iii star , respectively . in s1040 h@xmath0
seems slightly filled in compared to a g8iv and a k0iii star .
this is also the case for s1242 compared to a g0iv - v star , but we note that no classification for this star is found in literature . for s1082 ,
no matching spectrum is available in this wavelength region .    filling in of the lines in the b triplet
( @xmath15 5167.33 , 5172.70 and 5183.62  ) and in the d doublet ( @xmath15 5889.95 and 5895.92
 ) is visible in some active stars .
the presence of a d@xmath73 ( @xmath72 5876.56 ) absorption or emission feature can also indicate activity ( see discussion in montes &  martn ( 1998 ) and references therein ) .
however , in none of the stars we see filled in b and d lines .
neither do we see a clear d@xmath73 feature . for s1082 ( b and d ) and s1113 ( b , d and d@xmath73 )
we find no suitable library stars for these features .
the two stars below the subgiant branch , s1063 and s1113 , both show relatively strong h&k emission , and are the only two stars in our sample showing h@xmath0 in emission , shown in fig .  [ subha ] .
we use the orbital solutions for both objects to try and identify the star responsible for these emission lines .
the velocities of both components in s1113 , and of one component in s1063 are indicated in fig .
[ subha ] .
the h@xmath0 line profile of s1063 is asymmetric , showing emission which is blue - shifted with respect to the absorption .
the location of the absorption line is compatible with the velocity of the primary , which dominates the flux ; the emission is probably due to the secondary . remarkably , the h&k emission peak is at the velocity of the primary .
this suggests that the h@xmath0 emission is not chromospheric in nature .
the h@xmath0 emission of s1063 does not show the double peak that is known to indicate accretion disk emission ( horne &  marsh , 1986 ) .    in s1113
the h@xmath0 emission profile is symmetric and broad , with full width at continuum level of 15  .
the emission peak is centered on the more massive star , which contributes 82%  of the total flux ( table  [ fluxes ] ) .
this suggests that the h@xmath0 emission is due to the primary .
the h&k emission shows marginal evidence for a double peak , suggesting that both stars contribute to the chromospheric emission . in fig .
[ subca ] we indicate the expected position of the h&k lines for both stars . for the phase observed , their peaks overlap in the crosscorrelation function . in the figure we use the velocities resulting from fitting the order that gives the cleanest crosscorrelation .      the h&k emission in the wide binaries s1072 and s1237 is only marginally significant .
the level of their x - ray and emission is more appropriate for active main - sequence stars ( fig.[xca ] ) .
one might speculate that it is due to the invisible companion of the giant detected in the crosscorrelation ; even if this were the case , we would not understand why this companion would be chromospherically active .
we conclude that we do not understand why these two stars are x - ray sources .
we find no indication for a faint secondary in the crosscorrelation profile of s1237 ; at the time of observation , the spectra of two equally massive stars as suggested by janes &  smith ( 1984 ) would be separated by 4 to 7 km s@xmath5 ( derived from the ephemeris in mathieu et al .
since the secondary is 1.6 magnitude fainter in v than the primary , we think that this small separation is compatible with finding a single peak in the crosscorrelation .
s1242 is chromospherically active , as shown by its h&k emission .
we suggest that this activity , which also explains the x - rays , is due to rapid rotation induced by tidal interaction at periastron , which tries to bring the subgiant into corotation with the orbit at periastron . if we assume that the observed period of photometric variability is the rotation period we derive an inclination of @xmath7 9@xmath74 using our maximum value of @xmath48 and the estimated radius
. this would be in agreement with a companion at the high end of the range 0.140.94 @xmath75 allowed by the mass function ( mathieu et al .
1990 ) .
our detection of clear chromospheric emission indicates that the x - ray emission of s1040 is due to the giant .
the white dwarf has a low temperature , and is unlikely to contribute to the x - ray flux .
we find a rather slow rotational velocity for the giant , 3 km s@xmath5 .
gilliland et al .
( 1991 ) detected a periodicity of 7.97 days in the visual flux ( b and v bandpasses ) of s1040 , with an amplitude of 0.012mag . if this is the rotation period of the giant , the radius of 5.1@xmath76 ( landsman et al .
1997 ) implies an equatorial rotation velocity of @xmath77 km s@xmath5 .
this is compatible with the velocity measured with our crosscorrelation , @xmath48 , for an inclination @xmath78 .
this inclination has an a priori probability less than 0.5% ; and it implies an unacceptably high mass for the white dwarf , from the measured mass function @xmath79 .
we conclude that the 8days period can not be the rotation period of the giant .
it is doubtful that the white dwarf can be responsible , as its contribution to the @xmath2 and @xmath1 flux is small .
the h@xmath0 absorption profile of the blue straggler s1082 is variable . if we consider the most symmetric spectrum profile , that of 00:01 ut , as the unperturbed profile of the primary , we find that the changes are due to extra emission .
this is illustrated in fig .
[ s1082 ha ] .
we suggest that this variation is due to the subluminous companion , possibly to a wind of that star .
we have also investigated the presence of a broad shallow depression underlying the d lines ( near @xmath72 5895  ) and the triplet ( near @xmath72 7775  ) as found by mathys ( 1991 ) .
we find that this broad component is variable , as illustrated in fig .
[ s1082na ] .
mathys ( 1991 ) suggests that the broad component originates in the subluminous companion .
this companion outshines the primary by a factor six at @xmath72 1500  and thus is presumably hot ( landsman et al . 1998 ) .
we note that the star can not be too hot or it would not show neutral lines .
in this paper we have tried to find an explanation for the x - ray emission of seven sources in .    for s1242 and s1040
we have concluded from the h&k emission cores that magnetic activity is responsible for the x - rays .
this is supported by filling in of h@xmath0 ( see e.g. montes et al .
1997 ; eker et al . 1995 ) . in s1242
, activity is likely to be triggered by interaction at periastron in the eccentric orbit .
this is also reflected in the period of photometric variability . for s1040
, the reason for activity is less clear .
the explanation could involve mass transfer from the precursor of the white dwarf to the giant and the latter s subsequent expansion during the giant phase . as was already noted by landsman et al .
( 1997 ) , a similar system is , a binary of a white dwarf of @xmath80 k ( simon et al .
1985 ) and a g5iii giant in an 56.8 days circular orbit .
the @xmath48 of that giant is also low , 4 km s@xmath5 , and the long photometric period of 77.2 days implies asynchronous rotation ( strassmeier et al .
the x - ray luminosity for is @xmath81 erg s@xmath5 in the 0.24 kev band as measured with einstein by walter &  bowyer ( 1981 ) , somewhat higher than the luminosity of s1040 .
( with the coronal model discussed by belloni et al .
( 1998 ) , the countrate for s1040 corresponds to @xmath82 erg s@xmath5 in the 0.24 kev band . )
walter &  bowyer attribute the x - rays to coronal activity of the giant .
the h&k emission cores in s1063 and s1113 are very strong . in s1113
we might even see emission of both stars .
due to the shape of the h@xmath0 emission we can not conclude with certainty that the x - rays arise in an active corona and not in a disk or stream .
the wings in the emission peak of s1113 are very broad .
however , montes et al .
( 1997 ) have demonstrated that the excess emission in the h@xmath0 lines of the more active binaries is sometimes a composite of a narrow and broad component , the latter having a full width at half maximum of up to 470 km s@xmath5 .
they ascribe this broad component to microflaring accompanied by large scale motions .
we note the similarity between s1113 and , a well known extremely active binary of a g5iv ( @xmath48 = 13 km s@xmath5 ) and k1iv ( @xmath83 km s@xmath5 ) star in a 2.84 days circular orbit and a mass ratio 0.79 ( strassmeier et al . 1993 ) ; the mass ratio of s1113 is 0.70 ( mathieu et al .
1998 , in preparation ) . from the countrate of in table  [ tabros ] we find @xmath84 erg s@xmath5 in the 0.12.4 kev band using the same model as in belloni et al .
( 1998 ) with @xmath85 , which is comparable to the luminosity of s1113 in the same band @xmath86 erg s@xmath5 ( belloni et al .
the h@xmath0 emission of s1063 is more difficult to explain .
as this system is not double - lined , h@xmath0 emisson by the ( invisible ) secondary star would have to be strong to rise above the continuum of the primary .    in the binaries
s1072 and s1237 we see no h@xmath0 emission while the level of h&k emission is low in comparison with active stars of the same luminosity class .
we have no explanation for this . for s1072 ,
an option is a wrong identification of the x - ray source with an optical counterpart .
belloni et al .
( 1998 ) give a probability of 43% that one or two of their twelve identifications of an x - ray source with a binary in is due to chance .
no h&k emission is seen in the spectrum of the blue straggler s1082 .
possibly , the x - ray emission has to do with the hot , subluminous secondary that could also cause the photometric variability and whose signature we might have seen in the h@xmath0 line .
the authors wish to thank g. geertsema for her help during our observations , p. groot for providing his program to compute projected rotational velocities with the fourier - bessel transformation method and m. van kerkwijk for comments on the manuscript .
mvdb is supported by the netherlands organization for scientific research ( nwo ) .
carter , d. , benn , c.  r. , rutten , r. g.  m. , breare , j.  m. , rudd , p.  j. , king , d.  l. , clegg , r. e.  s. , dhillon , v.  s. , arribas , s. , rasilla , j .- l . , garcia , a. , jenkins , c.  r. , charles , p.  a. 1993 , isis users manual , http://www.ing.iac.es                                                                                zimmermann , h. , becker , w. , belloni , t. , dbereiner , s. , izzo , c. , kahabka , p. , schwentker , o. 1994 , exsas user s guide : extended scientific analysis system to evaluate data from the astronomical x - ray satellite rosat , technical report 257 , mpe