document: tidal interactions play an important role in dynamical processes in the two - body problem in close binary systems : star  star ( binary star ) or star planet . 
   they can lead to such phenomena as synchronization and orbital circularization ( hut 1981 ; zahn 1977 ) as well as to the tidal capture ( press and teukolsky 1977 ) or disruption of an object ( star ) ( ivanov and novikov 2001 ) and the fall of the object onto the star ( rasio et al . 
   1996 ; penev et al . 2012 ; bolmont and mathis 2016 ) . in this paper 
   we consider the tidal interaction of two bodies : a star and a point source . 
   the point source can be both a star ( a neutron star , a white dwarf , etc . ) and a planet . 
   below we will call the point source a planet by implying that this can also be a star . 
   since the evolution time scales of the eccentricity or semimajor axis strongly depend on the orbital period of the binary system and for some stellar models can take values up to 108 yr or more for periods of about 5 days ( ivanov et al . 2013 ; chernov et al . 2013 ) , the evolution of the star itself should be taken into account on such time scales . 
   as the star evolves , the orbital parameters change due to tidal interactions . in this paper 
   we investigate the dynamical tides by taking into account the stellar evolution . for our study 
   we chose three types of stars with masses of one , one and a half , and two solar masses . 
   we consider all stars without allowance for their rotation and magnetic field and touch on the stellar physics itself superficially , as far as this problem requires . 
   the star of one solar mass at lifetimes @xmath3 yr closely corresponds to our sun and has a radiative core and a convective envelope on the main sequence . 
   the other two stars of one and a half and two solar masses are more massive and have a more complex structure . 
   these stars have a convective core and a radiative envelope on the main sequence ( a more precise structure is presented below ) . 
   the problem of determining the tidal evolution is reduced to the problem of determining the normal modes of stellar perturbations and to calculating the energy and angular momentum exchange in the star.planet system ( ivanov and papaloizou 2004 , 2010 ; papaloizou and ivanov 2010 ; lanza and mathis 2016 ) . 
   the low - frequency g - modes of the stellar oscillations play an important role in the theory of dynamical tides . 
   the tidal interactions are fairly intense at small periastron distances of the planet . for a periastron distance @xmath4 au , the dimensionless excitation frequency is @xmath5 , which corresponds to g - modes . 
   a large number of exoplanets in stellar systems have been discovered in the last few years owing to the kepler , superwasp , and other observational programs . in particular , 
   short - period massive planets with an orbital period of a few days , the so - called hot jupiters , have been detected ( winn 2015 ) . as a rule , 
   the hot jupiters have low eccentricities , which points to the importance of tidal interactions ( ogilvie 2014 ) . 
   the results of this paper can be directly applied to some of such systems . 
   for example , the system ybp1194 is a solar twin ( brucalassi et al . 
   a planet with a mass of 0.34@xmath6 and a period of only 6.9 days revolves around this star . 
   for such a short - period planet the dynamical tides must be fairly intense and must affect the orbital evolution . 
   predictions about the subsequent evolution of this planet can be made by analyzing this system . 
   one of the results of this evolution is the fall of the planet onto the star . 
   the possibility of such a fall has been considered in many papers ( see , e.g. , rasio et al . 
   1996 ; penev et al . 2012 
   ; weinberg et al . 2012 ; essick and weinberg 2016 ) . rasio et al . 
   ( 1996 ) considered the possibility of the fall of the planet onto the star due to quasi - static tides and provided a plot for solar - like stars that shows the threshold , as a function of planetary mass and orbital period , below which the planet falls onto the star . 
   penev et al . 
   ( 2012 ) considered tides with a constant tidal quality factor @xmath7 specified phenomenologically . in reality , this factor will depend on the planet s orbital period ( ivanov et al . 
   2013 ) and stellar age . 
   essick and weinberg ( 2016 ) took into account the energy dissipation due to nonlinear interaction of modes with one another . 
   in contrast to our approach ( ivanov et al . 
   2013 ; chernov et al . 
   ( 2013 ) , the simultaneous solution of a large number of ordinary differential equations for each stellar model is suggested , with only solar - type stars having been considered . 
   in this paper we consider the evolution of stars with masses of one , one and a half , and two solar masses . 
   data on the stars are presented in tables 1 3 . 
   a novelty of this study is a consistent allowance for the stellar evolution . for each moment of the star 
   s lifetime we calculated the spectra of normal modes , the overlap integrals ( press and teukolsky 1977 ) , which are a measure of the intensity of tidal interactions ( for a generalization to the case of a rotating star , see papaloizou and ivanov 2005 ; ivanov and papaloizou 2007 ) , and the time scales of the orbital parameters ( the tidal quality factor @xmath7 ) . 
   the overlap integrals are directly related to the tidal resonance coefficients that were introduced by cowling ( 1941 ) ( see also rocca 1987 ; zahn 1970 ) . 
   all of the quantities marked by a tilde are dimensionless ; the normalization is presented in appendix a. 
   such a quantity as the overlap integral q is of great importance in the theory of dynamical tides . 
   it is specified by the expression ( press and teukolsky 1977 ; zahn 1970 ) @xmath8 where m is the stellar mass , r is the stellar radius , and @xmath9 is the dimensionless overlap integral ; the remaining quantities and their dimensionless forms are defined in appendix a. the overlap integral serves as a measure of the intensity of the tidal forces and plays a crucial role in the physics of dynamical tides ( press and teukolsky 1977 ; zahn 1970 , 1975 , 1977 ) . here 
   , we will consider only the quadrupole part of the tidal forces @xmath10 , because it is this part that makes the greatest contribution to the overlap integral ( press and teukolsky 1977 ) . 
   the change in orbital semimajor axis and eccentricity with time is determined from the energy and angularmomentum conservation laws and is specified by the following formulas ( ivanov et al . 
   2013 ) : @xmath11 here , @xmath12 is the orbital eccentricity , @xmath13 is the semimajor axis , and the time scales of the change in semimajor axis @xmath14 and eccentricity @xmath15 in the case of low eccentricities ( @xmath16 ) are specified by the relations ( ivanov et al . 
   2013 ) @xmath17 where @xmath18 is the planetary mass . 
   the rate of change of energy @xmath19 , @xmath20 in the case of a dense spectrum and moderately large viscosities is specified by eqs . 
   ( 51 ) from ivanov et al . 
   the spectrum is deemed to be dense enough if @xmath21 . for simplicity 
   , we will consider here only the time scales of the changes in semimajor axis . 
   using formulas from subsection 7.1 of the paper by ivanov et al . 
   ( 2013 ) , we obtain @xmath22 where @xmath23 is the ratio of the planetary mass to the stellar mass , @xmath24 is the difference between adjacent frequencies , @xmath25 is the frequency number , and @xmath26 ( ivanov et al . 
   . the orbital period of the planet around the star is specified by the relation @xmath27 below , using eqs . 
   ( [ ta ] ) and ( [ porb ] ) and solving eq . 
   ( [ tate ] ) , we will present the change in the planet s orbital period as a function of time by taking into account the stellar evolution .    as is clear from eq . 
   ( [ overlapsintegral ] ) , the eigenfunctions and , accordingly , eigenfrequencies of the stellar oscillations should be known to calculate the overlap integral . 
   there exist three types of eigenfrequencies for a nonrotating star : the so - called p- , f- , and g - modes . 
   we will be concerned only with the low - frequency g- modes . 
   the properties of the g - modes are determined by the brunt - v@xmath28is@xmath29l@xmath29 frequency . 
   it is specified as follows ( christensen - dalsgaard 1998 ) : @xmath30 where @xmath31 is the pressure , @xmath32 is the density , @xmath33 is the gravitational acceleration , @xmath34 is the adiabatic index . 
   however , this definition of the brunt - v@xmath28is@xmath29l@xmath29 frequency is difficult to apply to many stars . 
   this is related to the numerical errors , to the calculation of the derivative of the density , and to the fact that there exist regions in stars where both terms in parentheses can be of the same order of magnitude . 
   therefore , a different definition is used ( brassard et al . 1991 ) : @xmath35,\end{aligned}\ ] ] where @xmath36 the temperature gradient @xmath37 @xmath38 is the mass fraction of atoms of type i , and @xmath39 it is convenient to divide the brunt - v@xmath28is@xmath29l@xmath29 frequency into two parts : the structural part ( structure terms ) brunt - v@xmath28is@xmath29l@xmath29 @xmath40\end{aligned}\ ] ] and the compositional part @xmath41 the compositional part is related to the change in stellar chemical composition due to nuclear reactions , and it is of great importance at the boundary of the convective and radiative regions . 
   another important stellar characteristic is the acoustic frequency defined by the formula @xmath42    the low - frequency g - modes are excited at @xmath43 . 
   the high - frequency p - modes are excited at @xmath44 . 
   we will not consider the p - modes , because they play a secondary role in the theory of dynamical tides . 
   thus , the problem is reduced to finding the eigenfunctions and eigenfrequencies of the stellar oscillations , calculating the overlap integrals , and calculating the characteristic orbital parameters as a function of time . 
   the equations that describe the eigenfrequencies and eigenfunctions of a star in the adiabatic and cowling approximations are presented in appendix a. the stars were modeled with the mesa software package ( paxton et al . 
   2011 , 2013 , 2015 ) ; the data obtained were then interpolated by the steffen ( 1990 ) method to two million points . using these data 
   , we solved the eigenfrequency and eigenfunction problem and calculated the overlap integrals . 
   the fourth - order runge 
   kutta method was used to solve the differential equations . 
   the data for the stars and the overlap integrals were compared with those from ivanov et al . 
   ( 2013 ) and chernov et al . 
   in this section we present the results of our calculations of the change in the planet s orbital period for three types of stars with masses of one , one and a half , and two solar masses . 
   the course of stellar evolution for these models is shown on the hertzsprung russell diagram depicted in fig . 
   data for each star are presented in the tables in appendix b. the main phase of stellar evolution is the main - sequence phase . 
   it is at this phase that any star spends most of its lifetime . 
   this phase is indicated in fig . 1 by the nearly horizontal straight line .          in this section we present the results of our calculations of the overlap integrals for various lifetimes of the star of one solar mass and the change in the planet s orbital period with time . 
   the star was modeled from @xmath45 to @xmath46 yr . at 
   times @xmath47 yr this star closely corresponds to our sun , and below we will call this star the sun . at the initial time the sun was completely convective , hydrogen burning has not yet begun , and its metallicity is z = 0.02 . 
   figure 2 shows a dimensionless dependence of the solar density on radius at the initial and subsequent times . 
   is@xmath29l@xmath29 frequency ( b ) for the sun versus radius for various times : the dashed , dotted , solid , and dash  dotted lines are for the ages @xmath45 , @xmath48 , @xmath49 , and @xmath46 yr , respectively.,title="fig:",scaledwidth=45.0% ] is@xmath29l@xmath29 frequency ( b ) for the sun versus radius for various times : the dashed , dotted , solid , and dash  dotted lines are for the ages @xmath45 , @xmath48 , @xmath49 , and @xmath46 yr , respectively.,title="fig:",scaledwidth=45.0% ]    is@xmath29l@xmath29 frequency versus radius for the solar age @xmath49 : the solid , dashed , and dash - dotted lines indicate the total brunt - v@xmath28is@xmath29l@xmath29 frequency @xmath50 , the structural part @xmath51 , and the compositional part @xmath52 , respectively.,scaledwidth=45.0% ]    the course of stellar evolution and the excitation of eigenmodes can be described as follows . at times 
   @xmath53 yr the sun was completely convective ( @xmath54 ) and , therefore , no low - frequency g - modes were excited , but only the high - frequency f- and p - modes were excited . subsequently , 
   as the sun evolves ( at times @xmath55 yr ) a radiative region ( radiative core ) where the brunt - v@xmath28is@xmath29l@xmath29 frequency is greater than zero ( @xmath56 ) appears and , consequently , the excitation of low - frequency g - modes begins . 
   only the structural part makes a major contribution to the brunt - v@xmath28is@xmath29l@xmath29 frequency ; the compositional part is zero . 
   this is because at such lifetimes of the sun hydrogen burning makes a minor contribution to the energetics and evolution dynamics of the sun . 
   in the course of subsequent evolution the radiative region expands , the number of lowfrequency g - modes increases , and their spectrum becomes dense . 
   the low - frequency g - modes begin to play an increasingly important role . at times 
   @xmath57 yr hydrogen burning begins to make an increasingly large contribution to the total energy of the star and , consequently , the compositional part of the brunt - v@xmath28is@xmath29l@xmath29 frequency also begins to contribute to the total frequency . 
   in the course of subsequent evolution , as hydrogen burns , this contribution becomes progressively larger . 
   one of the problems being investigated here is to take into account the evolution of the brunt - v@xmath28is@xmath29l@xmath29 frequency and the influence of this evolution on the tidal interactions . 
   the radiative region expands approximately to a radius @xmath58 . for the present sun 
   the boundary between the radiative and convective regions is @xmath59 . in fig . 
   2 the brunt - v@xmath28is@xmath29l@xmath29 frequency is plotted against radius for various stellar evolution times . 
   it can be seen from the plot that the brunt - v@xmath28is@xmath29l@xmath29 frequency is subject to significant changes as the solar age increases ; both the frequency spectrum and its density change accordingly , which can directly affect the tidal forces . in fig . 
   3 the brunt - v@xmath28is@xmath29l@xmath29 frequency is plotted against radius for the present sun . 
   it follows from this plot that the brunt - v@xmath28is@xmath29l@xmath29 frequency exhibits a double - humped structure . 
   the first and second humps are associated with the compositional and structural parts of the brunt - v@xmath28is@xmath29l@xmath29 frequency , respectively . at a solar age 
   @xmath60 yr the two contributions become comparable .    the overlap integrals ( 1 ) are plotted against frequency in the cowling approximation in fig . 
   4 . for a solar age @xmath53 yr , 
   as has been said above , only the f- and p - modes exist ; therefore , there is no low - frequency part of the spectrum . 
   it can be seen from the plot that the absolute value of the overlap integrals decreases in the low - frequency part of the spectrum as the sun evolves . 
   this decrease is gradual , because the stellar structure does not change in the evolution time of the sun . on the main sequence 
   the sun always has a radiative core and a convective envelope . 
   in the high - frequency part of the spectrum , as the stellar age increases , the absolute value of the overlap integrals increases and a characteristic structure appears in the spectrumin the form of kinks .    , @xmath48 , @xmath49 , and @xmath46 yr , respectively.,scaledwidth=45.0% ]    figure 5 plots the change in the orbital period ( in days ) of a planet with a mass of one jupiter mass around the sun with time ( in years ) . the solar evolution models for which the overlap integrals were calculated are presented in table 1 . at lifetimes of the sun @xmath61 
   yr the spectrum at frequencies @xmath62 is insufficiently dense . at other lifetimes of the sun 
   the spectrum may be deemed sufficiently dense with a good accuracy for all frequencies . 
   therefore , eq . ( 4 ) may be used . 
   equation ( 2 ) was integrated over time from @xmath63 to @xmath64 yr . 
   the initial orbital periods of the planet around the sun were taken to be @xmath65 days ( solid curves ) . 
   it follows from fig . 
   5 that in the main - sequence lifetime of the sun a planet with a mass of one jupiter mass with the initial orbital period @xmath66 days will reduce its orbital period to @xmath67 days due to dynamical tides . 
   the time scale of the change in semimajor axis for this orbital period @xmath67 days and this solar age @xmath64 yr is @xmath68 yr . 
   consequently , it can be said with confidence that this planet will fall onto the sun within several tens of millions of years . 
   the same analysis showed that a planet with the initial orbital period @xmath69 days would reduce its orbital period to @xmath70 days and would fall onto the sun only within approximately half a billion years . 
   for the initial orbital period @xmath69 days fig . 5 compares the changes in orbital period with allowance for the solar evolution ( solid curve ) , without allowance for the solar evolution ( dash  dotted curve ) , and the one calculated using eq . 
   ( 25 ) from essick and weinberg ( 2016 ) ( dotted curve ) . to calculate the change in orbital period without allowance for the solar evolution 
   , we took the time scale at the solar age @xmath0 yr . 
   all curves qualitatively coincide , and the influence of solar evolution may be disregarded in the estimation . 
   this is because hydrogen burning in the sun is gradual , and no change in structure occurs . 
   quantitatively , the dotted curve from essick and weinberg ( 2016 ) gives smaller changes in the planet s orbital period . 
   this may be related to both the neglect of the change in the time scales of the orbital parameters as the sun evolves and to a different approach in this paper to the problem of the energy dissipation of tidally excited modes associated with their nonlinear interactions .          in this subsection 
   we present the results of our calculations of the overlap integrals and the change in orbital period for the star with a mass of one and a half solar masses . 
   the star is considered from the time the protostar begins to collapse . 
   this time is quite arbitrary and is @xmath72 yr ; the initial radius is @xmath73 . 
   the initial metallicity of the star is @xmath74 . at lifetimes of the star up to @xmath75 yr it is completely convective ( @xmath54 ) . in such a star no lowfrequency g - modes 
   are excited , and only the highfrequency f- and p - modes are present ; therefore , the effective energy and angular momentum transfer from the orbit to the star will be suppressed . 
   as the stellar age increases further , hydrogen burning begins at the stellar center and a radiative core appears ( @xmath56 ) , which begins to expand , encompassing progressively newer layers . 
   the star has a well - defined convective envelope and a radiative core . 
   the radiative region expands up to @xmath76 for the age @xmath77 yr . a fairly dense spectrum of low - frequency g - modes appears very rapidly in such a star 
   .    at times @xmath78 yr the stellar structure changes completely . 
   a convective core appears in the star , while its envelope becomes completely radiative . 
   the star abruptly contracts from a radius @xmath79 to @xmath80 . 
   the evolution time scales of the orbital parameters ( 2 ) increase sharply . 
   the star passes from a solar - type star that has a radiative core and a convective envelope to a star with a convective core and a radiative envelope . 
   is@xmath29l@xmath29 frequency ( b ) for the star of one and a half solar masses versus radius for various times : the dashed , dotted , solid , and dash  dotted lines are for the ages @xmath81 , @xmath82 , @xmath83 , and @xmath84 yr , respectively.,title="fig:",scaledwidth=45.0% ] is@xmath29l@xmath29 frequency ( b ) for the star of one and a half solar masses versus radius for various times : the dashed , dotted , solid , and dash  dotted lines are for the ages @xmath81 , @xmath82 , @xmath83 , and @xmath84 yr , respectively.,title="fig:",scaledwidth=45.0% ]    in fig . 
   6 the stellar density and brunt - v@xmath28is@xmath29l@xmath29 frequency are plotted against radius for various times . from the figure 
   we see how the central density increased as the star contracted . 
   it can be seen from the figure for the brunt - v@xmath28is@xmath29l@xmath29 frequency that the star was initially completely convective ( dashed curve ) and subsequently a solar - type one ( the envelope was convective , and the core was radiative ; the dotted curve ) . at the subsequent times 
   the star has a more complex structure , the core becomes convective , and the envelope is completely radiative ( solid curve ) . 
   a thin convective envelope then again appears in its envelope ( dash  dotted curve ) . 
   is@xmath29l@xmath29 frequency versus radius for the star of one and a half solar masses for the age @xmath84 yr : the solid , dashed , and dash - dotted lines indicate the total brunt - v@xmath28is@xmath29l@xmath29 frequency @xmath50 , the structural part @xmath51 , and the compositional part @xmath52 , respectively.,scaledwidth=45.0% ]    figure 7 shows the brunt - v@xmath28is@xmath29l@xmath29 frequency as well as its structural and compositional parts . it follows from the figure that the compositional part makes the greatest contribution at the boundary between the convective and radiative regions .    , @xmath82 , @xmath83 , and @xmath84 yr , respectively.,scaledwidth=48.0% ]    the overlap integrals are plotted against frequency in fig . 
   . the dashed line does not extend to the lowfrequency part of the spectrum , because the star at this time is completely convective ( @xmath54 ) . 
   it can also be seen from the figure that the absolute value of the overlap integrals rapidly drops by more than two orders of magnitude as the star s lifetime increases . 
   such catastrophic changes are due to the change in stellar structure . 
   in the high - frequency part the absolute value of the overlap integrals increases by an order of magnitude , and characteristic kinks appear in the spectrum . 
   figure 9 plots the change in the orbital period ( in days ) of a planet with a mass of one jupiter mass around the star of one and a half solar masses with time ( in years ) . 
   the stellar evolutionmodels for which the overlap integrals were calculated are presented in table 2 . 
   equation ( 2 ) was integrated over time from @xmath85 to @xmath86 yr . 
   for such a stellar lifetime interval the g - mode spectrum is sufficiently dense . 
   the initial orbital periods of the planet around the sun are @xmath87 days ( solid curves ) . 
   it follows from the figure that all planets with a mass equal to one jupiter mass and with an initial orbital period of 2 days or shorter will fall onto the star in the main - sequence lifetime of the star . 
   a planet with the initial orbital period @xmath88 days will fall onto the star , while a planet with the initial orbital period @xmath89 days will reduce its orbital period to @xmath90 days in the main - sequence lifetime of the star . 
   the dashed curve indicates the change in orbital period without allowance for the stellar evolution . 
   the time scale of the orbital change was taken at the star s lifetime @xmath91 yr . 
   it can be seen from fig . 
   9 that the planet does not fall onto the star in the mainsequence lifetime of the star without allowance for its evolution . 
   the last star that we will consider is the star with a mass of two solar masses . 
   the time the protostar begins to collapse is @xmath92 yr , the initial stellar radius is @xmath93 , and the stellar metallicity is @xmath74 . at initial times up to @xmath94 
   yr the star is completely convective ( @xmath54 ) . 
   no low - frequency g - modes are excited ; only the high - frequency f- and p - modes are present . 
   the dynamical tides play no significant role in such a star . 
   however , as the star evolves further , at times @xmath94 yr a radiative region ( radiative core ) appears in the star , which begins to expand . 
   the low - frequency g - modes are well excited in this region ( @xmath56 ) , and the dynamical tides begin to play an increasingly significant role as the star evolves further . 
   such a star resembles solartype stars that have a radiative core and a convective envelope . 
   the radiative region expands quite rapidly , and this region almost reaches the stellar surface @xmath95 already at times @xmath96 yr . the star becomes almost completely radiative ; the brunt - v@xmath28is@xmath29l@xmath29 frequency in the entire region becomes positive ( @xmath56 ) in the entire region , except for the stellar surface . 
   a sufficiently dense g - mode spectrum is generated in such a star . in fig . 
   10 the stellar density and brunt - v@xmath28is@xmath29l@xmath29 frequency are plotted against radius for various times . 
   is@xmath29l@xmath29 frequency ( b ) for the star of two solar masses versus radius for various times : the dashed , dotted , solid , and dash  dotted lines are for the ages @xmath97 , @xmath98 , @xmath99 , and @xmath100 , respectively.,title="fig:",scaledwidth=45.0% ] is@xmath29l@xmath29 frequency ( b ) for the star of two solar masses versus radius for various times : the dashed , dotted , solid , and dash  dotted lines are for the ages @xmath97 , @xmath98 , @xmath99 , and @xmath100 , respectively.,title="fig:",scaledwidth=45.0% ]    as the lifetime increases further ( @xmath101 yr ) , a convective core ( @xmath54 ) appears in the star , which begins to expand as the star evolves . 
   the compositional part of the brunt - v@xmath28is@xmath29l@xmath29 frequency @xmath52 begins to play an increasingly large role at the boundary between the convective and radiative regions ( fig . 
   .    the stellar structure does not undergo significant changes as the lifetime increases further . 
   only the convective core evolves . 
   initially , the convective region expands , and then its expansion is replaced by contraction ( fig . 
   is@xmath29l@xmath29 frequency versus radius for the star of two solar masses for the age @xmath100 yr : the solid , dashed , and dash - dotted lines indicate the total brunt - v@xmath28is@xmath29l@xmath29 frequency @xmath50 , the structural part @xmath51 , and the compositional part @xmath52 , respectively.,scaledwidth=45.0% ]    , @xmath98 , @xmath99 , and @xmath100 yr , respectively.,scaledwidth=45.0% ]    the overlap integrals are plotted against frequency in fig . 
   . the dashed line does not extend to the low - frequency part of the spectrum , because the star at this time is completely convective ( @xmath54 ) . as the star s lifetime increases , the absolute value of the overlap integral rapidly drops by more than an order of magnitude . 
   such catastrophic changes are due to the change in stellar structure . when the stellar structure undergoes no significant changes , the overlap integrals change gradually and insignificantly ( see fig . 
   12 , the solid and dash - dotted lines ) . 
   figure 13 plots the change in the orbital period ( in days ) of a planet with a mass of one jupiter mass around the star of two solar masses with time ( in years ) . 
   the evolution models of the star of two solar masses for which the overlap integrals were calculated are presented in table 3 . 
   equation ( 2 ) was integrated over time from @xmath102 to @xmath103 yr . 
   for such a stellar lifetime interval the g - mode spectrum is sufficiently dense . 
   the initial orbital periods of the planet around the sun were taken to be @xmath104 days ( solid curves ) . 
   it follows from the figure that all planets with a mass equal to one jupiter mass and with an initial orbital period of 2 days or shorter will fall onto the star in the main - sequence lifetime of the star ( fig . 
   a planet with the initial orbital period @xmath89 days will reduce its orbital period only to @xmath105 days in the main - sequence lifetime of the star . 
   the dashed curve indicates the change in orbital period without allowance for the stellar evolution . 
   the time scale of the orbital change was taken at the star s lifetime @xmath106 yr . as can be seen from fig . 
   13 , the planet does not fall onto the star in the main - sequence lifetime of the star without allowance for its evolution . 
   q        just as for the star of one and a half solar masses , the evolution time scales ( 3 ) of the planet decrease as the star begins to move off the main sequence . 
   this is because the hydrogen reserves at the stellar center begin to be depleted , and the convective core becomes a helium one . 
   this leads to an expansion of the star and , consequently , to an increase in the number of g - modes and the density of the spectrum , which leads to an increase in the pumping of orbital energy into the energy of stellar oscillations and to the fall of short - period planets onto the star . 
   we considered the change in the orbital period of a planet around a star with a mass of one , one and a half , and two solar masses . 
   the behavior of the overlap integrals for these models as a function of the main - sequence lifetime of the star was analyzed . 
   we showed that for the same star with a different structure during its evolution the overlap integrals could differ by two or more orders of magnitude . 
   we showed that all planets with a mass of one jupiter mass that revolve around the sun with an orbital period @xmath1 days or shorter should fall onto the stellar surface in the main - sequence lifetime of the star . 
   all planets with a mass of one jupiter mass that revolve around the stars with a mass of one and a half and two solar masses with an orbital period @xmath2 or shorter will fall onto the stellar surface in the main - sequence lifetime of the star . 
   such a fall may not occur without allowance for the evolution of the star itself . 
   these results qualitatively agree with the previous results from bolmont and mathis ( 2016 ) and penev ( 2012 ) . 
   calculations for the stellar models considered in this paper should be made for a more accurate quantitative estimate . 
   i am grateful to p.b . ivanov who read the paper for a number of valuable remarks and to j. papaloizou for a fruitful discussion . 
   this work was financially supported by the russian foundation for basic research ( project nos . 
   140200831 , 150208476 , 160201043 ) , grant no . 
   nsh-6595.2016.2 from the president of the russian federation for state support of leading scientific schools , and program 7 of the presidium of the russian academy of sciences . 
   the system of equations describing the adiabatic perturbations in stars is written in dimensionless form ( christensen - dalsgaard 1998 ) : @xmath107 where the following dimensionless variables are introduced : @xmath108 is the radius , @xmath109 is the radial perturbation , @xmath110 is the horizontal perturbation , @xmath111 and @xmath112 are the perturbations of the gravitational potential and its derivative , respectively . 
   the subscript @xmath113 refers to an unperturbed quantity , and the prime @xmath114 refers to a perturbed quantity . 
   a tilde over a quantity means that this quantity is dimensionless . 
   the normalization is done as follows : @xmath115 where @xmath116 is the sound speed squared , @xmath31 is the pressure , @xmath32 is the density , and @xmath50 is the brunt - v@xmath28is@xmath29l@xmath29 frequency squared . 
   four boundary conditions should be added to the four first - order differential equations . 
   two boundary conditions are specified at the stellar center , and the other two are specified on the stellar surface ( christensen - dalsgaard 1998 ) . 
   the boundary conditions at the stellar center @xmath117 are @xmath118 the boundary conditions on the stellar surface @xmath119 are @xmath120 there exists the cowling ( 1941 ) approximation , in which the perturbed gravitational potential is neglected , to investigate the low - frequency g - modes . 
   this approximation holds good and simplifies considerably the system of equations for the adiabatic perturbations ( 14 ) . in the cowling approximation 
   the system of equations ( 14 ) will be rewritten in a simpler form : @xmath121 in the cowling approximation the system of equations ( 18 ) ceases to depend explicitly on the stellar density . 
   thus , in the low - frequency limit the radial and horizontal displacements do not depend on the stellar density . 
   in contrast , the overlap integrals ( 1 ) directly depend on the density . in the cowling approximation not four but two 
   boundary conditions should be added to the system of equations ( 18 ) . 
   one boundary condition is specified at the stellar center , and the other one is specified on the stellar surface . 
   the boundary condition at the stellar center @xmath117 is @xmath122 the boundary condition on the stellar surface @xmath123 is @xmath124 to determine the eigenfrequency spectrum of the system of equations ( 14 ) or ( 18 ) , this system should be redefined . 
   this requires that the determinant composed of the eigenfunctions of the system of equations ( 14 ) or ( 18 ) be equal to zero . 
   this will hold only for some frequencies ( eigenfrequencies ) of the system of equations ( 14 ) or ( 18 ) , which determines the perturbation spectrum . 
   the tables present the models for which the overlap integrals ( 1 ) and the evolution time scales of the orbital parameters ( 2 ) were calculated . 
   the notation is as follows : @xmath125 is the time in years , @xmath126 is the radius measured in solar radii , @xmath127 is the luminosity measured in solar luminosities , @xmath128 is effective surface temperature measured in kelvins , h1 is the hydrogen mass fraction at the stellar center , and @xmath129 is the mean stellar density measured in @xmath130 .
summary: we investigate the change in the orbital period of a binary system due to dynamical tides by taking into account the evolution of a main - sequence star . 
   three stars with masses of one , one and a half , and two solar masses are considered . 
   a star of one solar mass at lifetimes @xmath0 yr closely corresponds to our sun . 
   we show that a planet of one jupiter mass revolving around a star of one solar mass will fall onto the star in the main - sequence lifetime of the star due to dynamical tides if the initial orbital period of the planet is less than @xmath1 days . 
   planets of one jupiter mass with an orbital period@xmath2 days or shorter will fall onto a star of one and a half and two solar masses in the mainsequence lifetime of the star .