Kinetic Energy and the Free Energy Principle in the Birth of Human Life

: The retrospective noninterventional study investigated the kinetic energy of video images of 18 fertilized eggs (7 were normal and 11 were abnormal) recorded by a time-lapse device leading up to the beginning of the first cleavage. The norm values of cytoplasmic particles were measured by the optical flow method. Three phase profiles for normal cases were found regarding the kinetic energy: 2.199 × 10 − 24 ± 2.076 × 10 − 24 , 2.369 × 10 − 24 ± 1.255 × 10 − 24 , and 1.078 × 10 − 24 ± 4.720 × 10 − 25 (J) for phases 1, 2, and 3, respectively. In phase 2, the energies were 2.369 × 10 − 24 ± 1.255 × 10 − 24 and 4.694 × 10 − 24 ± 2.996 × 10 − 24 (J) (mean ± SD, p = 0.0372), and the time required was 8.114 ± 2.937 and 6.018 ± 5.685 (H) ( p = 0.0413) for the normal and abnormal cases, respectively. The kinetic energy change was considered a condition for applying the free energy principle, which states that for any self-organized system to be in equilibrium in its environment, it must minimize its informational free energy. The kinetic energy, while interpreting it in terms of the free energy principle suggesting clinical usefulness, would further our understanding of the phenomenon of fertilized egg development with respect to the birth of


Introduction
Biological systems consume energy to carry out various vital functions [1].In the partitioning phenomenon during human cells' growth process, there are naturally potential changes in aspects such as kinetic energy, chemical energy, and thermal energy within the cell.Although measurements are extremely difficult with current measurement technology for cells, kinetic energy can be determined by observing the movement of microscopic intracytoplasmic particles within a fertilized egg using a sophisticated microscope that provides detailed images.It is therefore possible to measure kinetic energy's dynamics over the course of the growth process.Cytoplasmic dynamics measured by particle image velocimetry that was based on identifying matching sub-regions within sequential pairs of images in time [2] were reported regarding human oocytes [3,4] and mouse embryos [2,5,6].In a normal mouse oocyte, actin filament flow generates cytoplasmic streaming [7].However, to our knowledge, there is no principled biological explanation for the phenomenon.
The free energy principle [8] states that, for any self-organized system to be in equilibrium in its environment, it must minimize its informational free energy [9].In recent years, substantial research has been conducted on this principle for describing the brain activity mechanism.As the principle also explains the function of a single neuron [10], the hypothesis that the principle can be applied to fertilized eggs also warrants consideration.The lead author (Y.Miyagi) has indirectly estimated free energy values by analyzing fetal facial expressions with artificial intelligence for fetal brain activity to obtain chaotic dimensions [11,12].However, there are no reports on free energy evaluation in fertilized eggs.
If the hypothesis is correct, the free energy fluctuations, especially decreases, should be observed during the developmental stages.
In this study, using an original microscope for high-precision time-lapse observation immediately after fertilization [13][14][15], we investigated the characteristics of intracytoplasmic particle movement observed after sperm penetration in fertilized eggs that underwent normal first cleavage and compared them with abnormal cases that failed at normal cleavage.We then investigated whether the free energy principle could explain the dynamics of kinetic energy.

Patients
Informed consent was obtained from participants at Mio Fertility Clinic between 31 October 2004 and 12 September 2008, with completely deidentified data enrolled.The Institutional Review Board of Mio Fertility Clinic conducted this retrospective, noninterventional study (Institutional Review Board number: MFC-2004-6, 1 June 2004).We defined cases in which the first cleavage occurred within 56 h and a live birth was obtained later as the "normal" group, and those in which the first cleavage did not occur up to 56 h and embryo transfer could not be performed as the "abnormal" group.

Image Capturing Method
A novel cinematography system for time-lapse observations, developed by Mio and colleagues and described elsewhere [13][14][15] was used in this study.Briefly, an inverted microscope (IX-71; Olympus, Tokyo, Japan) with Nomarski differential interference contrast optics (Olympus) and a micromanipulator (Narishige, Tokyo, Japan), which was covered with a handcrafted chamber made of acrylic resin, were used.An air heater was placed in the corner of the chamber to maintain the optimal temperature.Our system also contained a small acrylic chamber (15 cm × 15 cm × 3 cm) surrounded by a small water bath on the stage of the microscope.A glass Petri dish with a microdrop of culture medium (5 µL) was placed in the center of the small chamber.Humidified CO 2 gas was infused into the chamber through the water bath.The volume of flowing CO 2 and the temperature within the chamber were adjusted to give the optimal values (temperature of 37 • C ± 0.3 • C, and pH of 7.45 ± 0.03).To obtain ideal culture conditions in the microdrop of culture medium that was covered with mineral oil (Sage, Pasadena, CA, USA), the following settings were used: CO 2 flow, 40 mL/min; air heater in the large chamber set to 38.0 • C; and the thermoplate on the microscope stage set to 41.8 • C. The inverted microscope was equipped with a CCD digital camera (Roper Scientific Photometrics, Tucson, AZ, USA), which was connected to the computer and display by MetaMorph (Universal Imaging Co., Downingtown, PA, USA).A glass Petri dish with a microdrop of culture medium was placed in the center of the small chamber.The oocytes were placed on the dish, and then the sperm were attached.Digital images were taken at 2 min intervals until the first cleavage began, or up to 56 h.
The image files were transferred offline to Medical Data Labo (Okayama, Japan).The analysis period was from the time the sperm reached the egg cell to the beginning of the first cleavage.

Vector Data Acquisition
The image was adjusted and cropped so that the fertilized egg was placed in the center of the 262 × 262 pixel square.The motion of particles measuring approximately 2 pixels at 706 fixed, evenly distributed grid points was extracted as vector time series data.
When the still image at time t was image(t), the difference in pixels between image(t) and image(t + ∆t) at each grid point was extracted to obtain the norm value using optical flow, which is the pattern of apparent motion of image objects between two consecutive frames caused by the movement of the object [16].From the norm value, the particle velocity v(t) (m/s) was calculated.The particle diameter was 0.8% of the egg cell diameter as measured on the image.

Kinetic Energy Calculation
By setting the oocyte diameter to 120 µm [17] and mass to 1.62 × 10 −12 kg, assuming the mass of an average human cell, the mass of the particle m was determined as 8.303 × 10 −19 kg.The kinetic energy value E m (t) at a certain time t was then obtained with the following equation: The observation period was divided into appropriate phases as needed in accordance with the change in kinetic energy.The phase period, kinetic energy value, kinetic energy value per hour, and respective regression functions were analyzed to compare the normal and abnormal groups.A regression line was obtained for each phase using the logarithm of the energy value as the dependent variable.

Application of the Free Energy Principle
Let W be the total workload and ∆F be the difference in free energy between the transient energy steady state at an appropriate phase from immediately after fertilization to the start of the first cleavage.As the two states are considered to follow a canonical distribution, the Jarzynski equality yields the following equation [18]: can be estimated.We investigated whether the following free energy principle can explain the important kinetic energy changes during the process of fertilizing eggs up to the first cleavage.
The variational free energy, F, in generating fetal expressions using the free energy principle is as follows: (5) x ∈ a, ∼ s = (s 1 , s 2 , . . .s t ), where a is actions, D KL is Kullback-Leibler divergence [19], m is a model, o t is situational information, P is generative density, Q is recognition density, s t is hidden states originally built into the fertilized egg, u is prediction of the result of causing an action [20], Ω is a set of information about the situation, µ is sufficient statistics [21], µ t is internal observational information of the fertilized egg itself, and x is actions that could be observed.
If the equilibrium state is temporarily reached and particle motion is barely observed, E m = E min when E min is the minimum value, then the variational free energy F will be the minimum value at this time.
That is min where therefore, We then considered the distance the vector travels at the observed fixed point in the image as the amount of action in the free energy principle and examined whether the fertilized egg development process followed this principle.We also developed an example model of variational free energy change with particle motion.

Statistical Analysis
Wolfram Language and Mathematica 12.3 (Wolfram Research, Champaign, IL, USA) was used for all analyses as well as statistical analyses: analysis of variance (ANOVA) test with Scheffé's method, Bartlett's test for variance, Chi-square test, Kruskal-Wallis test, linear regression analysis, Mann-Whitney test, t-test, and variance test.We set p < 0.05 as statistically significant.

Results
There were 7 normal and 11 abnormal cases enrolled.The maternal age of the normal and abnormal cases was a mean ± standard deviation (SD) of 31.86 ± 4.06 (range: 27-38) and 31.45 ± 4.29 (range: 27-39), respectively.Maternal ages were not significantly different (p = 0.845).The fraternal age of the normal and abnormal cases was 32.42 ± 4.54 (range: 27-41) and 36.18 ± 5.90 (range: 25-50), respectively.The fraternal ages were not significantly different (p = 0.061).The sterility periods (months) of the normal and abnormal cases were 29.29 ± 19.89 (range: 12-60) and 23.90 ± 16.05 (range: 10-60), respectively.Sterility periods were not significantly different (p = 0.359).The speculated reasons for sterility were two endometrioses, two male infertility, and three unknowns for the normal group, and three endometrioses, three male infertility, and five unknowns for the abnormal group.The speculated reasons for sterility were not significantly different (p = 0.110).All the assisted reproduction methods, if implemented, were IVF.
The motion of the particles in the fertilized egg was obtained as a vector, and energy changes over time.Figure 1 shows the samples of the actual profile and scheme, with three types we classified by the profile of the kinetic energy change.All normal cases first showed a slow transition (phase 1), followed by a sudden increase in kinetic energy values mid-transition, resulting in the energy peak that we defined as the start of phase 2, then a gradual decrease (phase 2), and, finally, another slow transition (phase 3).This was defined as the "cliff" type.In addition to the 3 cliff-type cases in the abnormal group, 10 cases were biphasic, with early phase 1 and late phase 1 defined as the "hill" type.The remaining case was a hill type that had moved to first cleavage without phase 3.
Figure 2 shows the time series evolution of the relationship among morphological characteristics, observed vectors, and kinetic energy values in the normal case.Brighter tones represent larger vectors.The coloration was uneven, and particle movement in the cytoplasm was not uniform.The time when the morphological features were seen did not coincide with the time of energy change, that is, the phase boundary.Classification of time variation of the kinetic energy of intracytoplasmic particles in a fertilized egg (top panel).Insemination was defined as the starting point, and analysis was conducted up to the start of the first cleavage.All normal cases first showed a slow transition (phase 1; red line), followed by a sudden increase in kinetic energy values mid-transition, then a gradual decrease (phase 2; blue line), and again a slow transition (phase 3; orange line).This was defined as the cliff type.In addition to the 3 cliff-type cases in the abnormal group, 10 cases were biphasic, with early phase 1 and late phase 1, which were defined as the hill type.The remaining case was a hill type Classification of time variation of the kinetic energy of intracytoplasmic particles in a fertilized egg (top panel).Insemination was defined as the starting point, and analysis was conducted up to the start of the first cleavage.All normal cases first showed a slow transition (phase 1; red line), followed by a sudden increase in kinetic energy values mid-transition, then a gradual decrease (phase 2; blue line), and again a slow transition (phase 3; orange line).This was defined as the cliff type.In addition to the 3 cliff-type cases in the abnormal group, 10 cases were biphasic, with early phase 1 and late phase 1, which were defined as the hill type.The remaining case was a hill type that had moved to first cleavage without phase 3. Actual examples of time variation of kinetic energy and regression lines are shown (lower three panels).The horizontal axis shows the elapsed time after insemination (hours).The vertical axis shows the kinetic energy calculated from norm values (J) and expressed in a logarithm.Cliff-type normal cases (upper panel), hill-type abnormal cases with early and late phase 1 (middle panel), and hill-type abnormal cases with early and late phase 1 but without phase 3 (lower panel) are shown.All normal cases were cliff-type.
that had moved to first cleavage without phase 3. Actual examples of time variation of kinetic energy and regression lines are shown (lower three panels).The horizontal axis shows the elapsed time after insemination (hours).The vertical axis shows the kinetic energy calculated from norm values (J) and expressed in a logarithm.Cliff-type normal cases (upper panel), hill-type abnormal cases with early and late phase 1 (middle panel), and hill-type abnormal cases with early and late phase 1 but without phase 3 (lower panel) are shown.All normal cases were cliff-type.
Figure 2 shows the time series evolution of the relationship among morphological characteristics, observed vectors, and kinetic energy values in the normal case.Brighter tones represent larger vectors.The coloration was uneven, and particle movement in the cytoplasm was not uniform.The time when the morphological features were seen did not coincide with the time of energy change, that is, the phase boundary.The coloration was uneven, and particle movement in the cytoplasm was not uniform.The time when the morphological features were seen did not coincide with the time of energy change; that is, the phase boundary.The fertilization cone, male pronucleus, female pronucleus, and halo appearance emerged during phase 2. PN, pronucleus.

Time
The transition time (hours) for each phase in normal cases was 7.429 ± 4.201 and 15.543 ± 5.113 (mean ± SD) for phases 1-2 and phases 2-3, respectively, and the time of the first cleavage start was 25.857 ± 6.203 (Table 1, Figure 3).Duration (hours) was 25.771 ± The coloration was uneven, and particle movement in the cytoplasm was not uniform.The time when the morphological features were seen did not coincide with the time of energy change; that is, the phase boundary.The fertilization cone, male pronucleus, female pronucleus, and halo appearance emerged during phase 2. PN, pronucleus.

Time
The transition time (hours) for each phase in normal cases was 7.429 ± 4.201 and 15.543 ± 5.113 (mean ± SD) for phases 1-2 and phases 2-3, respectively, and the time of the first cleavage start was 25.857 ± 6.203 (Table 1, Figure 3).Duration (hours) was 25.771 ± 6.172, 7.429 ± 4.201, 8.114 ± 2.937, and 10.314 ± 3.062, for the entire period and phases 1, 2, and 3, respectively.There were no significant differences in variance or mean among phases 1, 2, and 3, but the mean and median tended to become longer as the phase progressed.
Table 1.Transition duration for each phase in normal and abnormal cases (J).For normal cases, there were no significant differences in variance or mean among phases 1, 2, and 3, but both the mean and median tended to become longer as the phase progressed.For abnormal cases, there was a difference in variance (p = 0.0235) and a significant difference in mean (p = 0.00367) when analyzing five groups: phase 1, early phase 1, late phase 1, phase 2, and phase 3. Analysis of four groups-phase 1, early and late phase 1, phase 2, and phase 3-showed no difference in variance and a significant difference in mean (p = 0.0263), with phase 2 smaller than the other three groups (p < 0.05).Analysis of all three groups (phases 1, 2, and 3) showed no difference in variance and a significant difference in mean (p = 0.0195), with phase 2 smaller than the other two groups (p < 0.05).The average kinetic energy of the intracytoplasmic particles of each phase in the normal case (J) is shown (middle panel).For normal cases, the variance value was higher in phase 1 (p = 0.0065), but there was no significant difference in mean among the phases.For abnormal cases, there was a difference in variance (p = 0.0227) and a significant difference in mean (p = 0.0101).There was a difference in variance (p = 0.0170) and a significant difference in mean (p = 0.000347) among four groups (phase 1, early and late phase 1, phase 2, and phase 3), with early and late phase 1 larger than the other three groups (p < 0.05).There was a variance (p = 0.000263) and a significant difference in mean (p = 0.0138) for all three groups (phases 1, 2, and 3), with phase 1 larger than the other two groups (p < 0.05) and phase 3 smaller (p < 0.05).cases, duration (hours) was a mean ± standard deviation, median of 25.771 ± 6.172, 25, 7.429 ± 4.201, 6.5, 8.114 ± 2.937, 8.5, and 10.314 ± 3.062, 11.2 for the entire period and phases 1, 2, and 3, respectively.There were no significant differences in variance or mean among phases 1, 2, and 3, but the mean and median tended to become longer as the phase progressed.For abnormal cases, there was a difference in variance (p = 0.0235) and a significant difference in mean (p = 0.00367) when analyzing five groups: phase 1, early phase 1, late phase 1, phase 2, and phase 3. Analysis of four groups-phase 1, early and late phase 1, phase 2, and phase 3-showed no difference in variance and a significant difference in mean (p = 0.0263), with phase 2 smaller than the other three groups (p < 0.05).Analysis of all three groups (phases 1, 2, and 3) showed no difference in variance and a significant difference in mean (p = 0.0195), with phase 2 smaller than the other two groups (p < 0.05).

Total
The slope β 1 of the regression lines for each phase using the logarithm of the energy value as the objective variable in normal and abnormal cases is shown (lower panel).For normal cases, there was no difference in variance among the phases, but there was a difference in mean (p = 0.0015), and the value of phase 2 was significantly smaller (p < 0.05).For abnormal cases, there was a significant difference in variance among phases (p = 2.325 × 10 −10 ) and a difference in mean (p = 1.454 × 10 −6 ), with a significantly smaller mean in phase 2 (p < 0.05).

Kinetic Energy Value
Table 2 and Figure 3 show the average kinetic energy of each phase's intracytoplasmic particles in the normal case.
Table 2.The average kinetic energy of the intracytoplasmic particles of each phase in normal cases (J).In normal cases, the variance value was higher in phase 1 (p = 0.0065), but there was no significant difference in means among the phases.For abnormal cases, when analyzed among the five groups, there was a difference in variance (p = 0.0227) and a significant difference in mean (p = 0.0101).There was a difference in variance (p = 0.0170) and a significant difference in mean (p = 0.000347) among four groups (phase 1, early and late phase 1, phase 2, and phase 3), with early and late phase 1 larger than the other three groups (p < 0.05).There was a variance (p = 0.000263) and a significant difference in mean (p = 0.0138) for all three groups (phases 1, 2, and 3), with phase 1 larger than the other two groups (p < 0.05) and phase 3 smaller (p < 0.05).

Regression Function
A regression line was obtained for each phase using the logarithm of the energy value Y as the objective variable Ŷ with the intercept as β 0 and slope as β 1 (Table 3).
The β 1 (J/H) was −0.123 ± 0.182, −0.249 ± 0.075, and 0.070 ± 0.103 (mean ± standard error (SE)) for phases 1, 2, and 3, respectively (Figure 3).There was no difference in variance among the phases, but there was a difference in means (p = 0.0015), and the value of phase 2 was significantly smaller (p < 0.05).In other words, there was a significant decrease in kinetic energy from hyperkinetic energy in phase 2. The kinetic energy value calculated from the regression function rising from the end of phase 1 to the beginning of phase 2 was 4.071 × 10 −24 ± 3.573 × 10 −24 .
The case-count study results showed that energy in phase 1 decreased in 85.7% (6/7) and increased in 14.3% (1/7) of cases.There was no significant increase in energy in phase 1.There was a significant decrease in energy in phase 2 for 100% (7/7) of cases.Energy in phase 3 increased in 71.4% (5/7) of cases but decreased in 28.6% (2/7) of cases.There was a significant increase in energy in phase 3 for 0 cases, a significant decrease in 14.2% (1/7), and neither a decrease nor an increase in 85.7% (6/7).Table 3. Parameters of regression lines obtained for each phase using the logarithm of the energy value Y as the objective variable Ŷ with the intercept as β 0 and slope as β 1 .In normal cases, there was no difference in variance between the phases, but there was a difference in mean (p = 0.0015), the value of phase 2 was significantly smaller (p < 0.05), and there was a significant decrease in energy in phase 2 for 100% (7/7) of cases.In abnormal cases, there was a significant difference in variance between phases (p = 2.325 × 10 −10 ) and a difference in mean (p = 1.454 × 10 −6 ), with a significantly smaller value in phase 2 (p < 0.05).3 show the time required for the entire period, cliff type phase 1, early phase 1, late phase 1, phase 2, and phase 3 in the abnormal cases.Durations (hours) were 34.733 ± 9.206, 18.000 ± 14.309, 7.412 ± 5.609, 3.713 ± 2.691, 6.018 ± 5.685, and 17.286 ± 7.644 (mean ± SD) for the entire period, cliff-type phase 1, early phase 1, late phase 1, phase 2, and phase 3, respectively.Analysis of all three groups (phases 1, 2, and 3) showed no difference in variance and a significant difference in mean (p = 0.0195), with phase 2 smaller than the other two groups (p < 0.05).

Comparison of Normal and Abnormal Cases
Comparing the overall transit time, the normal and abnormal cases showed 25.857 ± 6.203 and 34.609 ± 9.210 (mean ± SD), respectively.There was no difference in variance between the two groups, with the normal cases reaching first cleavage earlier than the abnormal cases (t-test, p = 0.0425).The total energy values for normal and abnormal cases were 5.647 × 10 −24 ± 2.599 × 10 −24 and 1.601 × 10 −23 ± 9.598 × 10 −24 , respectively, with a large variance (p = 0.000492) for the abnormal cases.The normal cases tended to have lower total energy values than the abnormal cases, though there was no significant difference (p = 0.124).Comparing the total energy per hour, the normal and abnormal cases were 2.110 × 10 −25 ± 8.495 × 10 −26 , and 4.989 × 10 −25 ± 3.017 × 10 −25 , respectively, and the variance significantly differed (p = 0.000742), indicating the total energy per hour was lower in the normal cases than in the abnormal cases (p = 0.0463).

Kinetic Energy and Free Energy Principle
Figure 4 shows an example of the relationship among kinetic energy, sufficient statistics, and free energy calculated every 2 min in phase 2 under certain conditions.As the kinetic energy decreases, the sufficient statistics increase, and the free energy decreases.
When sufficient statics are maximized, kinetic energy and free energy are minimized.In phase 2, the change in kinetic energy was considered a sufficient condition for applying the free energy principle.0.0372). in phase 2 was −0.249 ± 0.075 and −0.770 ± 0.902 (J/H) for the normal and abnormal groups, respectively, with  less in the normal group (p = 0.0372).

Kinetic Energy and Free Energy Principle
Figure 4 shows an example of the relationship among kinetic energy, sufficient statistics, and free energy calculated every 2 min in phase 2 under certain conditions.As the kinetic energy decreases, the sufficient statistics increase, and the free energy decreases.When sufficient statics are maximized, kinetic energy and free energy are minimized.In phase 2, the change in kinetic energy was considered a sufficient condition for applying the free energy principle.An example of the relationship among kinetic energy, sufficient statistics, and free energy calculated every 2 min in phase 2 under certain conditions (lower panel).As the kinetic energy (blue curve) decreases, the sufficient statics (green curve) increase, and the free energy (red curve) decreases (upper panel).When sufficient statics are maximized, kinetic energy and free energy are minimized.In phase 2, the change in kinetic energy was considered a sufficient condition for applying the free energy principle.

Discussion
We have used the free energy principle to explain the biological principle of cytoplasmic dynamics using time-lapse image information.We quantitatively found three phases of kinetic energy change of the particles at the boundary of each elapsed time of 7.429 and 15.543 h in normal cases.Phase 1 is the preparation period for kinetic energy expression, and phase 2 is the period for kinetic energy expression, both of which have almost the same duration.In other words, the time for expression and preparation are almost the same, but the reason for this is unknown.The free energy principle could explain the kinetic energy's behavior in phase 2, where this energy was most active.The fact that Figure 4.An example of the relationship among kinetic energy, sufficient statistics, and free energy calculated every 2 min in phase 2 under certain conditions (lower panel).As the kinetic energy (blue curve) decreases, the sufficient statics (green curve) increase, and the free energy (red curve) decreases (upper panel).When sufficient statics are maximized, kinetic energy and free energy are minimized.In phase 2, the change in kinetic energy was considered a sufficient condition for applying the free energy principle.

Discussion
We have used the free energy principle to explain the biological principle of cytoplasmic dynamics using time-lapse image information.We quantitatively found three phases of kinetic energy change of the particles at the boundary of each elapsed time of 7.429 and 15.543 h in normal cases.Phase 1 is the preparation period for kinetic energy expression, and phase 2 is the period for kinetic energy expression, both of which have almost the same duration.In other words, the time for expression and preparation are almost the same, but the reason for this is unknown.The free energy principle could explain the kinetic energy's behavior in phase 2, where this energy was most active.The fact that sustained exercise was observed is consistent with the theory that continuous energy consumption is required for stabilizing the adapted state [22].
In phase 2 of normal cases, the kinetic energy value was the largest, though not significant (Table 2, Figure 3), and morphological changes such as male and female pronucleus formation, both pronuclei abutted, and halo appearance were observed.Regarding the clinical significance of the magnitude of kinetic energy and time course in relation to outcome, normal cases reached first cleavage earlier than abnormal cases (p = 0.0425).Kinetic energy profiles showed that all normal cases had a sudden increase in kinetic energy in phase 2, while 27.3% (3/11) of abnormal cases were cliff-type, but 72.7% (8/11) of abnormal cases had a slow transition to phase 2. Total kinetic energy tends to be lower in normal cases.Normal cases have less kinetic energy in phase 2. We considered phase 2 to be the most relevant period.Phase 1 seems to be an energetic preparatory state for phase 2, and when entering phase 2, kinetic energy suddenly rises within 2 min, followed by a gradual decrease to end phase 2. The process then enters phase 3, the final preparatory period for the first cleavage.The biochemical energy metabolism, including thermodynamics, inside the fertilized egg is unknown, but there is probably a major change in the metabolic pathway at the time of the phase 2 transition.In addition to kinetic energy, chemical and thermal energy may be present in the cell, but direct measurement of these energies is extremely difficult because sufficient heat dissipation is required for rapid or accurate development [22].
The findings indicate that when an abnormality in the fertilized egg prevents conception and implantation, the free energy principle theoretically suggests the existence of an abnormality in the generative model generative density, inferential model recognition density, internal observational data of the egg itself, or latent variables originally incorporated in the egg.Specifically, biochemical signal transduction abnormalities correspond to generative model abnormalities, abnormal protein synthesis due to DNA abnormalities and protein production disorders due to ribosomal receptor abnormalities correspond to internal observation data abnormalities, and inborn errors in latent variables are considered to correspond to pathological conditions such as functional abnormalities in which the egg cleavage is not set as recognition density.Therefore, fertilized eggs' developmental abnormalities could be explained from the viewpoint of the free energy principle.Recently, many embryo evaluation methods using artificial intelligence [23][24][25][26] and estimation of chromosomal aberrations [27,28] have been reported.However, as this research method is in its early stages, the measurement of various parameters related to time and kinetic energy offers potential for prenatal diagnosis, including embryonic development, prediction of live birth, and detection of chromosomal aberrations and genetic abnormalities.
This study does have limitations.The kinetic energy change was largely divided into three phases, though it is possible to further subdivide the time course and interpret that the energy changes up to the first egg cleavage should have more steps.The appropriate image acquisition interval is unknown.In normal cases, the transition to phase 2 occurs within 2 min and should be observed at least every ≤2 min.Second, the relationship between factors possibly associated with particle movement, such as maternal age, maternal complications, information on sperm, chromosomal abnormalities, and genetic abnormalities, needs examining.
The free energy principle has been widely reported in recent years regarding brain function.Though direct free energy measurement is difficult, we were able to estimate the free energy indirectly.The free energy principle, also known as the survival principle of living organisms, was, through this research, suggested as already existing from the moment of birth.The biology of human development encompasses kinetic energy changes and seems to have minimized the developmental system's informational free energy.

Conclusions
The phenomena up to the initiation of the first cleavage in human fertilized eggs are accompanied by dynamic changes in the kinetic energy of cytoplasmic particles.This suggests they may follow the free energy principle.The measurement of kinetic energy while interpreting it in terms of the free energy principle, suggesting clinical usefulness, would further our understanding of the phenomenon of fertilized egg development with respect to the birth of human life.The free energy principle has probably never been discussed in all organisms, not just humans, at the birth of life.Research on the physiological significance of the free energy principle in humans will begin and will be followed by research on its pathological and clinical significance in the future.

Figure 1 .
Figure 1.Classification of time variation of the kinetic energy of intracytoplasmic particles in a fertilized egg (top panel).Insemination was defined as the starting point, and analysis was conducted up to the start of the first cleavage.All normal cases first showed a slow transition (phase 1; red line), followed by a sudden increase in kinetic energy values mid-transition, then a gradual decrease (phase 2; blue line), and again a slow transition (phase 3; orange line).This was defined as the cliff type.In addition to the 3 cliff-type cases in the abnormal group, 10 cases were biphasic, with early phase 1 and late phase 1, which were defined as the hill type.The remaining case was a hill type

Figure 1 .
Figure 1.Classification of time variation of the kinetic energy of intracytoplasmic particles in a fertilized egg (top panel).Insemination was defined as the starting point, and analysis was conducted up to the start of the first cleavage.All normal cases first showed a slow transition (phase 1; red line), followed by a sudden increase in kinetic energy values mid-transition, then a gradual decrease (phase 2; blue line), and again a slow transition (phase 3; orange line).This was defined as the cliff type.In addition to the 3 cliff-type cases in the abnormal group, 10 cases were biphasic, with early phase 1 and

Figure 2 .
Figure 2. The time series evolution of the relationship among morphological characteristics, observed vectors, and kinetic energy values in a normal case.Brighter tones represent larger vectors.The coloration was uneven, and particle movement in the cytoplasm was not uniform.The time when the morphological features were seen did not coincide with the time of energy change; that is, the phase boundary.The fertilization cone, male pronucleus, female pronucleus, and halo appearance emerged during phase 2. PN, pronucleus.

Figure 2 .
Figure 2.The time series evolution of the relationship among morphological characteristics, observed vectors, and kinetic energy values in a normal case.Brighter tones represent larger vectors.The coloration was uneven, and particle movement in the cytoplasm was not uniform.The time when the morphological features were seen did not coincide with the time of energy change; that is, the phase boundary.The fertilization cone, male pronucleus, female pronucleus, and halo appearance emerged during phase 2. PN, pronucleus.

Figure 3 .Figure 3 .
Figure 3.The transition time for each phase in normal and abnormal cases (top panel).For normal cases, duration (hours) was a mean ± standard deviation, median of 25.771 ± 6.172, 25, 7.429 ± 4.201,Figure 3. The transition time for each phase in normal and abnormal cases (top panel).For normal

Figure 4 .
Figure 4.An example of the relationship among kinetic energy, sufficient statistics, and free energy calculated every 2 min in phase 2 under certain conditions (lower panel).As the kinetic energy (blue curve) decreases, the sufficient statics (green curve) increase, and the free energy (red curve) decreases (upper panel).When sufficient statics are maximized, kinetic energy and free energy are minimized.In phase 2, the change in kinetic energy was considered a sufficient condition for applying the free energy principle.