Smart Seismic Control System for High-Rise Buildings Using Large-Stroke Viscous Dampers Through Connection to Strong-Back Core Frame

An innovative structural control system is proposed for high-rise buildings. A damping layer is provided between a stiff upper core frame suspended from the top of the main building and a stiff lower core frame connected to the building foundation. As the ratio of stiffness of both core frames to that of the main building becomes larger, the relative displacement in the damping layer (damper deformation) approaches to the top floor displacement of the main building. The large displacement of the top floor displacement of the main building is taken full advantage in the proposed control system as most of the total displacement of the main building results from the damper deformation instead of interstory drift. Transformation of the multi-degree-of-freedom (MDOF) model into the single-degree-of-freedom (SDOF) model enables a simplified but rather accurate response evaluation for pulse-type and long-duration earthquake ground motions. The results of the time history response analysis of buildings including this control system are presented for various recorded ground motions. Finally, the effectiveness of the proposed structural control system is discussed from the viewpoint of earthquake input energy.


INTRODUCTION
After emerging in the 1980s and 1990s (Leipholz, 1980;Leipholz and Abdel-Rohman, 1986;Housner et al., 1997), the technique of structural control using active and passive control mechanisms has become main stream in structural engineering for tall and special buildings, e.g., base-isolated buildings (Hanson and Soong, 2001;Christopoulos and Filiatrault, 2006;Takewaki, 2009;Lagaros et al., 2012;Domenico et al., 2019). For civil and building structures, passive control plays a central role. This is because of increased demand for taller buildings and the introduction of new construction materials, accelerating the development of new techniques in the field of passive structural control (for example, Takewaki, 1997;Garcia, 2001;Lavan and Levy, 2005;Aydin et al., 2007;Silvestri and Trombetti, 2007;Whittle et al., 2012). However, the characteristics of near-fault and long-duration ground motions can have a significant effect on tall building response. The effects of near-fault ground motions and long-duration/long-period ground motions on tall and base-isolated buildings have been investigated in detail (Hall et al., 1995;Takewaki et al., 2011). The advantage of passive control systems is that they are robust to unexpected disturbances (Takewaki, 2007(Takewaki, , 2009).
Various passive structural control systems exist for tall buildings (Takewaki, 2009;Lagaros et al., 2012;Takewaki, 2015, 2017;Tani et al., 2017;Domenico et al., 2019). The most popular includes the interstory-type (Takewaki, 2009;Lagaros et al., 2012;Domenico et al., 2019) and soft first story-type (Tani et al., 2017). However, passive structural control systems, able to respond to both near-fault ground motions and long-duration/ long-period ground motions, are very limited (Murase et al., 2013;Takewaki, 2015, 2017;Hayashi et al., 2018). Structural control via passive control devices is difficult for near-fault ground motions because the earthquake energy input is made during very short time interval. For near-fault ground motions, base-isolation system attaining an un-resonant state can be an effective technique. However, base-isolation system requires large site areas to prevent impact to retaining walls and increases construction cost. Moreover, base-isolation systems are vulnerable to resonant long-duration/long-period ground motions.
In this paper, a new structural control system is proposed for high-rise buildings. In this system, a sub core frame is located along a main building and a damping layer is provided between a stiff upper core frame suspended from the top of the main building and a stiff lower core frame (strong-back core frame) attached to the building foundation. Although, the concept is not the same, the strong-back core frame system was investigated in the past (Lai and Mahin, 2015;Palermo et al., 2018). The previous strong-back core frame system aimed at avoiding the plastic deformation concentrations in a few stories by distributing the frame deformations to many stories. The large displacement of the top floor displacement of the main building is taken full advantage in the proposed control system. Transformation of the multi-degree-of-freedom (MDOF) model into the singledegree-of-freedom (SDOF) model enables a simplified but rather accurate response evaluation for pulse-type and long-duration earthquake ground motions which are represented by the double impulse and the multi impulse, respectively. While passive viscous dampers are effective for long-duration ground motions owing to the sufficient time for energy dissipation, they are not necessarily effective for near-fault ground motions. These properties are also investigated in the proposed passive control system. The results of the time history response analysis of buildings including this control system are presented for various recorded ground motions. Finally, the effectiveness of the proposed structural control system is discussed from the viewpoint of earthquake input energy in which the original energy transfer function plays a central role (Housner, 1959(Housner, , 1975Berg and Thomaides, 1960;Housner and Jennings, 1975;Zahrah and Hall, 1984;Uang and Bertero, 1990;Ordaz et al., 2003;Takewaki, 2004a,b).

PROPOSED PASSIVE CONTROL SYSTEM USING LARGE-STROKE VISCOUS DAMPERS THROUGH CONNECTION TO STRONG-BACK CORE FRAME System Overview
Consider a passive control system using large-stroke viscous dampers connected to a strong-back core frame, as shown in Figure 1A. This control system is proposed to reduce the seismic response of high-rise buildings. It is assumed that the main building is used for an apartment house and the lower strongback core frame is used for car parking. The upper strong-back core frame is used only for the stiffness element connecting the top of the main frame to the lower strong-back core frame. To enable the efficient use of passive dampers, a stiff core (strong-back core frame) is constructed and attached to the foundation. Another stiff core frame is hung from the top story stiff sub-assemblage of the main building. The height of the lower strong-back core frame is determined from the architectural user demand based on the number of parked cars. Because this lower strong-back core frame is made of a wall-type reinforced concrete structure with relatively small mass, the horizontal force demand is not significant. Since deformation, or story drift demand, is concentrated in the connecting story, a large stroke is required for viscous dampers. Oil dampers of large stroke typically used usually for base-isolation systems are employed in this system.

Simplification of Controlled Building System Into MDOF Model
The building structure including the proposed control system can be represented by a linear MDOF system as shown in Figure 1B. The main building structure is a reinforced concrete structure of 35 stories. The subsystem consists of two subassemblages, i.e., the lower strong-back core frame and the upper strong-back core frame as shown in Figure 1A. Since a sub frame of large stiffness is required in the upper strong-back, braced systems were used for the upper strong-back. Large-stroke viscous dampers (oil dampers) were installed at the connection layer of the lower and upper strong-back core frames, as shown in Figure 1A. The total damping coefficient of those large-stroke viscous dampers is denoted by C d . The fundamental natural period of the total system by complex eigenvalue analysis was 2.24(s). It is possible to decompose the overall system into three subassemblages with fixed boundaries to characterize the component properties, as shown in Figure 1B. This division was made for defining the stiffnessses of three parts. The masses of the main building and the strong-back core frames were determined from the data of a project building. The fundamental natural periods of these three subassemblages (main building, lower strong-back core frame and upper strong-back core frame) are 2.20, 1.06, 0.51(s). The structural damping is considered only in the main building and the damping ratio is 0.03. Eight oil dampers of 2500[kNs/m]-class are installed and the damping coefficients C d =0 ∼ 5.0 × 10 7 [Ns/m] are considered for parametric analysis. This model is treated as the original model. Because the stiffness of the upper and lower strong-back core frames affects the control performance of the proposed passive damper control system, another model with stiffer lower and upper strongback core frames (10 times larger than the original) is used in later section for comparison (sections Seismic Response of Buildings Including Proposed Passive Control System and Its Simplified Evaluation Method and Earthquake Input Energy and Energy Transfer Function). The original model is described as "the Original Model" and the other model is called "the Stiffer Model."

Eigenvalue Analysis of MDOF Model of Controlled Building System
Since the main building and the strong-back core frames have structural damping proportional to the stiffness and the passive viscous damper has an independent damping coefficient, the overall system has non-proportional damping. Complex eigenvalue analysis was performed for this model. Figure 2A shows the lowest-mode complex mode vector multiplied by the corresponding complex participation factor for three viscous damper damping coefficients for the Original Model. In upper figures, the absolute value, the real value and the imaginary values are plotted. In lower figures, the amplitudes and the phase angles of the lowest-mode complex mode vector at three points (top of the main frame, point just above the damper and point just below the damper) are shown. The distance between the point just above the damper and point just below the damper indicates the damper stroke. Figure 2B illustrates the lowest-mode complex eigenvector multiplied by the corresponding complex participation factor for three viscous damper damping coefficients for the Stiffer Model. It can be observed that the absolute value of the displacement of the subsystem becomes smaller as the damping coefficient becomes larger. In addition, the absolute value of the deformation of the connecting story becomes smaller after it becomes larger for some damper damping coefficient. This means that there exists an optimal damper damping coefficient that effectively reduces the deformation of the main building.

REDUCTION OF MDOF MODEL INTO SDOF MODEL USING CORRECTION FACTOR ON DAMPER DAMPING COEFFICIENT
Since the MDOF model requires complex computation for the idealized inputs (double impulse for near-fault ground motion and multi impulse for long-duration ground motion) and recorded ground motions, a simple SDOF model is introduced to clarify the intrinsic earthquake response properties of the buildings including the proposed control system. The response evaluation method using this SDOF model is called "a simplified evaluation method".
Let M main , k main , and c main denote the total mass of the main building, stiffness of the main building and structural damping coefficient of the main building as an SDOF model. In addition, let M upper , k upper , k lower denote the mass of the upper strong-back core frame, stiffness of the upper strong-back core frame, and stiffness of the lower strong-back core frame. The MDOF model is reduced to an equivalent SDOF model using the correction factor β on the damping coefficient C d of dampers. This correction factor is introduced to account for various factors that could cause differences in damper behavior in the MDOF and SDOF models. A similar procedure was adopted in the reference (Hayashi et al., 2018). Let C d * (= βC d ) denote the corrected damping coefficient. At first, a transformation method from the MDOF model into an equivalent SDOF model with a Maxwell-type spring-dashpot model was introduced as shown in Figure 3 with the help of the correction factor. The stiffness k sub of the sub frame can be derived assuming that k upper and k lower are in series. Then, the Maxwell-type springdashpot model is reduced to a Kelvin-Voigt model as shown in Figure 3. It was confirmed that, since a large deformation occurs at the connection between the upper core frame and the lower core frame and the mass of the main frame is much larger than both core frames, it is reasonable to disregard the mass of the lower core frame and regard the lower core frame as a spring. Table 1 presents the parameters of equivalent SDOF models of three subassemblages shown in Figure 1B for the Original Model.
In the transformation from the Maxwell-type spring-dashpot model to a Kelvin-Voigt model, the following relations on the frequency-dependent damping coefficient c e (ω) and the frequency-dependent stiffness k e (ω) are used.
Substitution of Equation (2) into Equation (1b) leads to the determination of k e . Subsequent substitution of k e into Equation (2) provides an approximate fundamental natural circular frequency ω e for the total system. k e (ω e ) and c e (ω e ) are used subsequently for time-history response analysis. Figure 4 shows the computation flow of the damping coefficient of the viscous damper in the equivalent SDOF model to ensure the validity of using the SDOF model in the response evaluation in relation to the MDOF model. For the appropriate comparison of damping performance with respect to varied damper damping coefficient between the SDOF and MDOF models, a damping ratio normalized by the peak value is introduced and denoted by h with an overbar. The difference between the lowest damping ratios of the SDOF model and the MDOF model is used to evaluate the correction factor β. In this process, the conditionh ′ Figure 5 relates the original damping coefficient C d of the viscous damper installed at the connection of the upper and lower strong-back core frames in the MDOF model and the damping coefficient C d * of the corresponding viscous damper in the equivalent SDOF model for the Original Model and the Stiffer Model. The plot of the 'calculated C d * ' was evaluated by using the equivalence of the normalized lowest damping ratios between the MDOF model and the SDOF model. The increase of damper damping coefficient does not always lead to the reduction of responses because excessively large dampers introduce often induces the addition of complex stiffness. In this case, the damping ratio exhibits a local maximum with respect to damper damping coefficient. The ratio of the damping ratio to such maximum value can be used to assess the performance of the dampers. The ratio of the damping ratio to such maximum value was introduced for the MDOF model and the SDOF model.  Figure 5 also demonstrates that, if the correction factor β is not introduced, i.e., β =1, a large error occurs in the response prediction by means of the SDOF model. In other words, a prediction error in the damper damping   Frontiers in Built Environment | www.frontiersin.org coefficient, 3.80 and 4.58 times, leads to a large error in the response evaluation.

SEISMIC RESPONSE OF BUILDINGS INCLUDING PROPOSED PASSIVE CONTROL SYSTEM AND ITS SIMPLIFIED EVALUATION METHOD
In this section, the seismic response of buildings including the proposed control system under idealized ground motions and recorded ground motions is investigated through the comparison with that of the simplified SDOF model introduced in the previous section. Kojima and Takewaki (2015a) demonstrated that one-cycle sinusoidal wave is a good representation of the main part of pulse-type near-fault ground motions. This double impulse is a simplified version of one-cycle sinusoidal input. When a SDOF is resonant to one-cycle sinusoidal input, the response to the double impulse provides a good estimate of the response to one-cycle sinusoidal input. However, when a MDOF is used, the correspondence between the responses to the double impulse and one-cycle sinusoidal input is not guaranteed, even if both inputs are resonant to the fundamental mode of the MDOF. For this reason, the comparison of the responses of the MDOF under the double impulse and one-cycle sinusoidal input is shown here.

Response of Buildings Including Proposed Control System Subjected to Double Impulse and Its Simplified Evaluation Method
Consider the seismic response of buildings including the proposed passive control system and its simplified evaluation method. As a representative of near-fault ground motions, the double impulse resonant to the building is introduced (Kojima and Takewaki, 2015a). The acceleration input of the double impulse can be expressed bÿ where δ(t) is the Dirac's delta function, V I is the initial input velocity and t 0 is the time interval of two impulses. The corresponding one-cycle sine wave can be described bÿ This one-cycle sine wave is introduced to ensure the validity of the double impulse in connection to actual near-fault ground motions. Some actual near-fault ground motions include an influential pulse-type motion which can be represented by onecycle sine wave (see Kojima and Takewaki, 2016). In Equation (4), ω p = π/t 0 is the input circular frequency and T 1 ′ = 2t 0 is the damped fundamental natural period of the building. V p1 is specified as V p1 = 1.2222V I to guarantee the equivalence of the Fourier amplitudes of the double impulse and the one-cycle sine wave even at the first peak (Kojima et al., 2017).
Let ω e ′ and h e denote the damped natural circular frequency and the damping ratio of the SDOF model. Since the response after the first impulse is a free vibration with zero initial displacement and initial velocity=−V I , the displacement of the SDOF model after the first impulse can be expressed by The displacement of the SDOF model is equal to the sum of the damper deformation u d and the deformation u sub of the lower and upper cores: Since the damper and the lower and upper core frames are in series, the force equivalence can be expressed by From Equation (5), (6a, b), the damper deformation can be obtained from The maximum displacement is obtained at time: t 1 = arctan( 1 − h e 2 /h e )/ω e ′ fromu e (t 1 ) = 0 (see Figure 6A). The time t 2 of the maximum deformation of the damper can be derived by substituting the conditionu d (t 2 ) = 0 into Equation (6): c eue (t 2 ) + k e u e (t 2 ) = 0 (8b) as In summary, the maximum displacement of the mass of the SDOF model and the maximum deformation of the damper can be described by u e (t 1 ) and u d (t 2 ) = u e (t 2 ), respectively. Consider a critical excitation problem for maximizing the response after the second impulse where the velocity input amplitude V I is fixed and the interval t 0 between two impulses is varied.
Since the critical input of the second impulse is taken into account, the displacement of the mass of the SDOF model just before the beginning of the second impulse is 0 and its velocity is V I exp(−πh e 1 − h e 2 ). Therefore, the vibration after the input of the second impulse corresponds to the free vibration from the initial displacement, =0, and the initial velocity, =V I {1 + exp(−πh e 1 − h e 2 )}. As a result, the maximum displacement Since the displacement of the SDOF model mass is different from the top displacement of the MDOF model (see Figure 6B)

Response of Buildings Including Proposed Control System Subjected to Multi Impulse and Its Simplified Evaluation Method
While the transient response was treated for the double impulse, the steady-state response is dealt with for the multi-impulse which can be used to approximate long-period, long-duration ground motion. Consider the multi impulse expressed bÿ This input was introduced by Kojima and Takewaki (2015b) for approximating long-duration sinusoidal input (Takewaki et al., 2011) When we consider the resonant input (the impulse is input at the zero displacement), the vibrations can be described by free vibration of the initial displacement, =0, and the initial velocity: Note, when the input cycle N becomes infinity, the initial velocity at later cycles converges to The maximum displacement of the equivalent SDOF model can be obtained as The maximum deformation of the damper is given by  the input of the double impulse. This characteristic is more significant in the Stiffer Model. An optimal damper damping coefficient exist in both the Original Model and the Stiffer Model. Note, the optimal damper damping coefficient can be obtained by using the SDOF model described in the next section.

Response of Buildings Including Proposed Control System Subjected to Recorded Ground Motions and Its Simplified Evaluation Method
The response characteristics of the building with the proposed control system under recorded ground motions are investigated in this section for various damper damping coefficient. Table 2 shows the employed recorded ground motions. The corresponding velocity response spectra are shown in Figure 9. Figure 10A presents the simplified evaluation of the maximum top displacement and the maximum damper deformation for the Original Model. The maximum top-story displacement and the maximum damper deformation are plotted with respect to the damper damping coefficient. Figure 10B indicates the corresponding evaluation for the Stiffer Model. While a slight difference exists in some ground motions, the proposed equivalent SDOF model provides a good approximate *Introduced by Japanese Governmental Agency for considering the long-duration/longperiod ground motion which is predicted during the next Nankai-Trouph earthquake.
FIGURE 9 | Velocity response spectra of input recorded ground motions.
of the response of the MDOF model. Furthermore, the optimal damper damping coefficient obtained for the multi impulse (Figure 8) minimizes the top displacement of the main building. This phenomenon can be seen clearly under the input of OS1 (long-period, long-duration ground motion). Figure 11A shows the maximum interstory drift angle for various damper damping coefficients for the Original Model and Figure 11B illustrates the corresponding one for the Stiffer Model. It can be seen that a large damper deformation is observed in the Stiffer Model. This means that the increase of stiffness of the strong-back core frame enhances the effectiveness of dampers. Furthermore, the optimal damper damping coefficient obtained for the multi impulse (Figure 8) reduces the maximum interstory drift angle of the main building most effectively. This phenomenon is significant for the input of OS1 (long-period, long-duration ground motion).
For the pulse-type ground motions (JMA Kobe NS, Rinaldi Sta. FN), the interstory drifts in the lower to middle stories are constant independent of the damper damping coefficients. The response reduction rates of the interstory drifts in the upper stories from the response of the model without damper are affected by the damper damping coefficients. In addition, such response reduction rates do not change much for the model with dampers larger than the optimal value. On the other hand, for the long-period, long-duration ground motion (OS1), the interstory drifts in the lower to middle stories are large and the response reduction rates of those are affected by damper quantities. The reduction rate is the largest for the optimal damper damping coefficient.

EARTHQUAKE INPUT ENERGY AND ENERGY TRANSFER FUNCTION
The earthquake input energy is an important measure for evaluating the seismic performance of structures (Housner, 1959;Uang and Bertero, 1990;Takewaki, 2004a,b). The earthquake input energy to an elastic structure can be expressed in the frequency domain (Ordaz et al., 2003;Takewaki, 2004a,b). When the Fourier transform of the ground accelerationü g (t) is denoted byÜ g (ω), the total input energy to the MDOF model can be expressed by where F MDOF (ω) is the energy transfer function defined by In Equation (18), H D (ω), M1 denote the displacement transfer function vector of the MDOF model, the mass matrix of the MDOF model, and the influence vector of 1 s, respectively. The energy dissipated by the dampers can be expressed by (ω) indicates the energy dissipation performance of the dampers. The area of the energy transfer function is meaningful for the random input and the maximum value of the energy transfer function is important for the long-duration sine wave because the Fourier spectrum of the long-duration sine wave has a sharp peak and the maximum value of the energy transfer function correlates with such a peak. Properties for the area and the maximum value of the energy transfer function for increasing damper quantities are disclosed below.
The total input energy to the equivalent SDOF model is: where F SDOF (ω) is the energy transfer function for the SDOF model defined by H SDOF D (ω) is derived by using the Maxwell model for the dashpotspring model (frequency-dependent spring and damping coefficient in the Kelvin-Voigt model). The energy dissipated by the dampers in the SDOF model can be expressed by where F SDOF D (ω) is the energy transfer function for the dampers in the SDOF model: (24) Figure 13A shows the area of the energy transfer function for the overall system and damper with respect to damper damping coefficient in the Original Model and the Stiffer Model. The area of the energy transfer function for the overall system is constant with respect to damper level regardless of the MDOF  model or the SDOF model (Takewaki, 2004a,b). The difference between the MDOF model and the SDOF model is due to the difference of mass (the SDOF model disregards the mass of the lower strong-back core frame). It can also be observed that, while the SDOF model underestimates the area of the energy transfer function of the total input energy in the MDOF model, it overestimates the area of the energy transfer function as an index for energy dissipation performance of the dampers in the MDOF model. There also exists an optimal damper damping coefficient maximizing the area of the energy transfer function for the damper, i.e., the most effective damper damping coefficient. Those damper quantities are shown in Figure 13A. These optimal damper damping coefficients correspond well to the optimal ones, shown in Figure 8, obtained for the minimum top displacement for the multi impulse. Figure 13B presents the maximum value of the energy transfer function with respect to the damper damping coefficient in the Original Model and the Stiffer Model. The SDOF model overestimates the maximum value of the energy transfer function for the total input energy and the damper energy consumption in the MDOF model in the Original Model. As the stiffness of strong-back core frame becomes larger, the energy consumption in the dampers governs the main part of the input energy.

LIMITATIONS OF PROPOSED CONTROL SYSTEM AND SIMPLIFICATION METHOD
To investigate the limitation of the proposed control system, further numerical computation was conducted for various parameters. It was clarified that (i) the stiffness ratio of strongback core frames to the main frame has to possess a large value, (ii) the foundation of the main frame and the lower strong-back core frame has to have sufficient stiffness, (iii) the proposed system is not effective enough for so-called long-period pulse-type motions of large velocity amplitude. As the stiffness ratio of the strong-back core frames to the main frame and the stiffness of the foundation become larger, the effectiveness of the proposed control system is enhanced. The long-period pulse-type motions of large velocity amplitude over 2.0[m/s] as recorded at Kumamoto (Japan) in 2016 are critical and could cause extremely large interstory drifts in the lower stories of the main frame. It was confirmed that the introduction of interstorytype passive control systems into lower stories of the main frame in addition to the proposed control system is effective for such critical ground motions. It was also clarified that the proposed simplification method does not work well in the cases mentioned above.

SUMMARIES AND CONCLUSIONS
A new damper deformation concentration-type structural control system has been proposed for pulse-type and long-duration earthquake ground motions. Summaries and conclusions are as follows: (1) A transformation method from an MDOF model into an SDOF model with a Maxwell-type spring-dashpot model was proposed, including a correction factor for compensating the difference in the damper behavior of both models. Then, the Maxwell-type spring-dashpot model was reduced to a Kelvin-Voigt model for a fixed frequency. The proposed control system and simplification method have some limitations: the stiffness ratio of the strong-back core frames to the main frame has to possess a large value, the foundation of the main frame and the lower strong-back core frame has to have sufficient stiffness, and the proposed system is not effective enough for so-called long-period pulse-type motions of large velocity amplitude over 2.0 [m/s].
(2) A simplified seismic response evaluation method to estimate the maximum deformation of the dampers was proposed for the reduced SDOF model, which was subjected to a resonant double impulse to represent near-fault pulse-type ground motions and a resonant multi impulse to represent long-duration, long-period ground motions. The closedform expression of the response under the critical double impulse and the critical multi impulse provides an accurate prediction of the response for the corresponding one-cycle and long-duration sine waves. By using this method, the maximum top-story displacement of the main building and the maximum deformation of dampers can be estimated. (3) The optimal damper damping coefficient minimizing the response of the main building can be obtained by changing the damper damping coefficient in the simplified response evaluation method. (4) The energy transfer functions characterizing the earthquake energy input were derived for the MDOF model and the SDOF model. The area and the maximum value of the energy transfer functions were employed to evaluate the effectiveness and efficiency of the dampers. The area and the maximum value of the energy transfer functions of the SDOF model can approximate those in the MDOF model. The optimal damper damping coefficients derived from the maximization of the area of the energy transfer function correspond well to the optimal ones obtained for the minimum top displacement for the multi impulse. (5) The response reduction performance from the model without dampers can be enhanced greatly by increasing the stiffness of the lower strong-back core frame attached to the foundation and the upper strong-back core frame hung from the top story stiff sub-assemblage of the main building.

DATA AVAILABILITY STATEMENT
The datasets generated for this study are available on request to the corresponding author.

AUTHOR CONTRIBUTIONS
AK formulated the problem, conducted the computation, and wrote the paper. TM conducted the computation and discussed the results. IT supervised the research and wrote the paper.

FUNDING
Part of the present work is supported by the Grant-in-Aid for Scientific Research (KAKENHI) of Japan Society for the Promotion of Science (Nos. 17K18922, 18H01584) and Takenaka Co.