high precision doppler spectroscopy is an important tool to investigate many astrophysical processes in the universe .
measurement of tiny doppler shifts from radial velocity ( rv ) can provide insight about the stellar interior ; reveal the existence of unseen planets around other stars ; help us probe the variability of fundamental constants of nature and provide useful information of binary and pulsating stars @xcite .
in particular , the radial velocity technique has been immensely successful in detecting hundreds of planets with different mass range and orbital periods @xcite . the main challenge and focus of contemporary research
is to find earth - analogs within the habitable zone of star where life sustaining conditions could prevail .
recently , the doppler data from eso s two precision rv instruments , namely harps and uves , was decisive in unveiling the presence of 1.3 earth - mass rocky planet orbiting proxima centauri - the nearest known star to the sun @xcite .
typically , doppler precision of @xmath3 cms@xmath2 or better is required to address many compelling science cases including the detection of terrestrial - size planets around sun - like stars @xcite .
leading rv observing campaigns are carried out with state - of - the - art echelle spectrograph , carefully designed and built to operate under tightly regulated thermovac conditions @xcite .
a high resolution spectrum of star in visible / nir band is recorded on a 2d array of ccd pixels . to extract the rv signal from time series data , spectrograph is wavelength calibrated with a laboratory source of well known emission or absorption lines . in calibration process several lines of known wavelengths are identified and a quadratic or high order polynomial fit is used to invert the pixel coordinates into accurate wavelength units .
two commonly used wavelength reference in astronomy are thorium - argon ( th - ar ) emission lamp and iodine absorption cell @xcite . despite averaging the shift over thousands of spectral lines , the best rv precision obtained from these sources is limited to @xmath4  ms@xmath2 .
the rv uncertainty of traditional sources arises from wavelength inaccuracies of individual lines , spectral contamination , aging of the lamp , limited spectral coverage , line blending effects , uneven distribution and varying intensity of lines across spectrograph band @xcite .
an ideal calibrator choice for enabling sub - ms@xmath2 rv precision in next generation spectrographs is a broadband laser frequency comb ( lfc ) .
the lfc consists of thousands of equispaced lines produced from a femtosecond mode - locked laser @xcite .
the entire comb spectrum is stabilized by anchoring to atomic clock with well established radio - frequency locking technique . since past decade different versions of lfcs
have been tested and validated @xcite .
new developments are still awaited to overcome a couple of technical and operational difficulties .
it might take a while for lfc to become a simple and affordable turn - key instrument for routine wavelength calibration at most observatories .
a simple and cost effective approach to produce reference spectra is to use cavity resonance lines of a fp etalon @xcite .
a plane parallel fp cavity when illuminated with a broadband white light source acts as a periodic wavelength filter .
the transmitted light consists of evenly distributed lines that result from multiple beam interference inside the cavity @xcite .
the cavity lines are separated by free spectra range @xmath5 , where @xmath6 is the light speed and @xmath7 is the cavity length .
many variants of fps ( e.g. solid , air spaced and fiber based ) have been successfully built and tested for astronomy calibration @xcite .
the positional accuracy of fp lines critically depends on stability of the cavity itself .
any noise affecting the cavity length would randomly displace the resonant lines thus directly impairing the achievable rv precision .
the impact of temperature and pressure fluctuations on the cavity length is generally eliminated by mounting the fp inside a thermally controlled evacuated chamber which is capable of achieving @xmath8  1 ms  @xmath2 short - term stability @xcite .
a passively stabilized fp apparatus producing reference spectra has no built  in mechanism to track systematic offsets or drifts in cavity length occurring over time scales of weeks , months or even years .
regardless of good passive isolation , cavity drift could still occur from temperature residuals , material aging or changes in optical properties of the coatings @xcite .
the radial velocity @xmath9 and the fractional length and frequency stability of the cavity can be expressed as @xmath10 .
for instance , a 1  nm change in a 1  cm long cavity corresponds to @xmath855  mhz frequency shift of cavity lines in visible ( @xmath11 nm ) region , which amounts to 30  ms@xmath2 rv uncertainty . to reach a rv precision @xmath12 cms@xmath2 , the cavity length should be stable to @xmath13 level or better
.    for long - term usability ( several months to years ) of astro - fp , research efforts are now being directed at controlling and tracking the time dependent cavity drift with advance techniques @xcite .
measurement of tiny dimensional changes in cavity length calls for advanced metrology based on laser frequency standards . a drift tracking or active length - control
requires cavity - locking to a frequency stabilized diode laser or vice - versa @xcite . in cavity - lock setup
, the frequency stability of the laser is expected to be transferred to the stability of the fp cavity lines .
the main task of frequency or phase locking mechanism is to keep cavity transmission point in resonance with frequency of the laser @xcite .    a short - length , low - finesse(@xmath14 ) cavity , while useful for generating broadband reference lines for astronomy spectrographs , is not ideal for laser metrology .
normally , a robust laser - lock requires a sharp discriminant signal which depends on the laser power and the cavity linewidth @xcite .
the linewidth @xmath15 of the cavity should be sufficiently narrow to produce a steeper error signal .
the inverse dependence of @xmath7 and @xmath16 on the linewidth @xmath17 implies a clear advantage of longer and high finesse cavities .
another practical difficultly of utilizing astro - fp in laser applications is the cavity geometry .
a confocal cavity design is more appropriate for locking purpose because it greatly simplifies beam alignment and mode - matching conditions for single frequency laser beam . on the contrary , a collimated white light beam
is easily coupled into a short - length planar cavity producing reference spectra @xcite .
these considerations show incompatible requirements of a cavity for laser metrology and astronomy calibration .
clearly , a single fp cavity is inadequate for serving dual purpose of ` _ _ wavelength calibration _ _ ' on one hand and ` _ _ drift tracking _ _ ' on other .
we propose to circumvent this problem by constructing two adjoining cavities of equal length but of different finesse on a single block of ultra low expansion ( ule ) glass . the low finesse _
astro - cavity _ would produce the calibration spectra while the high finesse _ lock - cavity _ would facilitate frequency locking to implement stability tracking mechanism .
the length of two cavities is defined by the length of monolithic spacer block . here
, the implicit assumption is that external perturbations , if any , act uniformly on the entire fp block , thus causing equal and identical length changes in both cavities .
in fact , drift measurements would only be accurate so long the relative length variation between two cavities is well below @xmath18 . under this premise alone the _ lock - cavity _ would function as a reliable proxy for measuring the stability of the _ astro - cavity_. departure from this scenario would result in intractable measurement errors .
mm on top are provided for the venting purpose .
( b ) the front view showing the spacer dimensions along with two cavities of diameter @xmath19  mm each and standard hr mirrors with @xmath20  mm each .
( c ) the isometric view .
the cavity axis is orientated along the @xmath21-direction and the lab vertical is along @xmath22-direction . ]    the main objective of this paper is to numerically examine the validity of our approach .
we identify the temperature and vibration induced noise conditions that impart unequal length changes in two adjoining cavities constructed on a monoblock spacer .
these studies are important to understand the level of performance one can realistically expect to achieve in the proposed device .
we use finite element method to numerically simulate the impact of small thermal gradients and orientation dependent vibration noise that could deforms two cavities in dissimilar manner .    from here on , the paper is organized as follows . the dual fp design and laser - lock setup for cavity drift tracking is discussed in section 2 . in section 3 , the mathematical formulism and fea model for cavity deformation is described .
the vibration analysis of horizontally mounted cavity is presented in section 4 .
a proposed vertical mounting configuration and numerical optimization procedure to reduce the vibration sensitivity of the fp device is given in section 5 .
thermo - elastic model is used to calculate the temperature sensitivity of the fp device .
effect of temperature gradient and cavity length offset is presented in section 6 .
numerical accuracy of fea results is discussed in section 7 and finally , the summary of the work is given in section 8 .
fabry - perot cavities are used in diverse range of applications in many areas of fundamental and applied research @xcite .
ultrastable cavities are also key to gravitational wave detection , testing of lorentz invariance and experiments in future space missions @xcite . for wavelength calibration , the fp cavity spacing and mirror reflectivity @xmath23
are selected based on spectrograph requirements .
the cavity spacing @xmath7 is chosen so that reference lines are easily resolved by spectrograph in the detector plane .
generally , a short cavity ( e.g. , @xmath24 cm , @xmath25 ghz ) is required to generate optimally spaced comb lines to calibrate a high resolution ( @xmath26 ) astronomy spectrograph .
the reflectivity ( @xmath27 ) of broadband coated mirrors yields a low finesse @xmath28 .
the resulting cavity linewidth @xmath29 is such that it is optimally sampled by the detector pixels . to produce reference lines that are 15 ghz apart
would require a planar cavity of length 1  cm .
the main component of the proposed device is a two - channel fp etalon .
it consists of two cylindrical cavities of diameter 12  mm each , bored adjacent to each other on a rectangular block ( @xmath30 mm@xmath31 ) of ultra low expansion ( ule ) glass .
the cavity mirrors are made of fused silica ( fs ) and optically contacted on either side of the spacer block .
we assume the optical contacts to be sufficiently strong and free from any imperfections .
this is a reasonable assumption as the modern techniques of joining optical elements have become more robust @xcite . a cad geometry of the device is shown in fig .  1
. the low - finesse _ astro - cavity _ ( fp1 ) is to provide astronomy calibration lines and high - finesse _ lock - cavity _ ( fp2 ) is to provide a light channel for robust laser - lock .
the cavity length @xmath32=1  cm is chosen to produce 15  ghz apart reference lines to calibrate a typical high resolution ( 100,000 ) echelle spectrographs . to maintain sufficient rigidity of the entire structure ,
a small 5  mm gap is provided on all sides between mirrors and spacer edges .
cavity mirrors are assumed to be flat .
replacing flat mirrors with curved mirrors will not alter our analysis in any major way , except a very small change in curvature due to differential thermal expansion .     generated by electro - optic modulator ( eom ) .
the rf electronics following the photodiode ( pd ) provides necessary feedback to keep one of the sideband frequencies , e.g. @xmath33 ( shown in red ) , locked to the peak of cavity transmission line ( shown in blue ) . for clarity ,
other ancillary components are not shown . ]
the conceptual principle of dual fp cavity - lock is illustrated in fig .
2 . the cavity drift can be tracked in several ways , but often requires two lasers @xcite .
however , the same goal can be accomplished using a simpler approach of _ sideband locking _ which needs only one laser @xcite . in this method
the frequency @xmath34 of the laser is pre - stabilized to rb - cell using standard saturation absorption @xcite . a tunable electro - optic modulator ( eom )
creates radio frequency ( rf ) sidebands @xmath35 around the central frequency @xmath34 of the laser .
one of the sidebands is then locked to the nearest resonance peak of fp2 cavity using modified pound - drever - hall ( pdh ) method .
the rf frequency @xmath36 of the oscillator driving eom is servo - locked to cavity .
this technique preserves the intrinsic stability of main laser beam @xmath34 that is locked to rubidium while exploiting the sideband tunability to track the cavity drift .    to estimate the drift , the frequency of the radio oscillator has to be readout independently by a frequency counter or rf spectrum analyzer .
a fractional length change ( @xmath37 ) about @xmath18 of the cavity block corresponds to rv error of 10  cms@xmath2 or a rf frequency shift of few hz in an eom operating at 1  ghz . measuring a ghz rf signal with an accuracy of few hz ( amounting to few ppb frequency resolution )
is a challenging task .
however , the desired precision can be reached with many cycles of repeated measurements .
alternatively , optical - beat measurement with two lasers yields much better sensitivity when the cavity is used in reflection mode , i.e. pdh locking @xcite .
for example , a fractional change @xmath18 of the cavity length would induce a readily measurable frequency excursion of @xmath38  khz in the optical - beat signal obtained with pdh setup using 780 nm laser .
the exact implementation details will be worked out in a separate study .
the resonance stability of fp cavity is determined by the stability of optical length between each mirror pairs .
different noise sources e.g. , seismic and acoustic vibrations , material inhomogeneities , fluctuations in temperature , pressure and refractive index of the air etc influence the cavity stability .
the impact of pressure on refractive index fluctuations is largely eliminated by placing the cavity inside a vacuum tank .
a properly designed vacuum system also attenuates acoustic noise .
our focus in this study is on temperature and vibration perturbations which affect the long- and short - term stability of the fabry - perot cavity . to achieve a desired tracking accuracy the length changes @xmath39 in astro - cavity fp1
must strongly correlate with the length changes @xmath40 in lock - cavity fp2 . how well @xmath39 and @xmath40 correlate would largely depends on cavity geometry and the way noise is transmitted into the experiment . normally , the thermal and acoustic noise couples into cavity through mounting and support system .
the main purpose of this paper is to estimate noise levels for which the differential length change in two cavities , @xmath41 becomes significant . in other words ,
when @xmath42 , the change in rf lock - point of fp2 would accurately reflect the length changes in fp1 .
uncorrelated cavity deformations would render the laser - lock technique ineffective .
the shape of an elastic body can deform under the application of a force that may arise from constant acceleration ( e.g. gravity loading ) , floor - vibrations , shocks or acoustic pickup . in most cases ,
acceleration experienced by a cavity is in the form of vibration perturbations .
these perturbations are time - dependent and contain a large range of frequencies that excite various mechanical eigenmodes of the cavity .
the exact nature of interaction depends on wavelength of excitation relative to the physical size of the cavity @xcite .
in high - frequency limit ( small wavelength ) , a full dynamic analysis is needed to obtain a complete information of local strain and displacement of the cavity points .
since the excitation wavelengths at higher frequencies are significantly smaller than the cavity size , the conventional numerical discretization methods for dynamic analysis are inefficient and require significant computer resources .
though from practical viewpoint , the high frequency vibrations are less of a problem . they can be effectively attenuated in lab using vibration isolation tables , suitable damping materials , soft cushioning and multi - stage suspension .    in low - frequency limit ( @xmath43 hz )
the dynamic problem of cavity deformation can be reduced to static analysis @xcite . in this regime
the excitation wavelength become longer than the cavity size and all points in the cavity respond to the applied force in unison manner .
the low - frequency collective motions ( compression or stretching ) of the particles can be approximated by quasi - static equilibrium resulting from application of a static force or a constant acceleration .
traditionally , the vibration or acceleration sensitivity of a cavity is expressed in normalized units of acceleration due to gravity @xmath44  [ ms@xmath45 . in this work
we also follow the same convention .
elastic solids change shape under the influence of external force .
the elastic deformations can be quantified by knowing the displacement of material points inside the body .
mathematically , elasticity problem can be formulated in terms of field equations which govern the relations among displacements , strains and stresses inside a loaded object .
the basic field equations for linear isotropic solids comprise @xcite :    * @xmath46__strain - displacement relation ( compatibility equations ) _ _ @xmath47 which completely specifies the strain @xmath48 in terms of displacement @xmath49 and its derivatives . + * @xmath46__equations of equilibrium _
_ @xmath50 which is a statement of newton s 2nd law for static cases where displacement has no time dependence and derivatives of stress tensor @xmath51 and body forces @xmath52 add to zero .
+ * @xmath46__constitutive equations _ _ ( hooke s law ) @xmath53 which establishes general stress - strain behaviour via material parameters , i.e. , young s modulus @xmath54 and poisson s ratio @xmath55 .    .
the entire structure is vertically compressed and horizontally expanded from original configuration .
( c ) the z - displacement and tilt of the cavity mirrors ( magnified by @xmath56 ) from their nominal positions .
the black outline indicates undeformed shape of the cavity .
the color scale represents the cavity displacement in @xmath21-direction . ]    but for simple geometries , the solution of general system of eqs.(1)-(3 ) by analytical methods is almost impossible .
sometimes further simplifications are made by developing a reduced set of field equations solely in terms of stresses or displacements .
one such formulism is to eliminate stresses and strains from eqs.(1)-(3 ) and express the displacement in a simplified form , known as navier s equation .
the vector form of navier s equation is @xcite : @xmath57 where @xmath58 is shear modulus and @xmath59 is _ _ lam__s constant
. equation ( 4 ) is usually solved with displacement boundary conditions of the form : @xmath60 , where @xmath61 denotes the boundary points on the surface of the body and @xmath62 are prescribed displacement values .    as a rough approximation ,
the cavity deformation problem is similar to stretching or compression of a prismatic bar when acted upon by an external force . for illustration
, we consider a simple rectangular block of length @xmath63 resting horizontally as shown in fig .
such a block will be compressed by its own weight with body forces @xmath64 and @xmath65 , where @xmath66 is the mass density and @xmath44 is acceleration due to gravity .
this problem can be solved by direct integration of field equations . to prevent the translation of the bar we fix the central point @xmath67 at the lower end of the bar , that is , @xmath68 and @xmath69 .
further , the rotation of the bar about @xmath67 and about an axis parallel to @xmath67 is prevented by setting @xmath70 at point @xmath67 .
final expressions for three displacement components can be written as : @xmath71\\    w & = &   \nu \frac{\rho g\,z}{e}(l - y)\end{aligned}\ ] ] where @xmath66 is the material density , @xmath44 is the acceleration , @xmath55 is poisson s ratio , @xmath63 is dimension of the bar along the direction of the force . from eqs.(5)-(7 ) and fig .
3(a ) we note that points on the y - axis are only vertically compressed , i.e. , @xmath72/2e$ ] .
for other points on the bar , the vertical contraction @xmath9 is accompanied by a finite stretching @xmath73 and @xmath74 in other two orthogonal directions .
the coupling of the displacement field is mediated by non - zero value of poisson s ratio .
the example of a prismatic bar illustrated in previous section has an analytical solution due to simple geometry .
the actual structure of the fp is quite complex .
it has two large cavities and four mirrors attached to the spacer .
the dimensions as well as mass of the mirrors are also comparable to the spacer .
using analytical means to find a closed - form solution for complicated geometry is extremely difficult task .
we , therefore , used finite element analysis ( fea ) to solve our 3d problem numerically .
the fea is by far the most widely used and versatile numerical technique for simulating deformable solids .     of horizontal ( x - cut ) and vertical ( y - cut ) probe points .
the displacement of any other point on the cavity mirrors is confined to the gray area between the curves .
@xmath75 denotes the displacement of the mirror centers . ]
we chose comsol multiphysics as our fea simulation software .
it provides a complete and integrated environment for creating , analyzing and visualizing multiphysics models . built - in cad tools and the _ structural mechanics _ module was used for constructing a 3d model of the fp cavity .
additional quantities and model inputs ( e.g. boundary constraints , body load , material type and properties etc . ) for stress and displacement analysis were defined in _
solid mechanics _ interfaces of the _ structural _ module .
the relevant material properties used in this study are listed in table 1 .
the mesh settings determine the resolution of the mesh element used for discretizing the model . in all simulations we used default tetrahedral mesh shape .
extra - fine mesh size was chosen to minimize the discretization error . for our study
, we selected the stationary solver since our load , stress and deformation do not vary with time .
the _ solid mechanics _
kernel numerically solves the navier s equations ( eqs.(4 ) ) and computes displacements , stresses , and strains over the entire fp domain .
the deformation of a material point ( @xmath76 ) in fp geometry is quantified by calculating the displacement ( @xmath77 ) from the undistorted shape . in subsequence analysis , our prime interest will be in displacement component @xmath78 which is along the axis of the cavity , i.e. @xmath21-direction .
.mechanical and thermal properties of the materials used in simulations [ cols= " < , < , < , < , < " , ]
the fp is mounted horizontally in @xmath79-plane as shown in fig .
1 . the cavity axis is along @xmath21-direction .
since the cavity length is short we have not considered any discrete support for horizontal mounting .
instead we let the spacer block rest on the bottom plane .
first we investigate a simple case where the fp unit is deformed under its own weight , i.e. , force of gravity @xmath44 acting uniformly on entire geometry in @xmath80 direction .
fixed boundary constrain @xmath81 is imposed on the bottom @xmath79-plane of the spacer to prevent rigid body motion in fea simulations .
cavity deformation in @xmath82-plane due to gravity load is illustrated in fig .
in accordance with eqs .
( 3)-(5 ) , the gravity - induced compression of the cavity in vertical direction also causes it to expand in @xmath21-direction , thus increasing the mirror separation . in order to numerically compute the cavity length changes
, we extract the displacement component @xmath74 for each mirror surface that is bonded to the spacer block .
the numerical probe points are defined along @xmath83 and @xmath22-direction ( see the fig .
4 inset ) on the interior surface of the mirrors forming cavities .
the numerical displacement @xmath74 of the probe points , namely the @xmath83-cut and @xmath22-cut , is shown in fig .
4 . in simulations
the origin of the coordinate system was chosen to coincides with the geometrical center of the fp spacer block . a near overlap between the fp1 and fp2 probe points
indicate that both cavities are deformed identically under the force of gravity .
it is also evident from fig .
4 that displacement @xmath74 is not uniform across the cavity cross - section .
the horizonal bulge shown by concave shaped red curves , is symmetric about the mirror centers where the displacement is maximum ( @xmath84  m ) and reduces to about @xmath85  m towards to cavity edge .
the displacement @xmath74 of the vertical probes , i.e. , the @xmath22-cut shown in blue , is asymmetric .
it is minimum ( @xmath86  m ) at the cavity top edge , attains maxima ( @xmath87  m ) at about 3  mm below the mirror centers and reduces slightly towards the edges .
this asymmetry in displacement is because of the unequal weight along vertical direction and a small @xmath22-tilt caused by outward buckling of the spacer and the mirrors .
the @xmath21-displacement of any other point on the mirror surface is bound by gray area between the blue and red curve .
since the displacement @xmath74 is not uniform across the reflecting face of the mirror we express the effective cavity length change @xmath88 for each cavity as the difference in displacement of the mirror centers .
finally , the quantity of interest for us is differential length change between two fp cavities , i.e. @xmath89 . for the present case , we find @xmath90pm and @xmath42 .
this is expected as the applied gravity load is symmetric and both cavities are deformed equally .     with gravity - like acceleration vector @xmath91 oriented in ( a ) @xmath92-plane , ( b ) @xmath82-plane and ( c ) @xmath93-plane of the dual cavity .
for each plane the angle @xmath94 and corresponding 2d projection of the fp is indicated in the inset figures .
the orientation angle @xmath94 was swept from @xmath95 in step of @xmath96 .
( d ) the differential length change @xmath97 between the two fp cavities as a function of sweeping angle . the y - scale to the right gives corresponding rv tracking error ( cms@xmath2 ) in each case . ]      now we examine cavity s response to a gravity - like static force that is orientation dependent .
this exercise is aimed at understanding the relative sensitivity of the fp to direction dependent noise perturbations . in simulations
, we swept the acceleration vector @xmath91 in different planes and computed differential length change @xmath98 for each orientation .
the acceleration @xmath91 oriented in a given @xmath99-plane , can be expressed in its components along @xmath100 and @xmath101 directions , i.e. , @xmath102 , where for consistency we have chosen @xmath103 . using parameter sweep feature of the comsol we simulated the impact of different orientation of @xmath91 in the @xmath99-plane .
cavity length change @xmath88 and @xmath98 as a function of orientation in different planes is shown in fig .
the impact of orientation dependent force can be understood by analyzing its response on cavity length separately along each direction .
a few noteworthy observations are :    * when the acceleration is along the @xmath83-direction ( @xmath104[@xmath105 in fig .
5(a ) and @xmath106[@xmath107 in fig .
5(c ) ) the sheared strain in fp is produced in such a way that fp1 cavity mirrors are deflected inward [ outward ] and fp2 cavity mirrors are deflected outward [ inward ] .
such response can be attributed to unequal shear strain produced because of the constrained boundary and geometrical asymmetry of the cavity in @xmath83-direction . a significant differential length change ( @xmath108pm )
occurs due to shortening of one cavity and elongation of the other .
this amounts to spurious rv tracking errors about 140 cms@xmath2@xmath109 or equivalently a frequency shift of 1.8 mhz@xmath109 at 782 nm laser . *
horizontally mounted dual fp is largely immune to relative length change in vertical vibrations .
when the acceleration is oriented in @xmath22-direction ( @xmath106[@xmath107 in fig .
5(a ) and @xmath104[@xmath105 in fig .
5(b ) ) it imparts equal length changes ( @xmath110  pm ) to each cavity .
since both cavities respond in same way , the relative length change between two cavities is nearly zero .
therefore , a tracking mechanism would measure a correct drift caused by vertical acceleration . * when the acceleration is in @xmath21-direction ( @xmath106[@xmath107 in fig .
5(b ) and @xmath104[@xmath105 in fig .
5(c ) ) each mirror is displaced along the cavity axis by same amount without causing any relative drift i.e. , @xmath111 .     and the middle mounting hole - pair is denoted by @xmath112 . the concentric through - holes ( dia @xmath113 1 mm ) are shown to provide clearance for the vertical suspension with wires .
]    v    from the proceeding analysis we conclude that proposed fp device is immune to differential length changes if the acceleration noise is oriented along the cavity length and in the vertical direction .
however , the noise coupling in @xmath83-direction can induced non - zero differential length changes ( @xmath114 ) between two cavities . in such case
, the instantaneous drift measured by laser - lock setup will not be same as the drift of astro - cavity . based on the type of noise
some strategies can be adopted to overcome this limitation .
for example , a typical noise spikes ( an impulse excitation ) only lasts for a fraction of second and temporarily deforms the elastic cavity to produce uncorrelated tracking error in drift measurement .
after the noise pulse has elapsed the elastic cavity would return to its original shape and laser - lock would resume the correct measurements . in such case ,
the actual drift has to be estimated from a series of measurements by carefully discarding the erroneous data .
if the source of noise is constant then setup has to be properly shielded .
for example , the measurement error due to differential length changes in the proposed device can be reduced below cms@xmath2 level by attenuating the noise by @xmath115 .
the vibration sensitivity of a cavity can be greatly reduced by mounting it vertically . in this scheme the cavity axis is aligned along gravity and fp is held close to the spacer mid - plane as illustrated in fig .
the major advantage of this configuration comes from symmetric distribution of mass above and below the cavity mid - plane .
the top and the bottom mirrors are equally displaced by force of gravity in downward direction thus resulting in near - zero change in cavity length .
a mounting mechanism can be devised to hold the cavity at suitable supporting positions -also called airy - points .
the cavity spacing can be made invariant by optimizing the size , location and depth of the supporting points .
mm and the optimized hole depth was found to be hd=4.76 mm .
( d)-(f ) mirror displacement at different diameter @xmath112 of the central hole pair . in these simulations
@xmath116 mm and hd=4.76 mm and the optimized diameter @xmath112 of the central hole pair was found to be 7.34  mm . ]
three pairs of mounting holes are partially drilled to provide discrete support - points at the under - side of the rectangular spacer .
these holes are drilled up to the cavity mid - plane as shown in fig .
. the flat bottomed - surface of each hole provides a smooth contact area for cavity to rest on mounting posts from below .
alternatively , the fp can be hanged in air by using suspension wires from the above .
the retaining disks attached to each end of the wire and smoothly inserted into the mounting holes will wear the cavity weight uniformly .
the concentric through - holes ( dia @xmath113 1 mm ) are only shown to provide necessary clearance for wires .
the diameter of the mounting holes at the corners ( @xmath117 , see fig .
6(b ) ) is chosen to be 5 mm to maintain a sufficient wall thickness from edges .      in vertical mounting airy - points
are not located exactly at the geometric mid - plane of the cavity .
this is because the compression in upper - half of the cavity is not exactly countered by stretching in lower - half .
in addition there is a weight imbalance around the supporting points .
the material removed to form mounting holes reduces the weight of the lower - half of the spacer block . in order to correctly locate vertical position of the airy points , we first varied the mounting hole depth ( hd ) in our simulations and determined the gravity induced displacement of top and bottom mirrors .
the cavity length change is inferred from the displacement of probe points located on the inner surfaces of mirrors as discussed in previous section .     which corresponds to fractional length change @xmath118 below @xmath119 .
]    three test cases with different displacement trends are shown in fig .
7(a)-(c ) .
when the supporting surface is below the cavity - mid plane , e.g. hd = 4 mm in fig .
7(a ) , the vertical displacement of the upper mirrors is higher as seen from corresponding @xmath83- and @xmath22-cuts .
the situation is reversed in fig .
7(c ) where the supporting location hd = 6 mm was moved above the cavity - mid plane .
the null point occurred at hd = 4.76 mm in fig .
7(b ) for which the displacement of top and bottom mirrors become nearly equal and the cavity length becomes invariant .
figure 8 shows the rms separation between the top and the bottom cavity mirrors at various hole depths .
the fractional length change at optimal hole - depth for both cavities is below @xmath120 .
if we assume manufacturing tolerance of @xmath121 m , the rms difference still remains below @xmath122 pm , indicating the robustness of the design against machining errors .
this tolerance limit is well within the current technology of precision manufacturing .     in vertical mounting scheme with acceleration in ( a ) @xmath92-plane , ( b ) @xmath82-plane and ( c ) @xmath93-plane .
coordinate axis are visually depicted in the inset images .
for each case the orientation angle @xmath94 was swept from @xmath95 in step of @xmath96 .
( d ) the differential length change @xmath97 between two fp cavities as a function of orientation angle @xmath94 .
the @xmath22-scale to the right gives corresponding radial velocity error ( cms@xmath2 ) . ]
shape of the profile along @xmath83- and @xmath22-directions in fig .
7(a)-(c ) also indicates that displacement is not uniform across the cavity cross - section .
a linear @xmath83-tilt is added onto a radially symmetric sag .
this tilt is caused by bending of spacer in the middle that bears larger mass concentration .
we enlarged the diameter @xmath112 to remove extra material for balancing and equalizing the weight distribution in middle part of the spacer .
figure 7(e)-(f ) demonstrates how the tilt component ( x - cut ) changes sign when @xmath112 was varied from 6 mm to 9 mm . at optimum diameter ( @xmath123 mm in fig .
7 ( e ) ) , the tilt component was largely eliminated and the rms difference between @xmath83 and @xmath22-radial profile was reduced below @xmath124 m.    having found an equilibrium plane for support , we also estimated the direction dependent noise sensitivity of the vertical mount . for fea simulations ,
the mechanical displacement of the mounting holes was constrained and the fp geometry was subjected to a sweeping acceleration vector in different planes .
the calculated cavity length variations @xmath88 and differential length change @xmath98 is plotted in fig .
9 . for acceleration
along @xmath83 , the vibration sensitivity to differential length change has reduced by a factor of seven compared to horizontally mounted case . for acceleration along @xmath22 and @xmath21-direction
, the rv error is below 1 cms@xmath2g@xmath2 .
since the vibration sensitivity is significantly reduced in vertical mounting scheme , the requirements for external noise isolation and damping are also relaxed .
thermally induced material expansion or contraction is also a major source of noise against the dimensional stability of optical materials .
a desired length stability of fp cavity can be assured by choosing appropriate spacer / mirror material with low thermal expansion coefficients . to minimize the impact of temperature fluctuations , the reference cavities are also operated near their zero - crossing temperature @xcite .
the entire setup is enclosed inside an evacuated and thermally shielded environment . in this section
we investigate effect of temperature fluctuations on our fp cavities .
to first approximation , the cavity length changes in fp1 and fp2 will strongly correlate if temperature variations within the enclosure are uniform . in that case
the tracking mechanism would expectedly register a correct drift .
however , small temperature gradients inside the enclosure will expose two cavities to a locally different environment , producing dissimilar length changes in fp1 and fp2 .
this would contaminate the tracking signal with intractable radial velocity errors .
was varied from 0 mk to 100 mk . ]
to analyze the temperature sensitivity of the fp we again used fea to solved the thermal - stress problem with comsol multiphysics . in our thermal model , realistic temperature gradients were simulated by imposing appropriate temperature conditions on selected boundaries of the spacer block .
a steady - state numerical solution then establishes a linearly varying temperature and stress fields between the selected sides .
temperature sensitivity also depends on the choice of the spacer and mirror material . in these studies ,
corning ule was used as cavity spacer while three glass - ceramics namely , zerodur , silica and corning ule were tested for mirror substrates .
two worse case scenario producing unequal cavity length changes are shown in fig .
10 . the temperature sensitivities
( @xmath125 per mk ) of the device can be determined from the slopes of the plots . among the chosen spacer - mirror pairs , the ule - fs cavity has the largest sensitivity @xmath126  mk  @xmath2 ( rv error @xmath113 3.81  cms  @xmath2  mk  @xmath2 ) while ule - ule has the least sensitivity @xmath127  mk  @xmath2 ( rv error @xmath113 0.27 cms@xmath2mk@xmath2 ) for temperature gradient shown along spacer edge in fig .
as see in fig .
10(b ) , the temperature sensitivity is almost halved for the gradient considered along diagonal direction .
the contribution of cte to cavity length change from mirror substrate is substantially large compared to the contribution from spacer material .
the main drawback of fused silica is its large cte which places stringent requirements on the design of thermal control system of the cavity .
for example , to keep the rv error of dual cavity fp below 1  cms@xmath2 level , a temperature uniformity at sub - mk level has to be ensured .
this requirement is relaxed to few - mk if both spacer - mirror pair are made with extremely low expansion materials such as zerodur and ule .
a new method of building temperature insensitive fs cavity mirrors is expected improve their performance in precision metrology @xcite . despite large cte
, fused silica is considered as best performing material for cavity mirrors .
this is because fs has very small mechanical losses at room temperature which helps to lower cavity s thermal noise floor which is about 2 - 3 times better compared to the ule or zerodur @xcite .      in ideal case of equal length cavities
i.e @xmath128 , temperature fluctuations will not produce any relative length change between two cavities . in practice
, it is not possible to manufacture cavities of equal length with arbitrary precision .
a small offset in lengths may arise from machining errors and/or imperfections in the optical bonding between the mirrors and the spacer . to study this effect we created cavities with different length offset and calculated their differential expansion from the thermal model .
the length offset was introduced by shortening the spacer thickness @xmath129 of fp1 without altering @xmath130 . in the fea simulations the temperature of the entire fp block
was varied uniformly in step of 1 mk .
figure 11 shows thermally induced differential length changes between two cavities for different length offsets .
clearly , a 0.1 k difference from the temperature set point can produce @xmath131 pm differential expansion ( the green curve in fig .
11 ) between two cavities if the initial length offset was @xmath132 m
. this amounts to @xmath133  cms@xmath2 rv tracking error .
these calculations reinforces the need of high thermal stability ( few mk ) and cavity length matching better than @xmath121 m for keeping the influence of temperature changes on rv error below 1  cms@xmath2 .
k to @xmath134 + 0.1 k in incremental step of 1 mk .
the spacer material used in simulation is ule and mirror material is fused silica .
various length offsets @xmath135 are given in @xmath58 m . ]
in fea , the numerical solution of a problem is approximated by partitioning the cad geometry into smaller parts .
the accuracy of the numerical solution is intrinsically linked to the mesh size . by refining the mesh ,
the solution becomes more accurate but this is achieved at the cost of finite computational resources and time . in a realistic fea model , the difference between exact and approximate solution should be minimized and error should not exceed the tolerance limit defined by some accepted criteria .
we used the default tetrahedral mesh settings with adaptive mesh refinement in our fea software to discretize the cavity model .
the refinement algorithm globally adjusts the mesh size to provide sufficient resolution around all regions of interest within the geometry .
there is no generic rule to verify the numerical accuracy of fea model -especially in complex problems where exact solution is not known a priory . in order to evaluate the discretization error due to finite mesh size
, we compare the relative displacement @xmath74 of fp1 and fp2 mirrors for a simple gravity loading case ( see fig .
4 ) already discussed in section 4.1 . from the application of uniform load and symmetry consideration , exact solution requires the cavity displacement in fp1 and fp2 to be equal .
any difference in displacement can be attributed to numerical inaccuracies originating from meshing errors .
figure 12 shows the displacement data ( @xmath83-cut ) plotted for input cavity mirrors at three different mesh densities .
we note that compared to the smaller mesh sizes , the magnitude of the displacement is slightly underestimated at _ ` coarse mesh'_. in addition , a small relative difference can also be seen in fig .
this difference reduces progressively with decreasing mesh size ( see fig .
12(b ) ) and becomes insignificant for ` _ _ extra - fine mesh _ _ ' shown in fig.12(c ) where two curves completely overlap . in each case
the numerical error is estimated from the rms difference between the two curves .
the discretization rms error for extra - fine mesh is @xmath136 m. this amounts to @xmath137 error for cavity length stability targeted at @xmath122 m level .
all fea simulations in this paper were carried out with extra - fine mesh .
another possible source of error is uncertainties in thermal and mechanical properties of the ceramic material listed in table 1 .
typically , the material parameters can have about @xmath138 deviation from their nominal values . we varied the material parameters in our simulation by @xmath139 but did not find any noticeable difference in vibration or temperature sensitivities .
the calculated deviation was of the same order as the discretization errors .
this is not unexpected since our fp is lightweight and cavity length is short . for the same level of uncertainties in poisson ratio ,
youngs modulus and material density , a vibration sensitivity of @xmath140 was achieved in longer cavities @xcite .
-data probes on the input mirrors of fp1 and fp2 for which the cavity displacement is numerically computed
. the default element size ( min , max ) for each case in also indicated . the green curve ( @xmath22-scale to the right )
shows the difference between the two plots . the rms error for coarse , fine and extra - fine
meshing is @xmath141 m , @xmath142 m and @xmath143 m , respectively . ]
the doppler survey for many astronomy programmes requires high resolution spectrograph with exceptional long term stability in rv precision . in this paper
we have proposed a double cavity fp design for aiding high precision doppler measurements in astronomy .
the proposed device has two parallel optical cavities embedded inside a 10  mm thick monolithic ule block .
a low - finesse _ astro - cavity _ is meant to generate multiline reference spectra for spectrograph calibration while a high finesse _ lock - cavity _ is needed to provide steep discriminant signal for tracking the cavity drift .
external noise can have severe detrimental effect on the performance and reliability of the tracking measurements .
an accurate measurement of drift is possible only if the noise interventions are common to both cavities .
we have used fea to analyze fp s sensitivity to mechanical vibrations and temperature fluctuations .
a 3d fea simulations were carried out for the fp structures comprising five domains ( 4 mirrors + 1 spacer ) . an optimized cavity design with proper
mounting and support system helps insulate the noise .
two mounting geometries were investigated for tracking errors resulting from differential length changes between the cavities .
the horizontal geometry showed vibrational sensitivity up to @xmath144@xmath109 ( rv errors @xmath145 cms@xmath2 ) for acceleration directed along the longer edge ( @xmath83-axis ) of the block .
the noise along @xmath22- and @xmath21-direction is common and does not contribute to differential length change . for achieving 1 cms@xmath2 rv stability of the device the horizontal mounting configuration may require a well thought strategy to insulate the vibration noise down to @xmath146 .
cavity lengths can be made invariant by mounting the fp vertically .
we worked out a vibration insensitive design by numerically optimizing the hole - depth and diameter of the mounting holes . in the optimized design
the tilt - component was eliminated and a clear null point was achieved for gravity induced displacement between the mirrors .
we also studied the impact of thermal instability and choice of mirror / spacer material on the tracking performance .
a steady - state thermal analysis was carried out to evaluate the effects of thermal gradient .
numerical simulations show that temperature gradient within the fp enclosure has to be accurately controlled and relative length offset between two cavities should be minimized to avoid differential length changes that may cause spurious drift and unintended tracking errors .
the proposed dual fp device can be easily fabricated with available machining precision of @xmath147 m and readily achievable thermal stability of few mk .
this work is also relevant to cavities used in laboratory experiments designed to test lorentz invariance where a pair of orthogonal fp cavities is employed to detect anisotropy in the speed of light @xcite .
the entire optical setup is put on a rotating table . a frequencies beat - note between two lasers , each locked to one of the two orthogonal cavities ,
is examined while the entire setup is continuously rotated on the turntable .
any possible anisotropy is expected to modulate the laser beat signal at twice the rotation rate .
however , the discrimination of the anisotropic signal from competing noise artifacts is not trivial . even though the length of two cavities is defined by a single ule block
, detection ambiguities may arise , e.g. , from periodic deformation of cavity length caused by gravitational and centrifugal forces .
the approach outlined in our studies might be useful to understand and model various systematics producing relative drift in cross - cavities .
rkb acknowledges the research support for this project ( reference no .
emr/2014/000941 ) from science and engineering research board ( serb ) , department of science and technology ( dst ) , india . a.r .
acknowledges the research support from european research council under the fp7 starting grant agreement number 279347 and funding from dfg grant re 1664/9 - 1 .
huke , p. , holzhuter , h. , & reiners , a. [ 2015 ] `` transfering rb+ hyperfine - structure stability to a fabry - perot resonator used as a frequency standard for astronomical spectrographs , '' in _ proc .
spie _ * 9526*,98260i-7 .
probst , r. a. , lo curto , g. , avila , g. , _ et al . _
[ 2014 ] `` a laser frequency comb featuring sub - cm / s precision for routine operation on harps , '' in _ spie astronomical telescopes & instrumentation , _ pp .
91471c-91471c .
schfer , s. , & reiners , a. [ 2012 ] `` two fabry - perot interferometers for high precision wavelength calibration in the near - infrared , '' in _ spie astronomical telescopes & instrumentation _ , pp . 844694 - 844694 .