\section{Introduction}

Mechanical shocks widely exist in many engineering applications, such as aerospace\cite{ECSS2015}, navy\cite{MIL-S-901D1989} and defence engineering\cite{810G}.
These mechanical shocks normally do not introduce damage to the main structures, but they may possibly cause major functional failures of on-board products, e.g., electronic and optical devices and other equipment, which would subsequently result in the partial or total loss of a mission.
Designed products should be tested accordingly in laboratories to ensure their reliability in mechanical shock environment during their service life.
However, a complete reproduction of field mechanical shock environment in a laboratory for product validation is almost impossible.
Unlike other simple impulsive or stationary vibration excitations, it is difficult to define the shock severity of a mechanical shock in time and frequency domains due to its complexity.
A metric is necessary for the measure of shock severity and the comparison between a real shock and the simulant shock reproduced in laboratory environments.
\end{fmtext}

\maketitle

The current metric, shock response spectrum (SRS), was proposed by Biot in 1932\cite{biot1932}.
SRS established a kind of measure for shock severity and shock equivalence in an elastic dynamics way.
It evaluates the modal responses using a standardized dynamical system (i.e. SDOF) and a selected maximum response quantity.
It has been shown that the maximum stress is proportional to the maximum particle velocity for a given material and model shape\cite{Hunt1960,Gaberson1969}.
For a given vibration mode, displacement, velocity and acceleration can be linked by the vibration frequency of the mode, and consequently, the stress can be correlated to the acceleration\cite{Gaberson2012}.
Therefore, the shock severity or its destructive potential is evaluated by the maximum structural response, which can be estimated by a given SRS using a response spectrum analysis method, e.g. absolute sum (ABSSUM) method\cite{Biot1933a, Biot1934a}.
The estimated maximum stress can be directly compared to material strength\cite{ECSS2015} to define the damage level, so that the equivalence among various shocks can be established by their SRS curves.
The initial purpose of this method was to study the effects of earthquakes on buildings and provide a methodology for seismic resistant design.
Later this method was generalized for the analysis of various types of mechanical shocks.
SRS was first adopted as the standard testing specification in 1983\cite{MIL-STD-810D1983} to compare relative shock severity between measured and reproduced shocks.
To date, SRS is still the state-of-the-art in measuring shock severity in testing standards in various engineering sectors, such as aerospace (e.g. NASA-STD-7003A\cite{NASA-STD-7003A2011}), defence (e.g. MIL-STD-810G\cite{810G}) and navy (e.g. MIL-S-901D\cite{MIL-S-901D1989}).
However, the shock equivalence based on SRS is questionable.
If the concerned structure is not specified, the maximum response of the structure cannot be evaluated appropriately using only SRS.
More specifically, the SRS method only provides the supremum for the maximum response range of a given structure.
For example, two shock signals with very close SRS curves can lead to over 40\% difference of the maximum acceleration responses for a multiple degree of freedom (MDOF) system\cite{christian2014}.
Especially for complex devices and equipment with rich contents of middle and high natural frequencies, SRS may not accurately represent the shock severity.

Other attempts have been made to give a better metric of shock severity.
For example, pseudo-velocity shock response spectrum (PVSRS) was suggested to substitute the maximum absolute acceleration\cite{Gaberson2012}.
However, the PVSRS is proposed as a rule of thumb without strict analytical justification and has not been widely accepted.
The fatigue damage spectrum (FDS) intends to combine a fatigue model with the shock response spectrum\cite{VanBaren2015}.
Some 3D contour spectra were proposed to represent the dynamical characteristics for a packaging system by changing the traditional SDOF responding model to 2-DOF\cite{Yu2012a}, 3-DOF\cite{Xu2012a} and tilted support\cite{Chen2014c} models.
These response spectra could be useful in a specific area, but their general applications are limited and they are still based on the same concept of SRS.
Based on the authors' best knowledge, there are no other comparable mechanical methods that may offer a better measure of shock severity than the measure based on SRS.

This study introduces the shock response matrix and shock response contour, from which a new measure of shock severity, i.e. shock severity infimum (SSI), can be extracted using singular value decomposition (SVD) method.
SSI provides the infimum of shock severity for a given shock environment, which together with SRS, offers a dual-spectra bounding the range of shock severity.
The proposed method is validated and illustrated by analysing the response of a simple structure based on the finite element method (FEM).

\section{Shock Response Matrix and Contour}\label{section_src}

When a shock represented by an acceleration-time history excitation is applied to a series of SDOF models, the temporal acceleration response of each SDOF oscillator is calculated.
These responses constitute a shock response matrix, whose two coordinates are time and natural frequency of a series of SDOF oscillators.
It should be noted that the proposed method is also applicable when other response quantities are adopted, e.g., displacement and velocity.

A schematic graph of shock response matrix $\bm{M}$ is presented based on the SRS calculation algorithm.
As shown in Fig.\ref{diagram_src}, time series data ($m$ samples) of temporal responses of all $n$ SDOF oscillators are obtained, which are assembled into an $m\times n$ real matrix $\bm{M}$.
As required in standards, time series sample number is usually much larger than frequency series sample number, and therefore, $m\gg n$ in this case.
\begin{figure}
	\centering
	\includegraphics[width=0.85\linewidth]{Figure/diagram_src}
	\caption{Diagram of shock response matrix $\bm{M}$}
	\label{diagram_src}
\end{figure}
To be in line with the conventional SRS method, the element-wise magnitude of matrix $\bm{M}$ is introduced and represented by a new $m\times n$ matrix $\bm{N}$, i.e.,
\begin{equation}\label{avoid_small}
n_{ij}=|m_{ij}|.
\end{equation}

To show the acceleration responses of all considered SDOFs, the contour plot of the matrix $\bm{N}$ is drawn, which is named shock response contour (SRC) in this paper.
The frequencies and amplitudes are presented in logarithm scale in order to show signals' low-frequency character.
Large negative values may be produced during the use of logarithm expression if there are some very small components in matrix $\bm{N}$.
To focus only on severe responses, the data range shown in the SRC is consistent with its corresponding SRS curve, e.g., from 240 m/s$^2$ to 200,000 m/s$^2$ in this study.

When SRC is projected to the acceleration-frequency plane, the spectrum of the maximum response of the projected SRC is SRS, which is a vector consisting of the maximum of each time series response in matrix $\bm{N}$, i.e.,
\begin{equation}\label{srs_definition}
\bm{v}_{\SRS}=[\max_{1\leq i\leq m} n_{ij}]^\top
\end{equation}
The main problem of SRS is the missing of phase information of peak acceleration responses, which can be visualized by the bright line in SRC.

\begin{figure}
	\centering
	\includegraphics[width=0.75\linewidth]{Figure/rvs}
	\caption{Acceleration-time history of `RVS'}
	\label{rvs}
\end{figure}

\begin{figure}
	\centering
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/rvs_src}
		\caption{SRC of shock `RVS'}
		\label{rvs_src}
	\end{subfigure}
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/rvs_cwt}
		\caption{CWT of shock `RVS'}
		\label{rvs_cwt}
	\end{subfigure}
	\caption{Comparison between SRC and CWT for `RVS' shock signal}
	\label{rvs_src_cwt}
\end{figure}

Fig.\ref{rvs} shows the acceleration-time history of shock `RVS', which is a near-field pyroshock measured in the separation stage of an unnamed re-entry rocket vehicle\cite{irvine2010}.
The SRC of `RVS' illustrated in Fig.\ref{rvs_src} can provide rich information for its temporal structure, and hence, could be potentially used as a time-frequency analysis method with considering the response of each SDOF oscillator to the given shock.
SRC can be compared to the frequently-used method for the analysis of a non-stationary signal in the time-frequency domain, continuous wavelet transform (CWT), which projects a signal on a series of zero-mean basis functions derived from an elementary function by dilations and translations.
The CWT scalogram of `RVS' with analytical Morlet wavelet is plotted in Fig.\ref{rvs_cwt} using Matlab, against the corresponding SRC in Fig.\ref{rvs_src}.
The pattern of colour scale is clearer in SRC since it reveals the features of the mechanical system, e.g. resonant response.
The relationship between SRC and CWT is comparable to the relationship between SRS and fast Fourier transform (FFT) in terms of their mechanical and mathematical representations.
Further studies with focuses on the application of SRC are necessary but are outside the scope of this paper.

\section{Feature Extraction Using Singular Value Decomposition}

The equivalence of shock severity can be measured by a better method based on the features of matrix $\bm{N}$.
SRS is the absolute maximum feature of matrix $\bm{N}$, but it totally ignores the phase information.
Beside absolute maximum, there are other available matrix analysis methods that can be adopted to extract frequency features from matrix $\bm{N}$.
In this paper, the singular value decomposition (SVD) will be applied.

\subsection{Singular Value Decomposition}

The SVD method\cite{Strang2016} can be used to decompose a rank $r$ matrix $\bm{N}$ into $r$ rank-one matrices $\bm{N}_k =[n_{ij}^k]_{m\times n} \ (k=1,2,\cdots r)$, which can be sorted by their singular values in a descending order, i.e.,
\begin{equation}\label{svd_decomposition}
\bm{N}=\bm{U}\bm{\Sigma}\bm{V}^\top=\sum_{k=1}^r \sigma_k \bm{u}_k \bm{v}_k^{\top}=\sum_{k=1}^r \bm{N}_k
\end{equation}
where $\bm{U}$ is an $m\times m$ orthogonal matrix, $\bm{\Sigma}$ is a diagonal $m \times n$ matrix with non-negative real numbers on the diagonal, and $\bm{V}$ is an $n\times n$ orthogonal matrix.
The $\bm{u}_k$ and $\bm{v}_k$ are the $k$th column of $\bm{U}$ and $\bm{V}$, corresponding to the $k$th singular vector in time and frequency domains, respectively, and $\sigma_k$ is the $k$th-order singular value.
The SRC decomposition and the contour plot of each decomposed component ($\bm{N}_k$) are shown in Fig.\ref{diagram_svd}.

\begin{figure}[t!]
	\centering
	\includegraphics[width=0.75\linewidth]{Figure/diagram_SVD2}
	\caption{Diagram of SRC decomposition}
	\label{diagram_svd}
\end{figure}

\subsection{Shock Severity Infimum}\label{Separability}

The matrix $\bm{N}$ can be generally considered as the discrete equivalence of a bivariate function $g(t,f)$ for the shock responses of SDOF oscillators, at regularly spaced discrete values of time $t$ and natural frequency $f$.
If the shock response bivariate function $g(t,f)$ can be factorized into two single-variate functions,
\begin{equation}\label{separable_function}
g(t,f)=u(t) \cdot v(f)
\end{equation}
The separated function in frequency domain $v(f)$ would provide a better description for the metric of shock severity because phase difference will not be an issue here.

However, the bivariate function corresponding to matrix $\bm{N}$ is normally not a separable function.
Therefore, this study intends to find a rank-one (separable) matrix $\widehat{\bm{N}}$ with the maximum fit to $\bm{N}$.
The mathematical expression of this minimization problem is given by
\begin{equation}\label{matrix_approximation_problem}
\text{minimize} \quad \text{over} \  \widehat{\bm{N}} \quad  \fnorm{\bm{N}-\widehat{\bm{N}}}  \quad \text{subject to} \quad \text{rank}(\widehat{\bm{N}})=1
\end{equation}
where $\fnorm{\cdot}$ is the Frobenius norm.
According to the Eckart–Young–Mirsky theorem\cite{Eckart1936}, $\widehat{\bm{N}}=\bm{N}_1$, where $\bm{N}_1$ is the first SVD component of $\bm{N}$.
Thus, $\bm{N}_1$ can be associated with the discretized separable function $g(t,f)$ in Eq.(\ref{separable_function}) in a matrix expression of $[g(t,f)]$.

To be in line with the definition of SRS, a new concept of shock severity infimum (SSI) of matrix $\bm{N}$ is introduced by projecting matrix $\bm{N}_1$ to the acceleration-frequency plane, i.e.,
\begin{equation}\label{n1}
\bm{v}_\SSI=[\max_{1\leq i\leq m} n^1_{ij}]^\top=\sigma_1 \maxnorm{\bm{u}_1} \bm{v}_1
\end{equation}
where $\maxnorm{\cdot}$ is the maximum norm.
The $\bm{u}_\SSI$ is defined correspondingly as $\bm{u}_1$ normalized by its maximum value, i.e.,
\begin{equation}
\bm{u}_\SSI= \dfrac{\bm{u}_1}{\|\bm{u}_1\|_\infty},
\end{equation}
so that $\bm{N}_1$ can be re-expressed as
\begin{equation}\label{separable_function_discrete}
\bm{N}_1=\bm{u}_{\SSI} \cdot \bm{v}_\SSI^\top=[u(t)] \cdot [v(f)]^\top=[g(t,f)].
\end{equation}
It should be noted that the SSI of matrix $\bm{N}$ is the SRS of matrix $\bm{N}_1$.

SSI is not a real description, but an approximation of shock severity metric, as SSI is extracted from $\bm{N}_1$ rather than $\bm{N}$.
Error of approximation comes from the difference between $\bm{N}_1$ and $\bm{N}$.
The relative residual error $\alpha$ of this approximation problem can be evaluated by
\begin{equation}\label{relative_residual_error}
\begin{split}
\alpha &= \bigg( \dfrac{\fnorm{\bm{N}-\bm{N}_1}}{\fnorm{\bm{N}}} \bigg) ^2\\
&= \bigg( \dfrac{ \fnorm{\sum_{k=2}^r \bm{N}_k}}{\fnorm{\bm{N}}} \bigg) ^2\\
&= \dfrac{\sum_{k=2}^{r} \sigma_k^2}{\sum_{k=1}^{r} \sigma_k^2}\\
&= 1-\dfrac{\sigma_1^2}{\sum_{k=1}^{r} \sigma_k^2}
\end{split}
\end{equation}
according to the property of SVD method and Forbenius norm.
$\alpha$ is a measure of the ``distance" of the system $\bm{N}$ from a separability function.
For example, $\alpha$ has been used as the degree of inseparability of a bivariate function in the study of the dynamic ripples in ferret primary auditory cortex\cite{Depireux2001}.
$\alpha$ value close to zero indicates good representation with just one decomposed component, i.e. just SSI can well estimate the shock severity.
In this case, SSI could possibly replace the current SRS technique.
$\alpha$ value approaching unity indicates an increase of inseparability, where only SSI can hardly estimate the shock severity.
For the latter case, both SRS and SSI should also be considered together as a dual spectra technique to bound the real shock severity.

\section{Metric of Shock Severity}

\subsection{Structural response based on modal analysis}

For a linear $n$-degree of freedom system, its $n\times n$ modal matrix of mass normalized eigenvector $\bm{\Phi}=[\phi_{ij}]$ can be calculated, where $\phi_{ij}$ is the modal displacement\cite[ch.11]{Clough2003a}.
Based on the modal analysis method, the actual structural acceleration response at coordinate $i$, i.e., $\bm{a}_i$, is considered as the superposed acceleration-time history response of a series of vibration modes\cite{alexander2009},
\begin{equation}
\bm{a}_i=\bm{M}\bm{x}_i
\end{equation}
where $\bm{x}_i=\bm{p}\circ \bm{\phi}_i$ is the Hadamard product of modal participation factor vector $\bm{p}$ and the eigenvector coefficient $\bm{\phi}_i$ at coordinate $i$.
The vector $\bm{\phi_i}$ is taken from the $i$th row of matrix $\bm{\Phi}$, and usually indicate the concerned location on structure.

The maximum absolute acceleration response has been used to calculate the maximum stress in a structure\cite[p.~477]{ECSS2015}, thus, the maximum acceleration response $\maxnorm{\bm{a}_i}$ can be used to represent the effect of shock severity.
It is evident that the maximum acceleration response calculated from matrix $\bm{N}$ is slightly higher than the maximum acceleration response calculated from matrix $\bm{M}$, according to their respective definitions,
\begin{equation}\label{inequality1}
\maxnorm{\bm{a}_i}=\maxnorm{\bm{M} \bm{x}_i} \lessapprox \maxnorm{\bm{N} \bm{x}}
\end{equation}
where $\bm{x}$ is the entrywise absolute values of vector $\bm{x}_i$.
However, inequality relation based on $\bm{M}$ is difficult for analysis.
Therefore, its corresponding absolute matrix $\bm{N}$ is discussed instead.

The time history response $\bm{N} \bm{x}$ is a linear combination of weighted singular vector $\bm{u}_k$ in time domain, i.e.,
\begin{equation}\label{time_combination}
\begin{split}
\bm{N} \bm{x} & =\bm{N}_1 \bm{x} + \sum_{k=2}^{r} \bm{N}_k \bm{x}\\
& =\sigma_1 \bm{u}_1 \bm{v}_1 ^\top \bm{x} + \sum_{k=2}^{r} \sigma_k \bm{u}_k \bm{v}_k ^\top \bm{x} \\
& = (\sigma_1 \bm{v}_1 ^\top \bm{x}) \bm{u}_1 + \sum_{k=2}^{r} (\sigma_k \bm{v}_k ^\top \bm{x}) \bm{u}_k \\
\end{split}
\end{equation}
$\bm{u}_1$ is the most important singular vector due to its dominant weight and usually the only positive singular vector of $\bm{N}$ in time domain.
By applying Eckart–Young–Mirsky theorem\cite{Strang2016} as shown in Appendix \ref{proof_trend}, the temporal vector $(\sigma_1 \bm{v}_1 ^\top \bm{x}) \bm{u}_1$ approximates the $\bm{N} \bm{x}$ in a least square sense, i.e.,
\begin{equation}\label{trend}
\| \bm{N}\bm{x}-(\sigma_1 \bm{v}_1 ^\top \bm{x}) \bm{u}_1 \|_2 \leq \sigma_2 \| \bm{x} \|_2.
\end{equation}
$(\sigma_1 \bm{v}_1 ^\top \bm{x}) \bm{u}_1$ usually depicts the average trend of the time history response of $\bm{N}\bm{x}$, while the rest of weighted singular vectors, $(\sigma_k \bm{v}_k ^\top \bm{x}) \bm{u}_k$ in time domain, contribute to local oscillations around $(\sigma_1 \bm{v}_1 ^\top \bm{x}) \bm{u}_1$.
The maximum of $\bm{N} \bm{x}$ is always larger than the maximum of $\bm{N}_1 \bm{x}$ according to calculation observations.
Therefore,
\begin{equation}\label{inequality2}
\bm{v}_\SSI^\top\bm{x}=\maxnorm{\bm{N}_1 \bm{x}} \leq \maxnorm{\bm{N} \bm{x}} \leq \bm{v}_\SRS^\top \bm{x}
\end{equation}
according to the proofs in Appendix \ref{proof_inequality2}.
The method to calculate the supremum of maximum response with SRS ($\bm{v}_\SRS$) is termed as absolute response spectrum analysis method\cite{ECSS2015,alexander2009}, which is the most conservative prediction according to Eq.(\ref{inequality2}).

The left inequality relationships in (\ref{inequality2}) is the main assumption based on numerous numerical calculation observations, which will be demonstrated by examples in Section \ref{case}.

\subsection{Margin between SSI and SRS}

The margin between infimum and supremum of a shock can be calculated by the ratio vector $\bm{l}$ in dB, i.e. calculating the amplitude ratio between SRS and SSI in logarithm scale,
\begin{equation}\label{margin_db}
l[j]= 20 \lg (\dfrac{v_\SRS[j]}{v_\SSI[j]}) \ \text{dB}
\end{equation}
where $j$ represents the $j$th discrete frequency.
This frequency-dependent margin can show the uncertainty of SRS-based shock severity over the entire frequency range.
A smaller margin $\bm{l}$ indicates better representativeness of shock severity by either SRS or SSI.
Instead, a larger margin $\bm{l}$ indicates the uncertainty of the SRS to represent shock severity, for which other complementary methods, e.g., SSI, is necessary for the description of the shock severity.

\subsection{Uncertainty in Shock Testing}

\begin{figure}
	\centering
	\includegraphics[width=0.75\linewidth]{Figure/diagram_measure}
	\caption{Measure of shock severity using SRS and SSI methods}
	\label{measure_severity}
\end{figure}

SRS is commonly used to evaluate the relative severity of shocks produced in field and laboratory environments in various industrial standards\cite{ECSS2015,810G,christian2014}.
Compared to the real maximum response, Fig.\ref{measure_severity} shows the measure of shock severity via SSI and SRS methods.
The real maximum response under two given shock excitations can not be calculated without the modal information of a specific structure.
Even if the SRS curves of two shocks are equivalent, it only indicates their severity equivalence in terms of the equivalence between their supremum responses.
The relationships between the magnitudes of their maximum responses to the two shocks are still unknown, which causes the uncertainties for the establishment of the equivalence between two shocks based on SRS.

\begin{figure}
	\centering
	\includegraphics[width=\linewidth]{Figure/fully_conservative}
	\caption{Fully conservative testing in a laboratory environment, (a) field shock, (b) laboratory shock with margin required by dual spectrum; (c) laboratory shock with 3dB/6dB margin. The subscript `F' and `L' represent the field shock and the laboratory shock, respectively.}
	\label{fully_conservative}
\end{figure}

The infimum and supremum of the maximum response to a shock can be represented by its SSI and SRS, respectively.
The range of all possible maximum responses can be perceived by the dual spectrum in an explicit way.
However, even the dual spectrum of both the field and laboratory testing shocks are the same, the response of equipment under field shock may still differ from the laboratory testing response, which may lead to over- or under-testing.
A fully conservative way to conduct laboratory shock testing is to make the SSI$_\text{L}$ of the laboratory testing shock to be equal to or higher than the SRS$_\text{F}$ of the field shock, as shown in Fig.\ref{fully_conservative}(b), i.e., the infimum of maximum response of the laboratory testing shock is no less than the supremum of maximum response of the field shock.
If equipment can pass such a test in a laboratory environment, then it shall also survive under the field shock.

As a comparison, a traditional method suggested by US Army 810G standard\cite{810G} is to apply a margin over the SRS$_\text{F}$ of field shocks, as shown in Fig.\ref{fully_conservative}(c), which considers the stochastic variability in the environment and the uncertainty in any predictive methods employed.
The margin is mainly suggested by engineering experience and statistical rationale, e.g., P95/50 rule\cite{Yunis2007}.
For mechanical shock, a 3dB or 6dB margin is normally added, depending on the degree of test level conservativeness desired\cite[p.~383]{810G}, while for pyroshock this margin is often set to 6dB.
A comparison between $\bm{l}$ margin and ordinary margin in a case study will be discussed in Section \ref{case}.

\section{Case Study}\label{case}

\subsection{Shock Signal Preparation and Analysis}

\begin{figure}[t]
	\centering
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/mis}
		\caption{`MIS'}
		\label{mis}
	\end{subfigure}
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/ras}
		\caption{`RAS'}
		\label{ras}
	\end{subfigure}
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/rss}
		\caption{`RSS'}
		\label{rss}
	\end{subfigure}
	\caption{Time history of three shock signals `MIS', `RAS' and `RSS'}
	\label{time_history_compare}
\end{figure}

Three more shocks are synthesised through different mechanisms but having almost the same SRS of the pyroshock `RVS' described in Section \ref{section_src}.
The mechanical impact shock `MIS' is synthesised with the net zero displacement filter method using a field shock generated by metal-metal impact in laboratory environment\cite{Yan2019}.
The random arrangement shock `RAS' and the reverse sine sweep shock `RSS' are synthesised with Tom Irvine's Matlab script\cite{irvine2018}.
The time histories and SRS curves (Q=10) of these shock signals are shown in Figs.\ref{time_history_compare} and \ref{srs_compare}, respectively.

\begin{figure}[htpb!]
	\centering
	\includegraphics[width=0.75\linewidth]{Figure/srs_compare}
	\caption{Comparison of SRS curves of four shock signals (`RVS', `RAS', `RSS' and `MIS')}
	\label{srs_compare}
\end{figure}

\begin{figure}
	\centering
	\begin{subfigure}{0.48\linewidth}
		\includegraphics[width=\linewidth]{Figure/mis_src}
		\caption{SRC of `MIS'}
		\label{mis_src}
	\end{subfigure}
	\begin{subfigure}{0.48\linewidth}
		\includegraphics[width=\linewidth]{Figure/ras_src}
		\caption{SRC of `RAS'}
		\label{ras_src}
	\end{subfigure}
	\begin{subfigure}{0.48\linewidth}
		\includegraphics[width=\linewidth]{Figure/rss_src}
		\caption{SRC of `RSS'}
		\label{rss_src}
	\end{subfigure}
	\caption{Comparison of SRC plots of various shock signals}
	\label{src_cwt}
\end{figure}

Although these shock signals have almost the same SRS curves, they represent three very different types of shock.
Their time-frequency features are also very different, as shown by their SRC (Q=10) in Figs.\ref{rvs_src_cwt}, and \ref{src_cwt}, respectively.
Therefore, their effects on mechanical systems should be generally different, which, however, can not be indicated by their SRS curves.

\begin{figure}[t]
	\centering
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/rvs_dual}
		\caption{`RVS'}
	\end{subfigure}
	\hfill
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/mis_dual}
		\caption{`MIS'}
	\end{subfigure}
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/ras_dual}
		\caption{`RAS'}
	\end{subfigure}
	\hfill
	\begin{subfigure}[b]{0.49\linewidth}
		\includegraphics[width=\linewidth]{Figure/rss_dual}
		\caption{`RSS'}
	\end{subfigure}
	\caption{Dual spectra and margin $\bm{l}$ of four different shock signals}
	\label{dual_compare}
\end{figure}

The proposed dual spectrum method can improve the ordinary SRS method by providing the infimum of maximum shock responses for any given structures.
Fig.\ref{dual_compare} shows the dual spectrum and margin $\bm{l}$ of these four shocks.
The upper solid lines are ordinary SRS curves, while the lower solid lines are the proposed SSI spectrum extracted from their corresponding SRCs using Eq.(\ref{svd_decomposition}).
The curves of margin $\bm{l}$ are depicted by the dashed line with their units on the right axis.

Although the same supremum is indicated by their respective SRS curves, these four shocks have different margins and SSI curves.
For convenience, the entire frequency domain is divided into low-, mid- and high-frequency ranges by taking the knee frequency as a reference, which is about 2000 Hz in this case.
Generally, the margins in low- and mid-frequency ranges are lower than the margins in high-frequency ranges.
The $\bm{l}$ value can be used to estimate the necessary amplification to conduct fully conservative testing in a laboratory environment.
In low- and mid-frequency ranges, the margins of shocks `MIS' and `RAS' are about 3dB and 6dB, respectively, which can be linked to the 3dB/6dB margin in the 810G standard.
While in the high-frequency range, these margin increase exponentially (i.e. the trend is nearly linear in logarithm coordinate in Fig.\ref{dual_compare}) to 9 dB.
The margin of RSS is much higher than those of other shocks, especially in mid- and high-frequency ranges.
The large margin for `RSS' indicates a large uncertainty if `RSS' signal is used for shock severity test even its SRS is almost the same as the SRS of a field shock signal.

\subsection{FEM Validation}

\begin{figure}
	\centering
	\includegraphics[width=0.75\linewidth]{Figure/cantilever_beam}
	\caption{Illustration of FEA model}
	\label{cantilever_beam}
\end{figure}

In order to compare the shock severity using various methods, the response of a simple structure, i.e., a cantilever beam, is investigated using commercial FEA as an example.
The transient modal dynamic analysis solver in Abaqus 6.14-3 is used in this simple case.
As shown in Fig.\ref{cantilever_beam}, acceleration response is measured at the centre of the free end of the cantilever beam, which is excited by the base motion (shocks) applied at the fixed end in the transverse section direction (i.e. x-direction).
This $150\times10\times10$ mm beam is made of aluminium, whose density, Young's modulus, Poisson's ratio and dilatation wave speed are 2700 kg/m$^3$, 68.9 GPa, 0.32 and 6043 m/s, respectively.
Lanczo's eigensolver is adopted to obtain the natural frequencies up to 25,600 Hz, corresponding to the maximum natural frequency in the SRS shown in Fig.\ref{srs_compare}.
Time increment at $5\times 10^{-6}$ s and mesh size of 2 mm are adopted to meet the simulation requirements\cite{YAN2018}.
To be consistent with the Q=10 damping during SRS calculation, the damping ratio of 0.05 is used in the modal dynamic analysis.
Table \ref{modal_information} shows the modal information of the beam in the transverse direction of the beam.
The accumulated effective mass of all calculated modes is 37.98 g, which is already 93.78\% of the total mass of 40.50 g.

\begin{table}
	\centering
	\begin{threeparttable}
		\caption{Modal information for the cantilever beam}
		\label{modal_information}
		\small{
			\begin{tabular}{lllll}
				\hline
				Mode No. & $f$ (Hz) & $\Gamma$ & $\Phi$ & $M_{\text{eff}}$ (kg) \\
				\hline
				1  & 355 & 0.0965    & 6.084      & 0.0093  \\
				2  & 355 & 0.1244     & 7.841      & 0.0154   \\
				3  & 2183 & -0.0306   & 3.409      & 0.0009 \\
				4  & 2183 & -0.0822   & 9.170      & 0.0067  \\
				5  & 4531 & 4.104-12  & 1.656      & 1.684e-23 \\
				6  & 5932 & -0.0499   & -9.237     & 0.0024  \\
				7  & 5932 & 0.0136    & 2.510      & 0.0002 \\
				8  & 8437 & 3.470-14  & -7.268e-06 & 1.204e-27 \\
				9  & 11168 & 0.0092   & -2.294     & 8.465e-05 \\
				10 & 11168 & -0.0361   & 9.010      & 0.0013  \\
				11 & 13591 & -3.508e-12 & -1.656     & 1.230e-23 \\
				12 & 17621 & -0.0287   & -8.837     & 0.0008 \\
				13 & 17621 & 0.0055   & 1.677      & 2.975e-05 \\
				14 & 22640 & -5.409e-12 & 1.65      & 2.926e-23 \\
				15 & 25020 & -0.0219   & 7.88      & 0.0005 \\
				16 & 25020 & 0.0101    & -3.644     & 0.0001 \\
				17 & 25288 & -4.333e-13 & -7.279e-05 & 1.877e-25\\
				\hline
			\end{tabular}
		}
		\begin{tablenotes}
			\item Note: $\Gamma$ is the modal participation factor; $\Phi$ is the mass normalized mode shape; $M_{\text{eff}}$ is the effective mass.
		\end{tablenotes}
	\end{threeparttable}
\end{table}

\begin{table}[htpb!]
	\centering
	\caption{Maximum acceleration responses (m/s$^2$) of the cantilever beam from various methods}
	\label{response}
	\begin{tabular}{ c c c c c c}
		\hline
		Shock\textbackslash Method & FEM & $\maxnorm{\bm{M} \bm{x}_i}$ & $\maxnorm{\bm{N}\bm{x}}$ & $\bm{v}_\SSI^\top \bm{x}$ & $\bm{v}_\SRS^\top \bm{x}$ \\
		\hline
		RVS & 1.90$\times 10^5$ & 1.82$\times 10^5$ & 2.11$\times 10^5$ & 1.60$\times 10^5$ & 2.94$\times 10^5$ \\
		\hline
		MIS & 1.78$\times 10^5$ & 1.78$\times 10^5$ & 2.00$\times 10^5$ & 1.63$\times 10^5$ & 2.94$\times 10^5$ \\
		\hline
		RAS & 2.24$\times 10^5$ & 2.15$\times 10^5$ & 2.29$\times 10^5$ & 1.51$\times 10^5$ & 3.25$\times 10^5$ \\
		\hline
		RSS	& 1.50$\times 10^5$ & 1.46$\times 10^5$ & 1.75$\times 10^5$ & 1.01$\times 10^5$ & 2.94$\times 10^5$ \\
		\hline
	\end{tabular}
\end{table}

The absolute maximum acceleration response of the cantilever beam at the measuring point under different shocks are calculated from various methods based on the modal information in Table \ref{modal_information}.
The results are recorded in Table \ref{response}, whose magnitude relationship diagrams are illustrated in Fig.\ref{bound}.
The acceleration-time histories in Fig.\ref{fem_mxi} obtained from FEM and $\bm{M}\bm{x}_i$ are consistently close, as expected, because they both give the actual acceleration responses in the structure.
Small differences between their maximum responses calculated from $\maxnorm{\bm{M} \bm{x}_i}$ and FEM mainly come from the inaccuracy of the data extraction from the linear interpolation of matrix $\bm{M}$.
Large response differences are observed among different shocks, as shown in Table \ref{response}, although these shocks have almost the same SRS.
Due to the phase difference, the actual maximum response excited by `RAS' is 49.33\% (3.48 dB) higher than that for `RSS', which demonstrates the uncertainty of SRS method.
The maximum response provided by $\maxnorm{\bm{N}\bm{x}}$ is slightly higher than that given by $\maxnorm{\bm{M}\bm{x}_i}$, as shown in Eq.(\ref{inequality1}).
In the diagram of magnitude relationship, the responses calculated by SSI and SRS methods successfully bound the actual maximum response as the infimum and supremum, which is consistent with Eq.(\ref{inequality2}).
Furthermore, the actual maximum response is closer to the infimum response determined by SSI than to the ordinary supremum response determined by SRS, especially for the pyroshock and mechanical shocks.

\begin{figure}
	\centering
	\includegraphics[width=0.9\linewidth]{Figure/bound_bar}
	\caption{Magnitude relationship diagram of maximum acceleration responses of the cantilever beam subjected to various shocks with the same SRS}
	\label{bound}
\end{figure}

\begin{figure}
	\centering
	\includegraphics[width=0.7\linewidth]{Figure/FEM_Mxi}
	\caption{Acceleration response of the cantilever beam under RVS calculated by $\bm{M}\bm{x}_i$ and FEM}
	\label{fem_mxi}
\end{figure}

The first 6 singular vectors ($\bm{u}_k$) of RVS in time domain normalized by their corresponding maximum norm ($\maxnorm{\bm{u}_k}$) are listed in Fig.\ref{time_subplot}, where the first normalized singular vector is the dominant time history component $\bm{u}_{\SSI}$.
The comparison of shock responses of various shocks calculated by $\bm{N}\bm{x}$ and $\bm{N}_1\bm{x}$ are illustrated in Fig.\ref{nx_n1x} along with their corresponding relative residual error $\alpha$.
Since $\bm{N}$ is the non-negative matrix, the linear combination of $\bm{N}$'s columns, i.e. $\bm{N}\bm{x}$, is also a non-negative vector.
It is evident that $\bm{N}_1 \bm{x}$ approaches $\bm{N} \bm{x}$ in terms of signal energy without retaining higher frequency oscillations, which are included in the rest of higher order matrix components.
The residual error $\alpha$ can indicate signal energy difference and its representativeness of $\bm{N}_1 \bm{x}$.
For shock events RVS, MIS and RAS with small $\alpha$ ($\sim 30$), it is reasonable to consider SRC as separable functions, and the response by $\bm{N}_1 \bm{x}$ can provide average trends; while for event RSS with a large $\alpha$ (=0.46), only $\bm{N}_1 \bm{x}$ is insufficient to characterise the trend of $\bm{N} \bm{x}$.

\begin{figure}[!htbp]
	\centering
	\includegraphics[width=\linewidth]{Figure/time_subplot}
	\caption{Weighted singular vectors of RVS in time domain}
	\label{time_subplot}
\end{figure}

\begin{figure}[htpb!]
	\centering
	\begin{subfigure}{0.45\linewidth}
		\includegraphics[width=\linewidth]{Figure/rvs_nx_n1x}
		\caption{RVS}
	\end{subfigure}
	\hfill
	\begin{subfigure}{0.45\linewidth}
		\includegraphics[width=\linewidth]{Figure/mis_nx_n1x}
		\caption{MIS}
	\end{subfigure}
	\begin{subfigure}{0.45\linewidth}
		\includegraphics[width=\linewidth]{Figure/ras_nx_n1x}
		\caption{RAS}
	\end{subfigure}
	\hfill
	\begin{subfigure}{0.45\linewidth}
		\includegraphics[width=\linewidth]{Figure/rss_nx_n1x}
		\caption{RSS}
	\end{subfigure}
	\caption{Comparison of shock responses calculated by $\bm{N}\bm{x}$ and $\bm{N}_1\bm{x}$}

	\label{nx_n1x}
\end{figure}

\section{Conclusions}

This study provides a further theoretical basis for the measure of shock severity.
The traditional shock response spectrum (SRS), as the supremum of the shock response, is insufficient to describe the severity of a shock signal.
Based on the responses of single degree of freedom (SDOF) oscillators to a given shock signal, shock response matrix (SRM) is introduced to record a complete set of responding data of SDOF oscillators in both time and frequency domains.
It is found that the SRC has rich information and can be used as a time-frequency analytical tool in the analyses of mechanical shock signals.
By applying singular value decomposition (SVD) theorem to the SRM, the first singular vector is extracted as shock severity infimum (SSI) based on the conventional SRS concept and calculation process.
The range of the severity of any given shock environment can be explicitly bounded by the dual spectra defined by SRS and SSI curves, which is illustrated by the FEA simulations of a representative structure.
To support the shock resistance design of devices and equipment, the proposed dual spectra can also be used to design more reliable shock tests in laboratory environment since the reproduced shocks in laboratory environments cannot always represent the field shock signals even they have similar SRS curves.
The dual spectra allow the definition of the SRS margin for a given shock specified in terms of its SRS, within which the simplified shock excitations in laboratory environments can be used to substitute field shock signals to certify the devices and equipment.

\enlargethispage{20pt}

\dataccess{The time histories of four shock signals, the source code of SSI method and the modal information of the cantilever beam are released on the GitHub of the author: \href{https://github.com/galois-yan/Shock_Severity_Infimum_Method}{https://github.com/galois-yan/Shock\_Severity\_Infimum\_Method}.}

\aucontribute{Y.Y. and Q.L. designed research; Y.Y. performed research; Y.Y. contributed new analytic tools; Y.Y. analysed data; Q.L. guided the research; Y.Y. and Q.L. wrote the paper.
Both authors agreement to be accountable for all aspects of the work.}

\competing{The authors declare no competing interests.}

\funding{The authors received no funding for this work.}


\section{Introduction and Preliminaries}\label{Introduction}

Let $\mathcal{H}$ denote the class of analytic functions in the unit disk $\mathbb{D}:=\{z\in\mathbb{C}:\, |z|<1\}$. Here  $\mathcal{H}$ is
a locally convex topological vector space endowed with the topology of uniform convergence over compact subsets of $\mathbb{D}$. Let $\mathcal{A}$ denote the class of functions
$f\in \mathcal{H}$ such that $f(0)=0$ and $f'(0)=1$. Let $\mathcal{S}$
denote the class of functions $f\in\mathcal{A}$ which are univalent ({\it i.e.}, one-to-one) in $\mathbb{D}$.
If $f\in\mathcal{S}$ then $f(z)$ has the following representation
\begin{equation}\label{p4_i001}
f(z)= z+\sum_{n=2}^{\infty}a_n z^n.
\end{equation}
A function $f\in\mathcal{A}$ is called starlike (convex respectively) if $f(\mathbb{D})$ is starlike with respect to the origin (convex respectively). Let $\mathcal{S}^*$ and $\mathcal{C}$ denote the class of starlike and convex functions in $\mathcal{S}$ respectively. It is well-known that a function $f\in\mathcal{A}$ is in $\mathcal{S}^*$ if,  and only if, ${\rm Re\,}\left(zf'(z)/f(z)\right)>0$ for $z\in\mathbb{D}$. Similarly, a function $f\in\mathcal{A}$ is in $\mathcal{C}$ if,  and only if, ${\rm Re\,}\left(1+zf''(z)/f'(z)\right)>0$ for $z\in\mathbb{D}$. From the above it is easy to see that $f\in\mathcal{C}$ if, and only if, $zf'\in\mathcal{S}^*$.
 Given $\alpha\in(-\pi/2,\pi/2)$ and $g\in \mathcal{S}^*$, a function $f\in \mathcal{A}$ is said to be close-to-convex with argument $\alpha$ and with respect to $g$ if
$$ {\rm Re\,} \left(e^{i\alpha}\frac{zf'(z)}{g(z)}\right)>0, \,\,\, z\in \mathbb{D}.$$
Let $\mathcal{K}_{\alpha}(g)$ denote the class of all such functions, and $$\mathcal{K}(g):=\bigcup\limits_{\alpha\in(-\pi/2,\pi/2)}\,\mathcal{K}_{\alpha}(g) \mbox{ and } \mathcal{K}_{\alpha}:= \bigcup\limits_{g\in\mathcal{S}^*}\,\mathcal{K}_{\alpha}(g)$$
be the classes of close-to-convex functions with respect to $g$, and close-to-convex function with argument $\alpha$, respectively. The class
$$\mathcal{K}:=\bigcup\limits_{\alpha\in(-\pi/2,\pi/2)}\,\mathcal{K}_{\alpha}=\bigcup\limits_{g\in \mathcal{S}^{*}}\, \mathcal{K}(g)$$
is the class of all close-to-convex functions. It is well-known that every close-to-convex function is univalent in $\mathbb{D}$. Geometrically, $f\in \mathcal{K}$ means that the complement of the image domain $f(\mathbb{D})$ is the union of non-intersecting half- lines. These standard classes are related by the proper inclusions $\mathcal{C}\subsetneq\mathcal{S^*} \subsetneq\mathcal{K}\subsetneq\mathcal{S}$. \\

For $0<\lambda\le1$, let $\mathcal{U}(\lambda)$ be the class of  functions $f\in\mathcal{A}$ satisfying
$$\left |f'(z)\left(\frac{z}{f(z)}\right)^{2}-1\right|<\lambda  \quad\mbox {  for }  z\in \mathbb{D}.
$$
Since $f'(z)(z/f(z))^{2} \ne 0$ in $\mathbb{D}$,  it follows that every function in the class $\mathcal{U}(\lambda)$ is non-vanishing in $\mathbb{D}\setminus \{0\}$. We set $\mathcal{U}:=\mathcal{U}(1)$.
It is known that functions  in $\mathcal{U}(\lambda)$ are locally univalent and functions in the class $\mathcal{U}$ are univalent (see \cite{Aksentev-1958}).
Furthermore, Aksentev \cite{ABU-Muhana-1992} and Ozaki and Nunokawa \cite{ozaki-1972}  have shown that  functions in $ \mathcal{U}(\lambda)$ are univalent, {\it i.e.},
$\mathcal{U}(\lambda) \subseteq \mathcal{S}$  for  $0<\lambda\le 1$.\\

It is worth to mention here some important definitions and results related to the subsets of $\mathcal{H}$.
\begin{defn}
A set $\mathcal{F}\subseteq \mathcal{H}$ is normal if each sequence $\{f_{n}\}$ in $\mathcal{F}$ has a subsequence $\{f_{n_{k}}\}$ which converges to a function $f \in \mathcal{H}$ uniformly on every compact subset of $\mathbb{D}$.
\end{defn}
\begin{defn}
A set $\mathcal{F}\subseteq\mathcal{H}$ is locally bounded if for each point $a\in\mathbb{D}$ there are constants $M$ and $r>0$ such that for all $f \in\mathcal{H}$,
$$|f(z)|\le M \mbox{ for } |z-a|<r.$$
That is, $\mathcal{F}$ is locally bounded if,  about each point $a\in\mathbb{D}$ there is a disk on which $\mathcal{F}$ is uniformly bounded.
\end{defn}
\begin{lem}\cite[Lemma 2.8, p. 153]{Conway-1973}
A set $\mathcal{F} \subseteq\mathcal{H}$ is locally bounded if, and only if, for each compact set $K\subset\mathbb{D}$ there is a constant $M$ such that $$|f(z)|\le M$$
for all $f\in\mathcal{F}$ and $z\in K$.
\end{lem}

\begin{thm} \cite[Montel's Theorem 2.9, p. 153]{Conway-1973} \label{Abhi-Vasu-thm}
A family $\mathcal{F} \subseteq\mathcal{H}$ is normal if, and only if, is locally bounded
\end{thm}
\begin{cor} \cite[Corollary 2.10, p. 154]{Conway-1973} \label{Abhi-Vasu-cor}
A set $\mathcal{F} \subseteq \mathcal{H}$ is compact if, and only if, it is closed and locally bounded.
\end{cor}

\begin{thm}\cite[Theorem 2.6 (Growth Theorem)]{Duren-book-1983}\label{Abhi-Vasu-P1-Growth}
For each $f\in \mathcal{S}$,
$$\frac{r}{(1+r)^{2}}\le |f(z)|\le \frac{r}{(1-r)^{2}},\,\,\,\,   |z|=r\,<\,1.$$
For each $z\in \mathbb{D}$, $z\ne0$, equality occurs if, and only if, $f$ is a suitable rotation of Koebe function.
\end{thm}
For $f$, $g \in \mathcal{H}$, we say that $f$ is subordinate to $g$, written as $f \prec g$ or $f(z) \prec g(z)$, if there
exists an analytic function $\omega: \mathbb{D} \rightarrow \mathbb{D}$ with $\omega(0)=0$ such that
$f(z)= g(\omega(z))$ for $z\in\mathbb{D}$. Furthermore, if $g$ is univalent in $\mathbb{D}$ then $f \prec g$ if, and only if, $f(0)= g(0)$ and $f(\mathbb{D})\subseteq g(\mathbb{D})$.

 If $\mathscr{G}\subseteq\mathcal{H}$, we use the notation $s(\mathscr{G})=\{f: f\prec g \mbox{ for some } g\in\mathscr{G}\}$.
 If $\mathscr{G}$ is a compact subset of $\mathcal{H}$ then it is
not difficult to show that $s(\mathscr{G})$ is compact subset of $\mathcal{H}$ (for
instance, see \cite[Lemma 5.19]{Hallenbeck-MacGregor-1984}).\\

Suppose $X$ is a linear topological vector space and $V \subseteq X$.  A point $x\in V$ is called an extreme point of $V$ if it has no representation of the form $x=ty+(1-t)z, 0<t<1$
as a proper convex combination of two distinct points $y$, $z\in V$. We denote $EV$ the set of extreme points of $V$. The convex hull of a set $V \subseteq X$ is the smallest convex set containing $V$. The closed convex
hull denoted by ${\overline{co}}V$ is defined as the intersection of all closed convex sets containing $V$.
That is,  the closed convex hull of $V$ is the smallest closed convex set containing $V$, which  is the closure of the convex hull of $V$. The Krein-Milman
Theorem asserts that every compact subset of a locally convex topological space is contained in the closed convex hull of its extreme
points (see, for instance, \cite{Dunford-1958}). For a general reference and for many important results on this topic, we refer to \cite{Hallenbeck-MacGregor-1984}. \\

As a first step for application of the knowledge of extreme point of these
classes Brickman {\it et al.} \cite{Brickman-Macgregor-Wilken-1971} pointed out the following general results.

\begin{customthm}{A}\label{p4-theorem001}
Let $\mathscr{G}$ be a compact subset of $\mathcal{H}$ and $J$ be a complex-valued continuous linear functional on $\mathcal{H}$. Then
$\max\{ {\rm Re\, }J(f): f\in \overline{co}\, \mathscr{G}\} = \max\{ {\rm Re\, }J(f): f\in \mathscr{G}\} = \max\{ {\rm Re\, }J(f): f\in E\overline{co}\, \mathscr{G}\}$.
\end{customthm}

\begin{defn}\label{Abhi-Vasu-P1-Defn}
If $\mathcal{F}$ is a convex subset of $\mathcal{H}$ and $J:\mathcal{H}\to\mathbb{R}$ then $J$ is called convex on $\mathcal{F}$ provided
that $J(tf+(1-t)g)\le tJ(f)+(1-t)J(g)$ whenever $f, g\in \mathcal{F}$ and $0\le t\le1$.
\end{defn}

\begin{customthm}{B}\label{p4-theorem005}
Let $\mathscr{G}$ be a compact subset of $\mathcal{H}$ and $J$ be a real-valued, continuous and convex functional
on $\overline{co}\, \mathscr{G}$. Then $\max\{ J(f): f\in \overline{co}\, \mathscr{G}\} = \max\{ J(f): f\in \mathscr{G}\} = \max\{ J(f): f\in E\overline{co}\, \mathscr{G}\}$.
\end{customthm}

The proof of these two results can be found in \cite[Theorem 4.5, Theorem 4.6]{Hallenbeck-MacGregor-1984}. In order to solve such linear
extremal problems over $\mathscr{G}$, it suffices to solve them over the smaller class $E\overline{co}\, \mathscr{G}$. This reduction
thereby becomes an effective technique for solving various linear extremal problems. Using this technique we solve the Zalcman conjecture for
the class $\mathcal{U}$.\\

In 1960s, L. Zalcman posed a conjecture that if a function $f\in\mathcal{S}$ is given by (\ref{p4_i001}) then
\begin{equation}\label{p4_i005}
|a_n^2-a_{2n-1}|\le (n-1)^2 \quad\mbox{ for } n\ge2,
\end{equation}
the equality holds only for the Koebe function $k(z)=z/(1-z)^2$ or its rotation.
It is important to note that the remarkable Zalcman conjecture implies the celebrated Bieberbach conjecture $|a_n|\le n$ for $f\in\mathcal{S}$
(see \cite{Brown-Tsao-1986}). A well-known consequence of the area theorem  shows that  (\ref{p4_i005}) holds good for
$n=2$ (see \cite{Duren-book-1983}). For $f\in\mathcal{S}$, Krushkal has proved the Zalcman conjecture for $n=3$ (see \cite{Krushkal-1995}) and recently
for $n=4,5,6$ (see \cite{Krushkal-2010}). For a simple and elegant proof of Zalcman conjecture for the case $n=3$, we refer to \cite{Krushkal-2010}.
The Zalcman conjecture for functions in the class $\mathcal{S}$ is still  open for $n>6$. However,  using complex geometry and universal Teichm\"{u}ller spaces Krushkal
has  proved it for all $n\ge 2$ in his unpublished work \cite{Krushkal-Unpublished}.\\

The Zalcman conjecture has been proved affirmatively for certain special subclasses of $\mathcal{S}$, such as starlike functions, typically real functions,
close-to-convex functions \cite{Brown-Tsao-1986,Ma-1988} and an observation also demonstrates that the Zalcman conjecture is asymptotically
true (see \cite{Efraimidis-Vukotic-preprint}). Recently, Abu Muhana {\it et al.} \cite{Muhanna-Li-Ponnusamy-2015} solved Zalcman conjecture for the class
$\mathcal{F}$ consists of  the family of  functions $f\in\mathcal{A}$ satisfying the condition
${\rm Re\, }(1+zf''(z)/f'(z))> -1/2$ for $z\in\mathbb{D}$. Functions in the class $\mathcal{F}$  are known to be convex in some direction (and hence close-to-convex and univalent)
in $\mathbb{D}$.  In 1986, Brown and Tsao \cite{Brown-Tsao-1986} proved the Zalcman conjecture for the starlike functions and typically real functions. In 1988, Ma \cite{Ma-1988}
proved that the Zalcman conjecture  for close-to-convex functions. For basic properties of starlike functions, typically real functions and close-to-convex functions
we refer to \cite{Duren-book-1983, Vasu-book}.\\

In 1999, Ma \cite{Ma-1999} proposed a generalized Zalcman conjecture for $f\in \mathcal{S}$ that for $n\ge 2, m\ge 2$, $$|a_{n}a_{m}-a_{n+m-1}|\le (n-1)(m-1),$$
which is still an open problem. Ma  \cite{Ma-1999} has proved this generalized Zalcman conjecture for classes $\mathcal{S}^{*}$ and $\mathcal{S}_{\mathbb{R}}$.  Here $\mathcal{S}_{\mathbb{R}}$  denote
the class of
all functions in $\mathcal{S}$ with real coefficients. In 2017, Ravichandran and Verma \cite{Ravichandran - 2017} proved it for the classes
of starlike and convex functions of given order and for the class of functions with bounded turning.\\

In the present paper, we prove the  Zalcman conjecture and generalized Zalcman conjecture   for the class $\mathcal{U}$ using extreme point theory.
We also  prove  the  Zalcman conjecture and generalized Zalcman conjecture  for the class $\mathcal{F}$ for the initial coefficients.
The organization of the paper is as follows. In Section \ref{Section-2} we prove that the class $\mathcal{U}(\lambda)$ for $0<\lambda\leq 1$ is compact.
In particular the class $\mathcal{U}$ is compact. In Section \ref{Section-3},  we will characterize the  closed convex hull of the class $\mathcal{U}$ and its extreme points.
Then by using extreme point theory,  we prove the  Zalcman conjecture in Section \ref{Section-3} and
generalized Zalcman conjecture    for the class $\mathcal{U}$ in Section \ref{Section-4}. We  prove  the  Zalcman conjecture and generalized Zalcman conjecture
 for the class $\mathcal{F}$ for the initial coefficients in Section \ref{Section-5}.\\

Before we prove our main results we recall some important results which will play vital role in  our proofs. In 2016, Obradovi\'{c} {\it et al.} \cite{Obradvoic-2016} prove the following interested result.

\begin{prop} \cite{Obradvoic-2016}  \label{Abhi-Vasu-P1-Prop-001}
 If $f\in \mathcal{U}(\lambda)$ for $0<\lambda\le 1$, then  for $z\in \mathbb{D}$,
\begin{equation*}\label{Abhi-Vasu-P1-Equ-001abc}
 \frac{f(z)}{z}\prec \frac{1}{(1-z)(1-\lambda z)}.
   \end{equation*}
\end{prop}

Let
$$\mathcal{R}:=\left\{F\in \mathcal{H}:\overline{co}s(F)=\left\{\int_{|x|=1}F(xz)\,d\mu(x):\mu\in \wedge\right\}\right\},
$$
where $\wedge$ denote the set of probability measure on $\partial\mathbb{D}$.  We recall the following well-known result of Hallenbeck {\it et. al}  \cite{Hallenbeck-1989}.

\begin{lem}\cite{Hallenbeck-1989}\label{Abhi-Vasu-P1-Lemma-Hallenbeck-1989}
$\displaystyle  \frac{1}{(1-z)^{\alpha+i\beta}} \in \mathcal{R}$ if, and only if,  $\alpha\ge 1$ and $\beta=0$.
\end{lem}

Let $\mathcal{P}$ denote the class of all analytic functions $p$ in $\mathbb{D}$ with $p(0)=1$ satisfying ${\rm Re\,}p(z)>0$ in $\mathbb{D}$.
Functions in the class $\mathcal{P}$ are called the $Carath\acute{e}odory$   functions and  can be expressed as
\begin{equation}\label{Abhi-Vasu-P1-E}
p(z)=1+\sum\limits_{n=1}^{\infty}c_{n}z^{n}.
\end{equation}
\begin{lem}\cite{Hallenbeck-MacGregor-1984}\label{Abhi-Vasu-P1_Lemma-Hallenbeck-1984}
 $p\in \mathcal{P}$ if, and only if, there is a probability measure $\mu$ on $\partial\mathbb{D}$ such that  $$p(z)=\int_{|x|=1}\frac{1+xz}{1-xz}\, d\mu(x).$$
\end{lem}
Equivalently, in view of the Lemma $\ref{Abhi-Vasu-P1_Lemma-Hallenbeck-1984}$, for $p\in \mathcal{P}$ given by $(\ref{Abhi-Vasu-P1-E})$ can be written as
\begin{equation}\label{Abhi-Vasu-eq}
p(z)=1+\sum\limits_{n=1}^{\infty}c_{n}z^{n}=\int\limits_{0}^{2\pi}\frac{1+e^{it}z}{1-e^{it}z}\,\, d\nu(t).
\end{equation}
On comparing both the sides of $(\ref{Abhi-Vasu-eq})$ we obtain
\begin{equation}\label{Abhi-Vasu-P1-Equ}
c_{n}=2\int\limits_{0}^{2\pi}e^{int}\, d\nu(t).
\end{equation}

\begin{lem}\cite[Lemma 2.3, p. 507]{Ravichandran - 2015}\label{Lemma-Ravichandran-2015}
If $p(z)=1+\sum\limits_{k=1}^{\infty}c_{k}z^{k} \in \mathcal{P}$, then for all $n, m\in \mathbb{N}$,
\[|\lambda c_{n}c_{m}-c_{n+m}|\le
		         \begin{cases}
		             2,  & 0\le\lambda\le1\\[2mm]
		             2|2\lambda-1|,  & \mbox{ elsewhere }.\\
		         \end{cases}
		              \]
If $0<\lambda<1$, the inequality is	sharp for the function $p(z)=(1+z^{n+m})/(1-z^{n+m})$. In other cases, the inequality is sharp for the function $p(z)=(1+z)/(1-z)$.
\end{lem}

\section{Compactness of the set $\mathcal{U}(\lambda)$ }\label{Section-2}
\begin{thm}\label{Compctness}
For $0<\lambda\le1$, the class $\mathcal{U}(\lambda)$ is compact.
\end{thm}
\begin{pf}
In the view of Corollary $\ref{Abhi-Vasu-cor}$, it is enough to show
that $\mathcal{U}(\lambda)$  is closed and locally bounded. Let $\{f_{n}\}$ be a sequence in $\mathcal{U}(\lambda)$ which converges to
$f$ uniformly on every compact subset of $\mathbb{D}$.  Clearly $f(0)=0$, $f'(0)=1$ and $f_{n}' \to f'$ uniformly on every compact subset of $\mathbb{D}$.
Let
$$g_{n}(z)=\left(\frac{z}{f_{n}(z)}\right)^{2}f_{n}'(z)-1 \quad\mbox { and }\,  g(z)=\left(\frac{z}{f(z)}\right)^{2}f'(z)-1.
$$We aim to show that $g_{n} \to g$ uniformly on every compact subset of $\mathbb{D}$. Let $$h_{n}(z)=\frac{f_{n}(z)}{z} \quad\mbox { and }\, h(z)=\frac{f(z)}{z}.$$ Then  $h_{n}(z)\ne0$ and $h(z)\ne0$, for $ z\in \mathbb{D}$ and $ n\in \mathbb{N}$.
Now we  prove that $h_{n} \to h$ uniformly on every compact subset of $\mathbb{D}$.
To show this, it is enough to show that $h_{n}\to h$ uniformly on $\overline{D}=\{z\in\mathbb{C}:|z|\le r\}$, $0<r<1$,  where $D=\{z:|z|< r\}$.  Let
$$M_{n}=\sup_{z\in \overline{D}}|h_{n}(z)-h(z)|=\sup_{z\in \overline{D}}\frac{|f_{n}(z)-f(z)|}{|z|}.$$ Since $h_{n}(z)-h(z)$ is analytic function in $D$ and continuous on $\overline{D}$, so by maximum modulus theorem
 $$\max\{|h_{n}(z)-h(z)|:z\in \overline{D} \}=\max\{|h_{n}(z)-h(z)|:z \in \partial D\}.$$
 Since $f_{n}$ converges to $f$ uniformly on $\overline{D}$ therefore, $M_{n}$ converges to $0$ as $n$ tends to $\infty$.
 Hence $h_{n}$ converges to $h$ uniformly on $\overline{D}$. Therefore, $h_{n}$ converges to $h$ uniformly on every compact subset of $\mathbb{D}$. Since $h_{n} \to h$ uniformly on every compact subset of $\mathbb{D}$ and $h_{n}(z)\ne0$, $h(z)\ne0$, for $ z\in \mathbb{D}$ and $ n\in \mathbb{N}$, it is not difficult to show that $1/h_{n}$ converges to $1/h$ uniformly on every compact subset of $\mathbb{D}$.\\

 If a sequence of continuous functions $f_{n}$ converge uniformly to a continuous function $f$ on some compact set, then $f_{n}$ is uniformly bounded on that compact set. In view of the above discussion,
  it is easy to see that
  $(1/h_{n})^2$ converges to $(1/h)^2$ uniformly on every compact subset of $\mathbb{D}.$
Also, $f_{n}'$ converges to  $f'$ uniformly on every compact subset of $\mathbb{D}$. Hence
  $(1/h_{n})^{2}f_{n}' $ converges to $(1/h)^{2}f' $ unifromaly on every compact subset of $ \mathbb{D}$.
 Therefore $g_{n}$ converges to $g$ uniformly on every compact subset of $\mathbb{D}$.\\

Since  $|g_n(z)|<\lambda$ for each $n\in\mathbb{N}$, we prove that   $|g(z)|\le\lambda$.
Suppose not, then there exists $z_{0} \in \mathbb{D}$ such that $|g(z_{0})|>\lambda$.
Let $\epsilon=|g(z_{0})|-\lambda$.  Then there exists $N\in \mathbb{N}$ such that
$$ |g_{n}(z)-g(z)|<\epsilon=|g(z_{0})|-\lambda
$$
for $ n\ge N$ and  $z\in \mathbb{D}$.  In particular,
$$|g_{N}(z_{0})-g(z_{0})|<|g(z_{0})|-\lambda.
$$
Therefore,
$$|g(z_{0})|-|g_{n}(z_{0})|\le   |g_{N}(z_{0})-g(z_{0})|<|g(z_{0})|-\lambda,
$$
shows that  $|g_{n}(z_{0})|>\lambda$ which is a contradiction and hence  $|g(z)|\le\lambda$ in $\mathbb{D}$.
If there exists some point $z_{0}\in \mathbb{D}$ such that $|g(z_{0})|=\lambda$ then by Maximum Modulus theorem,
 $g$ must be a constant function, which is a contradiction. Therefore, $|g(z)|<\lambda$ for $z\in \mathbb{D}$ and hence $f\in \mathcal{U}(\lambda)$. This shows that  $\mathcal{U}(\lambda)$ is closed.
In view of Theorem $\ref{Abhi-Vasu-P1-Growth}$, it is easy to observe that the class $\mathcal{S}$  is locally bounded. Since the class $\mathcal{U}(\lambda)\subseteq \mathcal{S}$ for $0<\lambda\leq 1$,  it follows that the class   $\mathcal{U}(\lambda)$ is also locally bounded.
Thus, $\mathcal{U}(\lambda)$ is compact.
\end{pf}

In particular,  for $\lambda=1$,  the class $\mathcal{U}$ is compact.

\section{Zalcman Conjecture for the class $\mathcal{U}$}\label{Section-3}

\begin{thm}\label{Close-to-convex-hull--Class-U}
$\overline{co}\, \mathcal{U}$ consists of all functions represented by
\begin{equation*}
f(z)=\int_{|x|=1}\frac{z}{(1-xz)^{2}}\, d\mu(x),
\end{equation*}
where $\mu \in \wedge$. Here $\wedge$ denotes the set of probability measure on $\partial{\mathbb{D}}$.
Further,  $E\overline{co}\, \mathcal{U}$ consists functions of the form
\begin{equation*}
f(z)=\frac{z}{(1-xz)^{2}},\quad |x|=1.
\end{equation*}
\end{thm}
\begin{pf}
Let $f\in \mathcal{U}$ then in view of Proposition $\ref{Abhi-Vasu-P1-Prop-001}$,  we have
$$\frac{f(z)}{z}\prec\frac{1}{(1-z)^{2}}.
$$
Let $F(z)=1/(1-z)^2$ then from Lemma $\ref{Abhi-Vasu-P1-Lemma-Hallenbeck-1989}$, we obtain
$$
\overline{co}s(F)=\left\{\int_{|x|=1}F(xz)\, d\mu(x):\mu\in \wedge\right\} \quad\mbox{ and } \quad E\overline{co}s(F)=\left\{F(xz):|x|=1 \right\}.
$$
Let
$$\mathcal{G}:=\left\{\int_{|x|=1}\frac{z}{(1-xz)^{2}}\, d\mu(x):\mu\in\wedge \right\}
$$
then our aim is to  prove that  $\overline{co}\, \mathcal{U}=\mathcal{G}$.
To prove this, we first prove that $\mathcal{G}$ is convex and compact. The fact that $\mathcal{G}$ is convex follows from the convexity of the set of probability measure on
$\partial\mathbb{D}$ {\it i.e.},  convexity of $\wedge$.  To prove that $\mathcal{G}$ is compact, in the view of Corrolary $\ref{Abhi-Vasu-cor}$, we will show that $\mathcal{G}$ is  closed and locally bounded.
 The fact that $\mathcal{G}$ is closed follows from the weak-star compactness of the set of probability measure on $\partial\mathbb{D}$.
To see that $\mathcal{G}$ is locally bounded,  let $|z|=r<1$ and $f\in \mathcal{G}$, then
 $$f(z)=\int_{|x|=1}\frac{z}{(1-xz)^{2}}\, d\mu(x)
 .$$
 Since $1-|xz|\le |1-xz|$, we have $(1-r)^{2}\le|1-xz|^{2}$ for $|x|=1$ and hence
 \begin{equation}\label{Abhi-Vasu}
 \frac{|z|}{|1-xz|^{2}}\le \frac{r}{(1-r)^{2}}.
 \end{equation}
In view of $(\ref{Abhi-Vasu})$,  we obtain
$$|f(z)|=\left|\int_{|x|=1}\frac{z}{(1-xz)^{2}}\, d\mu(x)\right|\le \frac{r}{(1-r)^{2}}.
$$
Therefore $\mathcal{G}$ is locally bounded.\\

If $f\in \mathcal{U}$ then $f(z)/z\in s(F)$ which  implies that $f(z)/z\in \overline{co}s(F)$. That is,
$$\frac{f(z)}{z}=\int_{|x|=1}\frac{1}{(1-xz)^{2}}\,\, d\mu(x)
$$
for some $\mu \in \wedge$ and hence,
$$f(z)=\int_{|x|=1}\frac{z}{(1-xz)^{2}}\,\, d\mu(x)
$$
for some $\mu \in \wedge$. Therefore  $f \in \mathcal{G}$ and hence   $\mathcal{U}\subseteq \mathcal{G}$.
 Since $\mathcal{G}$ is closed and convex and $\overline{co}\, \mathcal{U}$ is the smallest closed convex set containing $\mathcal{U}$, it follows that
 $\overline{co}\, \mathcal{U}\subseteq\mathcal{G}$.  Since
$$E\mathcal{G}=\left\{\frac{z}{(1-xz)^{2}}:|x|=1 \right\}
$$
and for each $x$ such that $|x|=1$, functions of the form $z/(1-xz)^2$ belong to $\mathcal{U}$, it follows that   $E\mathcal{G}\subseteq\mathcal{U}$.  Since $\mathcal{G}$ is compact and convex, we conclude that
$\overline{co}\, \mathcal{G}=\mathcal{G}$  and hence by Krein-Milman Theorem,  it follows that
$$\mathcal{G}\subseteq\overline{co}\, E\mathcal{G}\subseteq\overline{co}\, \mathcal{U}.
$$
Therefore,  $\mathcal{G}\subseteq\overline{co}\, \mathcal{U}$ and hence $\overline{co}\, \mathcal{U}=\mathcal{G}$ and
$$E\overline{co}\, \mathcal{U}=\left\{\frac{z}{(1-xz)^{2}}:|x|=1\right\}.
$$
This completes the proof.
\end{pf}

Let $f\in \mathcal{S}$ be given by $(\ref{p4_i001})$. Then for fixed $n\in \mathbb{N}$, define the functional $\phi:\mathcal{S}\to \mathbb{C}$ defined by, $\phi(f)=a_{n}^{2}-a_{2n-1}$.
 The rotations of $f\in\mathcal{S}$ be given by
$$g(z)=e^{-i\theta}f(e^{i\theta}z)=z+\sum\limits_{n=2}^{\infty}A_{n}z^{n},
$$
where $A_{n}=a_{n}e^{i(n-1)\theta}$.  Since $\mathcal{S}$ is rotationally invariant, $g\in\mathcal{S}$. A simple computation shows that
$$\phi(g)=A_{n}^{2}-A_{2n-1}=a_{n}^{2}e^{2i(n-1)\theta}-a_{2n-1}e^{2i(n-1)\theta}=e^{2i(n-1)\theta}\phi(f).
$$
This shows that $|\phi(g)|=|\phi(f)|$. That is,  $|a_{n}^{2}-a_{2n-1}|$ is invariant under rotations. Since $\mathcal{U}$ is rotationally invariant, maximizing $|a_{n}^{2}-a_{2n-1}|$ over $\mathcal{U}$ is equivalent to
maximizing ${\rm Re\,}(a_{n}^{2}-a_{2n-1})$ over $\mathcal{U}$. It is easy to show that
\begin{eqnarray}\label{Abhi-Vasu-array}
{\rm Re\,}(a_{n}^{2}-a_{2n-1}) & = & {\rm Re\,}(a_{n}^{2})-{\rm Re\,}(a_{2n-1})\\ \nonumber
&=&({\rm Re\,}(a_{n}))^{2}-({\rm Im\,}(a_{n}))^{2}-{\rm Re\,}(a_{2n-1})\\ \nonumber
& \le & ({\rm Re\,}(a_{n}))^{2}-{\rm Re\,}(a_{2n-1}).
\end{eqnarray}
 In view of $(\ref{Abhi-Vasu-array})$, we maximize $({\rm Re\,}(a_{n}))^{2}-{\rm Re\,}(a_{2n-1})$ over $\mathcal{U}$ to prove the Zalcman conjecture for the class $\mathcal{U}$.

\begin{thm}\label{Zalcman-Conjecture-Class-U}
Let  $f\in \mathcal{U}$  be given by $(\ref{p4_i001})$. Then  $|a_{n}^{2}-a_{2n-1}|\le (n-1)^{2}$ for $n \ge 2$.
This inequality is sharp with equality for the Koebe function and its rotations {\it i.e.}, functions of the form $f(z)=z/(1-xz)^{2}$ where $|x|=1$.
\end{thm}
\begin{pf}
Since $\mathcal{U}\subseteq\mathcal{S}$, for the case $n=2$, the proof of Zalcman Conjecture holds good (see \cite{Duren-book-1983}) as consequence of the area theorem.
For the case  $n=3$, the Zalcman Conjecture has been proved  by Krushkal \cite{Krushkal-2010} for the class  $\mathcal{S}$. Therefore, it suffices to prove the Zalcman Conjecture for $n\geq 4$ for the class $\mathcal{U}$.
For this, for fixed $n\in \mathbb{N}$ we define  the functional $J:\mathcal{S}\to \mathbb{R}$  by
$$J(f)=({\rm Re\,}(a_{n}))^{2}-{\rm Re\,}(a_{2n-1}).
$$
We first prove that  $J$ is  convex on $\overline{co}\, \mathcal{U}$. Let $f,g\in\overline{co}\, \mathcal{U}$ be given by
$f(z)=z+\sum_{n=2}^{\infty}a_{n}z^{n}$ and $g(z)=z+\sum_{n=2}^{\infty}b_{n}z^{n}$. In view of Definition $\ref{Abhi-Vasu-P1-Defn}$,   we show that
$J(tf+(1-t)g)\le tJ(f)+(1-t)J(g)$.  Let
 $tf(z)+(1-t)g(z)=z+\sum_{n=2}^{\infty}A_{n}z^{n}
 $,
 where $A_{n}=ta_{n}+(1-t)b_{n}$. A  computation shows that
\begin{eqnarray*}J(tf+(1-t)g)
&=&({\rm Re\,}(A_{n}))^{2}-{\rm Re\,}(A_{2n-1})\\
&=&({\rm Re\,}(ta_{n}+(1-t)b_{n}))^{2}-{\rm Re\,}(ta_{2n-1}+(1-t)b_{2n-1})\\
&=&({\rm Re\,}(ta_{n}+(1-t)b_{n}))^{2}-t{\rm Re\,}(a_{2n-1})-(1-t){\rm Re\,}(b_{2n-1})\\
&=& t^{2}({\rm Re\,}(a_{n}))^{2}+(1-t)^{2}({\rm Re\,}(b_{n}))^{2}+2t(1-t){\rm Re\,}(a_{n}){\rm Re\,}(b_{n})\\
&&\hspace*{0.25cm}  -t({\rm Re\,}(a_{n}))^{2}-(1-t)({\rm Re\,}(b_{n}))^{2}+t({\rm Re\,}(a_{n}))^{2}\\
&&\hspace*{0.25cm} +(1-t)({\rm Re\,}(b_{n}))^{2}-t{\rm Re\,}(a_{2n-1})-(1-t){\rm Re\,}(b_{2n-1})\\
&=& t^{2}({\rm Re\,}(a_{n}))^{2}-t({\rm Re\,}(a_{n}))^{2}(1-t)^{2}({\rm Re\,}(b_{n}))^{2}-(1-t)({\rm Re\,}(b_{n}))^{2}\\
&&\hspace*{0.25cm}+2t(1-t){\rm Re\,}(a_{n}){\rm Re\,}(b_{n}) t\left(({\rm Re\,}(a_{n}))^{2}-{\rm Re\,}(a_{2n-1})\right)\\
&&\hspace*{0.25cm} +(1-t)\left(({\rm Re\,}(b_{n})^{2}-{\rm Re\,}(b_{2n-1})\right)\\
&=& t(t-1)({\rm Re\,}(a_{n}))^{2}+t(t-1)({\rm Re\,}(b_{n}))^{2}-2t(t-1){\rm Re\,}(a_{n}){\rm Re\,}(b_{n})\\
&&\hspace*{0.25cm} +tJ(f)+(1-t)J(g) \\
&=& t(t-1)({\rm Re\,}(a_{n})-{\rm Re\,}(b_{n}))^{2}+tJ(f)+(1-t)J(g) \\
&=& tJ(f)+(1-t)J(g)-t(1-t)({\rm Re\,}(a_{n})-{\rm Re\,}(b_{n}))^{2}\\
&\le &  tJ(f)+(1-t)J(g).
\end{eqnarray*}
and hence $J$ is a convex functional. In view of Theorem \ref{p4-theorem005}, Theorem $\ref{Compctness}$ and Theorem $\ref{Close-to-convex-hull--Class-U}$,  we consider the function $f_0$ of the form
$$f_{0}(z)=\frac{z}{(1-xz)^{2}}=z+\sum\limits_{n=2}^{\infty}A_{n}z^{n},
$$
where $|x|=1$ and $A_{n}=nx^{n-1}$. Therefore, $\phi(f_{0})=A_{n}^{2}-A_{2n-1}=(n^{2}-2n+1)x^{2n-2}$ and hence $|\phi(f_{0})|=(n-1)^{2}$. For $x=e^{i\theta}$, a simple computation shows that
\begin{eqnarray*}
J(f_{0})&=&({\rm Re\,}(a_{n}))^{2}-{\rm Re\,}(a_{2n-1})\\
&=&n^{2}({\rm Re\,}(x^{n-1}))^{2}-(2n-1){\rm Re\,}(x^{2n-2})\\
&=&n^{2} \cos^{2}(n-1)\theta-(2n-1)\cos(2n-1)\theta\\
&=&n^{2} \cos^{2}(n-1)-(2n-1)(2\cos^{2}(n-1)\theta-1)\\
&=&\cos^{2}(n-1)\theta(n^{2}-4n+2)+2n-1\\
&\le & n^{2}-4n+2+2n-1 \qquad(\mbox{ since } n^{2}-4n+2>0 \mbox{ for } n\ge 4.) \\
&=& n^{2}-2n+1=(n-1)^{2}.
\end{eqnarray*}
 Therefore,  $({\rm Re\,}(a_{n}))^{2}-{\rm Re\,}(a_{2n-1})\le (n-1)^{2}$ and hence
${\rm Re\,}(a_{n}^{2}-a_{2n-1})\le (n-1)^{4}$. This implies that $|a_{n}^{2}-a_{2n-1}|\le(n-1)^{2}$ and the equality holds for the function
of the form $f_{0}(z)=z/(1-xz)^{2}$, where $|x|=1$.  This completes the proof.
\end{pf}

\section{Generalized Zalcman Conjecture for the class $\mathcal{U}$}\label{Section-4}

Let  $p(z)=1+\sum_{k=1}^{\infty}c_{k}z^{k} \in \mathcal{P}$ then from $(\ref{Abhi-Vasu-P1-Equ})$, it is easy to see that for all $n, m\in \mathbb{N}$,
$$\lambda c_{n-1}c_{m-1}-c_{n+m-1}=2\left(2\lambda\int\limits_{0}^{2\pi}e^{i(n-1)t}d\nu(t)\int\limits_{0}^{2\pi}e^{i(m-1)t}d\nu(t)-\int\limits_{0}^{2\pi}e^{i(n+m-2)t}d\nu(t)\right).
$$
By using Lemma \ref{Lemma-Ravichandran-2015},  we obtain
\begin{eqnarray}\label{Abhi-Vasu-P1-Equry}
&& \left|2\lambda\int\limits_{0}^{2\pi}e^{i(n-1)t}d\nu(t)\int\limits_{0}^{2\pi}e^{i(m-1)t}d\nu(t)-\int\limits_{0}^{2\pi}e^{i(n+m-2)t}d\nu(t)\right|\\[5mm]
&&\hspace{5cm} \le
\begin{cases}
 1,  & 0\le\lambda\le1\\[5mm]
 |2\lambda-1|,  & \mbox{ elsewhere }.\nonumber
\end{cases}
\end{eqnarray}

We now prove the generalized Zalcman conjecture for the class $\mathcal{U}$.
\begin{thm}\label{Generalized-Zalcman-Conjecture-Class-U}
Let  $f\in \overline{co}\, \mathcal{U}$  be given by $(\ref{p4_i001})$. Then for $n,m \ge 2$
\[|a_{n}a_{m}-a_{n+m-1}|\le\begin{cases}
   n+m-1,&\mbox{ if }(n,m)\mbox{ is } (2,n), (m,2), (3,3), (3,4), (4,3)\\[5mm]
   (n-1)(m-1),&\mbox{ otherwise. }
   \end{cases}
   \]
   The second  inequality is sharp and the equality holds for the Koebe function and its rotations.
\end{thm}
\begin{pf}
Let $f\in \overline{co}\,\mathcal{U}$ be given by $f(z)=z+\sum_{n=2}^{\infty}a_{n}z^{n}$. Then from Theorem $\ref{Close-to-convex-hull--Class-U}$,
there exists a probability measure $\mu$ on $\partial\mathbb{D}$ such that
$$f(z)=\int_{|x|=1}\frac{z}{(1-xz)^{2}}\, d\mu(x).
$$
Equivalently, there exists a probability measure $\nu$ on $[0,2\pi]$ such that
$$f(z)=\int\limits_{0}^{2\pi}\frac{z}{(1-e^{it}z)^{2}}\, d\nu(t)$$
which can be written as,
\begin{equation}\label{Abhi-Vasu-P1-Equ2}
z+\sum\limits_{n=2}^{\infty}a_{n}z^{n}=z+\sum\limits_{n=2}^{\infty}n\left(\int\limits_{0}^{2\pi}e^{i(n-1)t}d\nu(t)\right)z^{n}.
\end{equation}
By comparing both the sides of $(\ref{Abhi-Vasu-P1-Equ2})$, we obtain
\begin{equation}\label{Abhi-Vasu-P1-Equ3}
a_{n}=n\int\limits_{0}^{2\pi}e^{i(n-1)t}\, d\nu(t).
\end{equation}
Using $(\ref{Abhi-Vasu-P1-Equ3})$, we obtain $$a_{n}a_{m}-a_{n+m-1}=nm\int\limits_{0}^{2\pi}e^{i(n-1)t}\, d\nu(t)\int\limits_{0}^{2\pi}e^{i(m-1)t}\, d\nu(t)-(n+m-1)\int\limits_{0}^{2\pi}e^{i(n+m-2)t}\, d\nu(t)$$
which can be written as
\begin{eqnarray}\label{Abhi-Vasu-P1-Equ4}
& a_{n}a_{m}-a_{n+m-1}=& n+m-1\left[\frac{nm}{n+m-1}\int\limits_{0}^{2\pi}e^{i(n-1)t}d\nu(t)\int\limits_{0}^{2\pi}e^{i(m-1)t}d\nu(t)\right.\\ \nonumber
&&\qquad -\left.\int\limits_{0}^{2\pi}e^{i(n+m-2)t}d\nu(t)\right].
\end{eqnarray}
Comparing $(\ref{Abhi-Vasu-P1-Equry})$ and $(\ref{Abhi-Vasu-P1-Equ4})$,  we obtain
\begin{equation}\label{Abhi-Vasu-P1-Equ5}
\lambda=\frac{nm}{2n+2m-2}\, .
\end{equation}
For $\lambda\le1$, we see that $(\ref{Abhi-Vasu-P1-Equ5})$ implies $nm-2n-2m+2\le0$ which is equivalent to $(n-2)(m-2)\le2$.
If any of $n,m$ is $2$ and for the pairs $(n,m)=(3,3), (3,4), (4,3)$ the inequality $(n-2)(m-2)\le2$ holds. In all other choices of $n,m$ we get $\lambda>1$. Therefore, we have
 \begin{eqnarray}\label{Abhi-Vasu-P1-Equry1}
&& \left|\frac{nm}{n+m-1}\int\limits_{0}^{2\pi}e^{i(n-1)t}d\nu(t)\int\limits_{0}^{2\pi}e^{i(m-1)t}d\nu(t)-\int\limits_{0}^{2\pi}e^{i(n+m-2)t}d\nu(t)\right|\\[5mm]
 &&\hspace{5cm}  \le
 \begin{cases}
   1,  & (2,m), (n,2), (3,3), (3,4), (4,3) \\[5mm]

  \frac{(n-1)(m-1)}{n+m-1},  & \mbox{ other pairs of } (n,m). \nonumber
  \end{cases}
  \end{eqnarray}
   Therefore from $(\ref{Abhi-Vasu-P1-Equ4})$ and $(\ref{Abhi-Vasu-P1-Equry1})$,  we obtain
    \[|a_{n}a_{m}-a_{n+m-1}|\le\begin{cases}
   n+m-1,&\mbox{ if } (n,m) \mbox{ is } (2,n), (m,2), (3,3), (3,4), (4,3)\\[5mm]
   (n-1)(m-1),&\mbox{ otherwise.}
   \end{cases}
   \]
   The second inequality is sharp and equality holds for the Koebe function and its rotations.
   \end{pf}

   \section{Zalcman Conjecture for certain subclass of Close-to-convex functions}\label{Section-5}
      In this section,  we consider the class of close-to-convex functions with argument $0$ and with respect to Koebe function $k(z)=z/(1-z)^{2}$. More precisely,  let
      $$\mathcal{F}=\{f\in \mathcal{A}: {\rm Re\,}(1-z)^{2}f'(z)>0, \quad z\in \mathbb{D}\}.
      $$
      Clearly the functions in $\mathcal{F}$ are convex in the positive direction of the real axis. The region of variability of the class $\mathcal{F}$ has been studied  by Ponnusamy {\it et al} in
      \cite{Samy-Vasu-Yanagihara}. In 2017, Ali and Vasudevarao \cite{Vasu-2017} obtained the sharp logarithmic  coefficient for functions in the class $\mathcal{F}$. \\

 Our aim is to solve the Zalcman conjecture  for $n=2$ and  generalized Zalcman Conjecture  for $n=2, m=3$.
      \begin{thm}\label{Abhi-Vasu-P1-subclass}
      Let $f\in\mathcal{F}$ given by $(\ref{p4_i001})$. Then
      \begin{enumerate}
      \item[(i)] $|a_{2}^{2}-a_{3}|\le 1$
      \item[(ii)] $|a_{2}a_{3}-a_{4}|\le 2$.
      \end{enumerate}
      These inequalities are sharp with equality for the Koebe function $k(z)=z/(1-z)^{2}$ and its rotations.
    \end{thm}
    \begin{pf}
    For $f\in\mathcal{F}$, let $g(z)=(1-z)^{2}f'(z)$. Clearly ${\rm Re\,}g(z)>0$. Since ${\rm Re\,}g(z)>0$ in $\mathbb{D}$,
    there exists an analytic function $\phi:\mathbb{D}\to\mathbb{D}$  such that
     \begin{equation}\label{Abhi-Vasu}
    \phi(z)=\frac{g(z)-1}{g(z)+1}.
    \end{equation}
    Clearly,  $\phi(0)=0$. Let
   $$\phi(z)=\sum\limits_{n=1}^{\infty}c_{n}z^{n}=c_{1}z+c_{2}z^{2}+c_{3}z^{3}+c_{4}z^{4}+\ldots.
   $$
   In view of the Schwarz Lemma, we get $|c_{1}|\le 1$. From $(\ref{p4_i001})$ we have
    $$g(z)=(1-z)^{2}f'(z)=(1+z^{2}-2z)(1+\sum\limits_{n=2}^{\infty}na_{n}z^{n-1}).$$
    From $(\ref{Abhi-Vasu})$, we have
    $$\phi(z)(g(z)+1)=g(z)-1.
    $$
    A simple computation shows that
    \begin{eqnarray}\label{Abhi-Vasu-00001}
    & a_{2}=& 1+c_{1}\\ \nonumber
   & a_{3}=& \frac{1}{3}\left(2c_{1}^{2}+4c_{1}+2c_{2}+3\right)\\ \nonumber
   &a_{4}=&\frac{1}{4}\left(2c_{1}^{3}+4c_{1}^{2}+4c_{1}c_{2}+6c_{1}+4c_{2}+2c_{3}+4\right).
    \end{eqnarray}

    Therefore using (\ref{Abhi-Vasu-00001}) we have
    \begin{eqnarray*}
    a_{2}^{2}-a_{3} &=& (1+c_{1})^{2}-\frac{1}{3}\left(2c_{1}^{2}+4c_{1}+2c_{2}+3\right)\\
    &=& \frac{1}{3}\left(3(1+2c_{1}+c_{1}^{2})-2c_{1}^{2}-4c_{1}-2c_{2}-3\right)\\
    &=& \frac{1}{3}\left(3+6c_{1}+3c_{1}^{2}-2c_{1}^{2}-4c_{1}-2c_{2}-3 \right)\\
    &=& \frac{1}{3}\left(c_{1}^{2}+2c_{1}-2c_{2}\right)
    \end{eqnarray*}
    which yields
    $$|a_{2}^{2}-a_{3}|=\frac{1}{3}\left(|c_{1}^{2}+2c_{1}-2c_{2}|\right)\le \frac{1}{3}\left(|c_{1}|^{2}+2|c_{1}|+2|c_{2}|\right).$$
    It is well-known that $|c_{n}|\le 1-|c_{1}|^{2}$ for $n\ge 2$. Therefore,
    \begin{equation}\label{Abhi-Vasu-P1-equation1}
    |a_{2}^{2}-a_{3}|\le \frac{|c_{1}|^{2}+2|c_{1}|+2(1-|c_{1}|^{2})}{3}=\frac{2|c_{1}|+2-|c_{1}|^{2}}{3}.
    \end{equation}
    Let $x=|c_{1}|$ and $|c_{1}|\le 1$. Let $h(x)=2x+2-x^{2}$, where $0\le x\le 1$. Clearly, $h$ is non-negative in $[0,1]$ and $h'(x)=2-2x$ is also non-negative in $[0,1]$.
    Therefore, $h$ is an increasing function and hence $h(1)=3$ is the maximum value. Hence $2|c_{1}|+2-|c_{1}|^{2}\le 3$. In view of this and  $(\ref{Abhi-Vasu-P1-equation1})$,  we obtain
    $$|a_{2}^{2}-a_{3}|\le 1.
    $$
    A simple computation using  (\ref{Abhi-Vasu-00001})  gives
    \begin{eqnarray*}
    a_{2}a_{3}-a_{4}&=&
    \frac{1}{3}\left((1+c_{1})(2c_{1}^{2}+4c_{1}+2c_{2}+3\right)-\frac{1}{4}\left(2c_{1}^{3}+4c_{1}^{2}+4c_{1}c_{2}+6c_{1}+4c_{2}\right.\\
    &&\qquad\left.+2c_{3}+4\right)\\
    &=&\frac{1}{3}\left(2c_{1}^{3}+6c_{1}^{2}+2c_{1}c_{2}+7c_{1}+2c_{2}+3\right)-\frac{1}{4}\left(2c_{1}^{3}+4c_{1}^{2}+4c_{1}c_{2}+6c_{1}+\right.\\
    &&\qquad\left. 4c_{2}+2c_{3}+4\right)\\
    &=&\frac{1}{12}\left(2c_{1}^{3}+12c_{1}^{2}-4c_{1}c_{2}+10c_{1}-4c_{2}-6c_{3}\right)\\[5mm]
    &=&\frac{1}{6}\left(c_{1}^{3}+6c_{1}^{2}-2c_{1}c_{2}+5c_{1}-2c_{2}-3c_{3}\right).
    \end{eqnarray*}
    Therefore,
    \begin{eqnarray*}
    |a_{2}a_{3}-a_{4}| &\le& \frac{1}{6}\left(|c_{1}|^{3}+6|c_{1}|^{2}+2|c_{1}||c_{2}|+5|c_{1}|+2|c_{2}|+3|c_{3}|\right)\\
    &\le& \frac{1}{6}\left(|c_{1}|^{3}+6|c_{1}|^{2}+2|c_{1}|(1-|c_{1}|^{2})+5|c_{1}|+2(1-|c_{1}|^{2})+3(1-|c_{1}|^{2})\right)\\
    &=& \frac{1}{6}\left(-|c_{1}|^{3}+|c_{1}|^{2}+7|c_{1}|+5\right).
    \end{eqnarray*}
    Thus we obtain
    \begin{equation}\label{Abhi-vasu-P1-equation2}
    |a_{2}a_{3}-a_{4}|\le \frac{1}{6}\left(-|c_{1}|^{3}+|c_{1}|^{2}+7|c_{1}|+5\right).
    \end{equation}
     Let $x=|c_{1}|$, $|c_{1}|\le 1$ and $L(x)=-x^{3}+x^{2}+7x+5$, where $0\le x\le 1$. Then $L'(x)=-3x^{2}+2x+7$. It is easy to see that both the functions
     $L$ and $L'$ are non-negative in $[0,1]$. Therefore,  $L$ is increasing in $[0,1]$ and $L(1)=12$ is the maximum value.
     Hence  $-|c_{1}|^{3}+|c_{1}|^{2}+7|c_{1}|+5\le 12$. Using $(\ref{Abhi-vasu-P1-equation2})$,  we obtain
    $$|a_{2}a_{3}-a_{4}|\le 2.$$
    This completes the proof.
 \end{pf}

\vspace{4mm}
\noindent\textbf{Acknowledgement:}
The first author thank SERB-MATRICS
and the second author thank PMRF-MHRD, Govt. of India for their support.