Paper deep dive
Breach in the Shield: Unveiling the Vulnerabilities of Large Language Models
Runpeng Dai, Run Yang, Fan Zhou, Hongtu Zhu
Models: Llama, Pythia-160M, Pythia-1B, Pythia-410M, Pythia-70M, Qwen2, Qwen-VL
Intelligence
Status: succeeded | Model: google/gemini-3.1-flash-lite-preview | Prompt: intel-v1 | Confidence: 95%
Last extracted: 3/12/2026, 6:12:52 PM
Summary
The paper introduces 'FI' (First order local Influence), a stability measure based on information geometry to quantify the sensitivity of LLMs and VLMs to input and parameter perturbations. The study demonstrates that a small subset of high-FI parameters or input dimensions significantly contributes to model brittleness, and mitigating these vulnerabilities improves model robustness during merging.
Entities (5)
Relation Signals (4)
High FI values ā contributesto ā Model Brittleness
confidence 95% Ā· A small subset of parameters or input dimensions with high FI values disproportionately contribute to model brittleness.
FI ā quantifiessensitivityof ā LLM
confidence 95% Ā· FI... quantifies the sensitivity of individual parameter and input dimensions.
FI ā quantifiessensitivityof ā VLM
confidence 95% Ā· Our extensive analysis across LLMs and VLMs... reveals that... FI effectively identifies vulnerabilities.
FI ā rootedin ā Information Geometry
confidence 90% Ā· FI, First order local Influence, which is rooted in information geometry
Cypher Suggestions (0)
No Cypher suggestions yet.
Abstract
Abstract:Large Language Models (LLMs) and Vision-Language Models (VLMs) have achieved impressive performance across a wide range of tasks, yet they remain vulnerable to carefully crafted perturbations. In this study, we seek to pinpoint the sources of this fragility by identifying parameters and input dimensions (pixels or token embeddings) that are susceptible to such perturbations. To this end, we propose a stability measure called \textbf{FI}, \textbf{F}irst order local \textbf{I}nfluence, which is rooted in information geometry and quantifies the sensitivity of individual parameter and input dimensions. Our extensive analysis across LLMs and VLMs (from 1.5B to 13B parameters) reveals that: (I) A small subset of parameters or input dimensions with high FI values disproportionately contribute to model brittleness. (II) Mitigating the influence of these vulnerable parameters during model merging leads to improved performance.
Tags
Links
- Source: https://arxiv.org/abs/2504.03714
- Canonical: https://arxiv.org/abs/2504.03714
Trouble viewing inline? Open PDF directly ā
Full Text
88,548 characters extracted from source content.
Expand or collapse full text
Breach in the Shield: Unveiling the Vulnerabilities of Large Language Models Runpeng Dai1ā Run Yang2ā Fan Zhou3 Hongtu Zhu1ā 1University of North Carolina at Chapel Hill 2BiliBili 3Shanghai University of Finance and Economics runpeng, htzhu@email.unc.edu yangrun@bilibili.com zhoufan@mail.shufe.edu.cn Abstract Large Language Models (LLMs) and Vision-Language Models (VLMs) have achieved impressive performance across a wide range of tasks, yet they remain vulnerable to carefully crafted perturbations. In this study, we seek to pinpoint the sources of this fragility by identifying parameters and input dimensions (pixels or token embeddings) that are susceptible to such perturbations. To this end, we propose a stability measure called FI, First order local Influence, which is rooted in information geometry and quantifies the sensitivity of individual parameter and input dimensions. Our extensive analysis across LLMs and VLMs (from 1.5B to 13B parameters) reveals that: (I) A small subset of parameters or input dimensions with high FI values disproportionately contribute to model brittleness. (I) Mitigating the influence of these vulnerable parameters during model merging leads to improved performance. Breach in the Shield: Unveiling the Vulnerabilities of Large Language Models Runpeng Dai1ā Run Yang2ā Fan Zhou3 Hongtu Zhu1ā 1University of North Carolina at Chapel Hill 2BiliBili 3Shanghai University of Finance and Economics runpeng, htzhu@email.unc.edu yangrun@bilibili.com zhoufan@mail.shufe.edu.cn 1 Introduction 111āBoth authors contributed equally to this work. Large Language Models (LLMs) and Vision Language Models (VLMs) such as GPT (Brown et al., 2020) and Llama (Touvron et al., 2023), have revolutionized the field of Natural Language Processing (NLP), exhibiting remarkable proficiency across a variety of tasks (Gong et al., 2024; Zheng et al., 2025b; Luo et al., 2025) and modalities (Bai et al., 2023; Liu et al., 2024b; Zheng et al., 2025a). These modern LLMs are massive in size, trained on vast amounts of data, and meticulously aligned to prevent from generating harmful content (Perez et al., 2022), leaking private information (Zhang et al., 2024), or exhibiting sexual or religious bias (Xie and Lukasiewicz, 2023). Despite the enthusiasm for these integrative approaches, a critical issue remains: LLMs remain susceptible to both external and internal perturbations, affecting their reliability and performance. Externally, LLMs are vulnerable to input perturbations, such as Embedding-Corrupted Prompts (Fort, 2023; Liu et al., 2024a). This susceptibility extends to visual inputs in VLMs, where adversarially optimized images can drastically alter model behavior (Qi et al., 2024). Beyond adversarial attacks, VLMs exhibit high sensitivity to perturbations in specific local regions of an imageāa common issue, as user-uploaded images often suffer from blurring, masking, or low resolution. The vulnerability is highlighted in our case study of the Qwen-VL model. As depicted in Figure 1, masking the ten most sensitive pixels, which are unrelated to the question, resulted in incorrect model outputs. Figure 1: A case study of the Qwen-VL model (Bai et al., 2023) on SCI-QA. The image on the far right visualizes the per-pixel FI values. Masking just 10 pixels with the highest FI values leads to a failure in producing the correct answer. Internally, LLM stability is further challenged by parameter perturbations, often introduced through model merging and quantization. While these techniques improve deployment efficiency by reducing inference costs (Frantar and Alistarh, 2023; Ashkboos et al., 2024), they can also induce hallucinations and degrade performance (Men et al., 2024; Yu et al., 2024; Li et al., 2024). However, our findings reveal that parameter susceptibility varies significantly. As Figure 2 illustrates, randomly dropping 5% of parameters has a minimal impact on performance. In contrast, zeroing out just 1% of the parameters identified by our measure can drastically reduce accuracy, even below random guessing levels. Figure 2: A case study of Qwen2 on MMLU-geography. This figure illustrates the decline in response accuracy as a greater percentage of parameters are zeroed. āFI-Highā targets parameters with the highest FI values, while āRandomā indicates random parameter removal. To pinpoint the sources of this fragility, we propose a novel stability measure called FI, First order local Influence, to quantitatively assess the stability of LLMs against perturbations. Specifically, we construct a perturbation manifold that encompasses all perturbed models, along with its associated geometric properties. Our stability measure quantifies the degree of local influence of a perturbation on a given objective function within this manifold, thereby reflecting the stability of individual LLM components. FIās versatility allows for effective stability assessment under both external and internal perturbations across various granularitiesāfrom individual parameters to input features like pixels and patches. FI effectively identifies vulnerabilities. Our extensive studies validate its effectiveness in pinpointing fragile pixels in VLM vision inputs, vulnerable embedding dimensions of tokens in LLMs (Section 4), and salient model parameters (Subsection 5.1). FI offers insights into improving model robustness. We further illustrate that understanding these vulnerabilities can lead to enhanced model resistance to perturbations. By focusing on model merging as an example, we show that safeguarding key parameters identified by high FI values can substantially reduce performance degradation during the merging process (Subsection 5.2). 2 Related Work Recent efforts to evaluate LLM stability typically adopt a coarse-grained approach, aiming to assess the overall robustness of models under various perturbations. One line of work investigates how stability is influenced by sampling parameters, such as temperature, which affect output variability during generation (Atil et al., 2024; Ouyang et al., 2025). Another direction studies model sensitivity to input or parameter perturbations. For instance, Liu et al. analyze input-level robustness using optimal transport to quantify a modelās response to distributional shifts in prompts. On the other hand, Peng et al. (2024) focus on parameter-space perturbations, demonstrating that LLMs remain robust to weight changes up to a certain threshold, beyond which performance significantly degrades. They estimate a modelās robustness tolerance by injecting random perturbations into model weights and evaluating the performance drop. Despite these contributions, fine-grained analysesāsuch as those examining the effect of individual input tokens, pixels, or specific model parametersāremain underexplored. Wei et al. (2024) take a step in this direction by leveraging pruning-based techniques, including SNIP (Lee et al., 2018) and Wanda (Sun et al., 2023), to identify critical neurons and low-rank structures that impact model safety and utility. However, there is still a lack of unified metrics or frameworks that assess stability with respect to both input- and parameter-level perturbations. Moreover, the downstream applications of such stability assessments remain largely unexamined. 3 Stability Measure of Large language models In this section, we propose a new metric called FI to quantify the stability of large language models against local perturbations. Considering the auto-regressive nature of LLMs, we first develop FI for single-step generation and discuss its theoretical and computational properties in detail. We then show how FI can be naturally extended to sequence generation tasks. Finally, we compare FI to existing stability measures, highlighting its unique advantages. 3.1 FI Metric Problem formulation. Consider an LLM parameterized by Īø, with input data x, which may consist of text or, for visual language models, a combination of text and images. Given x, the model generates a probability distribution over its vocabulary to predict the next token, which can be framed as a classification problem with K classes, where K represents the vocabulary size. However, vocabulary sizes are typically large (Bai et al., 2023; Dubey et al., 2024), and predictions are often concentrated on a small subset of tokens. Instead of using the entire vocabulary, it is more efficient to focus on a relevant subset based on the task. For example, in multiple-choice questions, probabilities are restricted to the choices "A", "B", "C", or "D". Classes can also be defined semantically, such as categorizing tokens as "neutral" or "notorious" in toxicity detection (Gehman et al., 2020). With appropriately defined classes, the predicted probability for class yā1,ā¦,K1ā¦yā\1,ā¦,K\y ā 1 , ⦠, K is denoted as Pā¢(y|x,Īø)conditionalP(y|x,Īø)P ( y | x , Īø ), satisfying āy=1KPā¢(y|x,Īø)=1superscriptsubscript1conditional1 _y=1^KP(y|x,Īø)=1āy = 1K P ( y | x , Īø ) = 1. Let ĻāādsuperscriptāĻ ^dĻ ā blackboard_Rd be a perturbation vector varies in an open subset Ī© Ī©. Ļ can be applied to a subset of the model parameters Īø and locations within the input data x. We denote the output of the perturbed model under this perturbation as Pā¢(y|x,Īø,Ļ)conditionalP(y|x,Īø,Ļ)P ( y | x , Īø , Ļ ). Perturbation Manifold and FI Since our primary interest lies in examining the behavior of Pā¢(y|x,Ļ,Īø)conditionalP(y|x,Ļ,Īø)P ( y | x , Ļ , Īø ) as a function of Ļ near Ļ0=0subscript00 _0=0Ļ0 = 0, we shift focus from Īø to Ļ. We introduce the perturbation manifold as defined in Zhu et al. (2007) and Zhu et al. (2011). Definition 3.1. Define the d-dimensional perturbation manifold ā³=P(y|x,Īø,Ļ):ĻāĪ©M=\P(y|x,Īø,Ļ):Ļā \M = P ( y | x , Īø , Ļ ) : Ļ ā Ī© , which encompasses all perturbed models. Assume that for all ĻāĪ©Ļā Ļ ā Ī©, the perturbed models Pā¢(y=i|x,Īø,Ļ)i=1Ksuperscriptsubscriptconditional1\P(y=i|x,Īø,Ļ)\_i=1^K P ( y = i | x , Īø , Ļ ) i = 1K are positive and sufficiently smooth. The tangent space TĻsubscriptT_ĻTitalic_Ļ of ā³MM at Ļ is spanned by the partial derivatives of the log-likelihood function āā¢(Ļ|y,x,Īø)=logā”Pā¢(y|x,Īø,Ļ)āconditionalconditional (Ļ|y,x,Īø)= P(y|x,Īø,Ļ)ā ( Ļ | y , x , Īø ) = log P ( y | x , Īø , Ļ ) with respect to Ļ, specifically TĻ=spanā¢āĻiā¢āā¢(Ļ|y,x,Īø)i=1dsubscriptspansuperscriptsubscriptsubscriptāconditional1T_Ļ=span\ ā _i (Ļ|y,x,% Īø)\_i=1^dTitalic_Ļ = span divide start_ARG ā end_ARG start_ARG ā Ļitalic_i end_ARG ā ( Ļ | y , x , Īø ) i = 1d. The metric gĻsubscriptg_Ļgitalic_Ļ on ā³MM can be defined with the metric tensor GĻsubscriptG_ĻGitalic_Ļ. Consider two tangent vectors at Ļ given by vjā¢(Ļ)=hjā¤ā¢āĻāā¢(Ļ|y,x,Īø)āTĻsubscriptsuperscriptsubscriptātopsubscriptāconditionalsubscriptv_j(Ļ)=h_j _Ļ (Ļ|y,x,Īø)ā T_Ļvitalic_j ( Ļ ) = hitalic_j⤠āĻ ā ( Ļ | y , x , Īø ) ā Titalic_Ļ, where h1subscriptā1h_1h1 h2subscriptā2h_2h2 are the weights on the basis. Their inner product is defined as: āØv1ā¢(Ļ),v2ā¢(Ļ)ā©gĻ=āy=1Kv1ā¢(Ļ)ā¢v2ā¢(Ļ)ā¢Pā¢(y|x,Īø).subscriptsubscript1subscript2subscriptsuperscriptsubscript1subscript1subscript2conditional v_1(Ļ),v_2(Ļ) _g_Ļ= _y=1^Kv_1(% Ļ)v_2(Ļ)P(y|x,Īø).⨠v1 ( Ļ ) , v2 ( Ļ ) ā©g start_POSTSUBSCRIPT Ļ end_POSTSUBSCRIPT = āy = 1K v1 ( Ļ ) v2 ( Ļ ) P ( y | x , Īø ) . The metric tensor GĻsubscriptG_ĻGitalic_Ļ is given by: GĻ=āy=1KāĻāā¢(Ļ|y,x,Īø)ā¢āĻā¤āā¢(Ļ|y,x,Īø)ā¢Pā¢(y|x,Īø,Ļ).subscriptsuperscriptsubscript1subscriptāconditionalsubscriptsuperscripttopāconditionalconditionalG_Ļ= _y=1^K _Ļ (Ļ|y,x,Īø)ā^% _Ļ (Ļ|y,x,Īø)P(y|x,Īø,Ļ).Gitalic_Ļ = āy = 1K āĻ ā ( Ļ | y , x , Īø ) āā¤Ļ ā ( Ļ | y , x , Īø ) P ( y | x , Īø , Ļ ) . Subsequently, the norm of vjā¢(Ļ)subscriptv_j(Ļ)vitalic_j ( Ļ ) under metric gĻsubscriptg_Ļgitalic_Ļ is āvjāgĻ=hjā¤ā¢GĻā¢hjsubscriptnormsubscriptsubscriptsuperscriptsubscriptātopsubscriptsubscriptā\|v_j\|_g_Ļ= h_j G_Ļh_jā„ vitalic_j ā„g start_POSTSUBSCRIPT Ļ end_POSTSUBSCRIPT = square-root start_ARG hitalic_j⤠Gitalic_Ļ hitalic_j end_ARG. Let Cā¢(t)=Pā¢(y|x,Īø,Ļā¢(t))conditionalC(t)=P(y|x,Īø,Ļ(t))C ( t ) = P ( y | x , Īø , Ļ ( t ) ) be a smooth curve on the manifold ā³MM connecting two points Ļ1=Ļā¢(t1)subscript1subscript1 _1=Ļ(t_1)Ļ1 = Ļ ( t1 ) and Ļ2=Ļā¢(t2)subscript2subscript2 _2=Ļ(t_2)Ļ2 = Ļ ( t2 ). Then, the distance between Ļ1subscript1 _1Ļ1 and Ļ2subscript2 _2Ļ2 along the curve Cā¢(t)C(t)C ( t ) is given by: SCā¢(Ļ1,Ļ2)subscriptsubscript1subscript2 S_C (Ļ_1,Ļ_2 )Sitalic_C ( Ļ1 , Ļ2 ) =ā«t1t2ā„ātlogP(y|x,Īø,Ļ(t))ā„gĻā¢t = _t_1^t_2 \| _t P(y|x,Īø,Ļ% (t))\|_g_Ļ\,dt= ā«t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPTt2 square-root start_ARG ā„ āt log P ( y | x , Īø , Ļ ( t ) ) ā„g start_POSTSUBSCRIPT Ļ end_POSTSUBSCRIPT end_ARG d t =ā«t1t2dā¢Ļā¢(t)Tdā¢tā¢GĻā¢(t)ā¢dā¢Ļā¢(t)dā¢tā¢t.absentsuperscriptsubscriptsubscript1subscript2superscriptsubscriptdifferential-d = _t_1^t_2 dĻ(t)^TdtG_Ļ(t% ) dĻ(t)dt\,dt.= ā«t start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPTt2 square-root start_ARG divide start_ARG d Ļ ( t )T end_ARG start_ARG d t end_ARG Gitalic_Ļ ( t ) divide start_ARG d Ļ ( t ) end_ARG start_ARG d t end_ARG end_ARG d t . With the Perturbation manifold ā³MM and respective metric gĻsubscriptg_Ļgitalic_Ļ defined, we are ready to propose the metric that quantifies the stability of large language models (LLMs) against various types of local perturbations. Let fā¢(Ļ)f(Ļ)f ( Ļ ) be the objective function of interest for sensitivity analysis, in our case being ālogā”Pā¢(ypā¢rā¢eā¢d|x,Īø,Ļ)conditionalsubscript- P(y_pred|x,Īø,Ļ)- log P ( yitalic_p r e d | x , Īø , Ļ ), we can define the following (first-order) local influence metric FI: Definition 3.2. Given the perturbation manifold ā³MM and its metric, the first-order local stability measure of fā¢(Ļ)f(Ļ)f ( Ļ ) at Ļā¢(0)=Ļ00subscript0Ļ(0)= _0Ļ ( 0 ) = Ļ0 is defined as Ļā¢(Ļ0)=maxCā¢limtā0[fā¢(Ļā¢(t))āfā¢(Ļā¢(0))]2SC2ā¢(Ļā¢(t),Ļā¢(0)).subscriptsubscript0subscriptsubscriptā0superscriptdelimited-[]02superscriptsubscript20FI_Ļ ( _0 )= _C _tā 0% [f(Ļ(t))-f(Ļ(0))]^2S_C^2(Ļ(t),Ļ(0)).FIitalic_Ļ ( Ļ0 ) = maxitalic_C limitalic_t ā 0 divide start_ARG [ f ( Ļ ( t ) ) - f ( Ļ ( 0 ) ) ]2 end_ARG start_ARG Sitalic_C2 ( Ļ ( t ) , Ļ ( 0 ) ) end_ARG . (1) The ratio in Equation 1 measures the amount of change introduced to the objective function relative to the distance of the perturbation on the perturbation manifold. Thus, Equation 1 can be naturally interpreted as the maximum local ratio of change among all possible perturbation curves Cā¢(t)C(t)C ( t ). Computation of FI. As we will show, Theorem A.1 on diffeomorphic reparameterization invariance enables us to derive an easy-to-compute solution for Equation 1, while addressing the low-dimensionality problem inherent in LLMs. Theorem 3.3. If GĻsubscriptG_ĻGitalic_Ļ is positive definite, the FI measure has the following closed-form: Ļā¢(Ļ0)=āfā¢(Ļ0)TGĻ0ā1ā¢āfā¢(Ļ0),subscriptsubscript0superscriptsubscriptāsubscript0subscriptsuperscript1subscript0subscriptāsubscript0FI_Ļ (Ļ_0 )=ā_f (% Ļ_0 )^TG^-1_ _0ā_f (Ļ_0 ),FIitalic_Ļ ( Ļ0 ) = āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT G- 1Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT , (2) where āfā¢(Ļ0)=āfā¢(Ļ)āĻ|Ļ=Ļ0.subscriptāsubscript0evaluated-atsubscript0ā_f (Ļ_0 )= ā f(Ļ)ā% Ļ |_Ļ=Ļ_0.āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT = divide start_ARG ā f ( Ļ ) end_ARG start_ARG ā Ļ end_ARG |Ļ = Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT . The detailed proof of Theorem 3.3 can be found in Appendix A.5. It is important to note that the closed form of FI in Theorem 3.3 depends on the positive definiteness of GĻsubscriptG_ĻGitalic_Ļ, which is not always guaranteed. This is due to the fact that the parameters in LLMs are often high-dimensional tensors with low-rank structures (Kaushal et al., 2023). We apply the invariant Theorem A.1 by transforming Ļ to a vector ν such that Gν=KsubscriptsubscriptG_ν=I_KGitalic_ν = Iitalic_K, where K is an integer. Specifically, we notice that GĻ0=B0Tā¢B0subscriptsubscript0superscriptsubscript0subscript0G_ _0=B_0^TB_0Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = B0italic_T B0, where B0=[Pā¢(y=i|x,Īø,Ļ)1/2ā¢āĻāā¢(Ļ|y=i,x,Īø)]i⩽K.subscript0subscriptdelimited-[]superscriptconditional12subscriptāconditionalB_0= [P(y=i|x,Īø,Ļ)^1/2 _Ļ (Ļ|y=i,x,% Īø) ]_i K.B0 = [ P ( y = i | x , Īø , Ļ )1 / 2 āĻ ā ( Ļ | y = i , x , Īø ) ]i ⩽ K . Let r0=rankā¢(GĻ0)subscript0ranksubscriptsubscript0r_0=rank(G_ _0)r0 = rank ( Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ), we apply the compact SVD to B0āāpĆKsubscript0superscriptāB_0 ^pĆ KB0 ā blackboard_Rp Ć K, which yields B0=V0ā¢Ī0ā¢U0subscript0subscript0subscriptĪ0subscript0B_0=V_0 _0U_0B0 = V0 Ī0 U0, where V0āāpĆr0subscript0superscriptāsubscript0V_0 ^pĆ r_0V0 ā blackboard_Rp Ć r0 and U0āār0ĆKsubscript0superscriptāsubscript0U_0 ^r_0Ć KU0 ā blackboard_Rr0 Ć K are semi-orthogonal matrices and Ī0āār0Ćr0subscriptĪ0superscriptāsubscript0subscript0 _0 ^r_0Ć r_0Ī0 ā blackboard_Rr0 Ć r0 is a diagonal matrix. Under the transformation ν=Ī0ā¢V0Tā¢ĻsubscriptĪ0superscriptsubscript0ν= _0V_0^TĻν = Ī0 V0italic_T Ļ, we have Ļā¢(Ļ0)=νā¢(ν0)subscriptsubscript0subscriptsubscript0FI_Ļ (Ļ_0 )=FI_ν (ν% _0 )FIitalic_Ļ ( Ļ0 ) = FIitalic_ν ( ν0 ) , which can be expressed as āfā¢(Ļ0)ā¤(V0R0)ā¤Ī0ā2(V0R0)āfā¢(Ļ0), ā_f (Ļ_0 ) (V_0R_0% ) _0^-2 (V_0R_0 )ā_f (% Ļ_0 ),āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT⤠( V0 R0 )⤠Ī0- 2 ( V0 R0 ) āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT , where the equality holds by applying the chain rule to GνsubscriptG_νGitalic_ν. FI for sequence generation. Sequence generation is essentially multiple rounds of next-token generation, where the l-th token y(l)superscripty^(l)y( l ) is generated given the initial input z and previously generated tokens (l)=y(1),ā¦,y(lā1)superscriptsuperscript1ā¦superscript1 y^(l)=\y^(1),ā¦,y^(l-1)\italic_y( l ) = y( 1 ) , ⦠, y( l - 1 ) . We define the FI measure for generating the l-th token y(l)superscripty^(l)y( l ) given the initial input z by averaging out the randomness from the preceding steps lā¢(z)=(l)ā¢[ā¢(z,(l),Īø,Ļ)|z]subscriptsubscriptsuperscriptdelimited-[]conditionalsuperscriptFI_l(z)=E_ y^(l)[FI(\z, y^(l)\,% Īø,Ļ)|z]FIitalic_l ( z ) = blackboard_Eitalic_y( l ) [ FI ( z , italic_y( l ) , Īø , Ļ ) | z ]. To formulate an overall measure for sequence generation, we aggregate these per-token FI measures. Since sequences generated by LLMs can vary in length, we propose two methods to handle this heterogeneity. The first approach sets a fixed horizon L and computes the mean FI over these rounds seqLā¢(z)=1Lā¢āl=1Llā¢(z).subscriptsuperscriptseq1superscriptsubscript1subscriptFI^L_seq(z)= 1L _l=1^LFI_l(z).FIitalic_Lseq ( z ) = divide start_ARG 1 end_ARG start_ARG L end_ARG āl = 1L FIitalic_l ( z ) . (3) Alternatively, inspired by the concept of average discounted rewards in reinforcement learning (Liu et al., 2018), we consider sequences of potentially infinite length and propose a discounted FI measure with discount factor γ seqā,γā¢(z)=(1āγ)ā¢āl=0āγlā lā¢(z).subscriptsuperscriptseq1superscriptsubscript0ā superscriptsubscriptFI^ā,γ_seq(z)=(1-γ) _l=0^ā% γ^lĀ·FI_l(z).FIā , γseq ( z ) = ( 1 - γ ) āl = 0ā γitalic_l ā FIitalic_l ( z ) . By taking the expectation over the distribution of z, we obtain the average FI for sequence generation in both cases Pzā¢[seqLā¢(z)]subscriptsubscriptdelimited-[]subscriptsuperscriptseqE_P_z[FI^L_seq(z)]blackboard_EP start_POSTSUBSCRIPT z end_POSTSUBSCRIPT [ FIitalic_Lseq ( z ) ] and Pzā¢[seqā,γā¢(z)]subscriptsubscriptdelimited-[]subscriptsuperscriptseqE_P_z[FI^ā,γ_seq(z)]blackboard_EP start_POSTSUBSCRIPT z end_POSTSUBSCRIPT [ FIā , γseq ( z ) ], respectively. 3.2 Other Measures & Discussion We note that several alternative methods can also serve as stability measures for LLMs. We provide their explicit formulations and compare them with FI. Jacobian Norm (Novak et al., 2018): āāĻfā¢(ypā¢rā¢eā¢d,Ļ)ā2subscriptnormsubscriptsubscript2\| _Ļf(y_pred,Ļ)\|_2ā„ āĻ f ( yitalic_p r e d , Ļ ) ā„2 SNIP (Lee et al., 2018): āĻāāĻfā¢(ypā¢rā¢eā¢d,Ļ)ā2subscriptnormdirect-productsubscriptsubscript2\|Ļ _Ļf(y_pred,Ļ)\|_2ā„ Ļ ā āĻ f ( yitalic_p r e d , Ļ ) ā„2 Both measures focuses solely on ypā¢rā¢eā¢dsubscripty_predyitalic_p r e d, while neglecting the probabilities assigned to other choices. For example, consider two output distributions: (0.9, 0.05, 0.05, 0.02) and (0.3, 0.25, 0.25, 0.2). In both cases, the model selects option A. However, the second distribution is more unstable, as a small perturbation in the probabilities could lead to a different prediction. In contrast, FI measure accounts for both the probability and gradient across all possible choices. Saliency map (Simonyan et al., 2013): 0ā¢ifā¢āfā¢(ypā¢rā¢eā¢d,Ļ)āĻā¢<0⢠or ā¢āyā ypā¢rā¢eā¢dāfā¢(y,Ļ)āĻ>ā¢0āāfā¢(ypā¢rā¢eā¢d,Ļ)āĻā¢āyā ypā¢rā¢eā¢dāfā¢(y,Ļ)āĻcases0ifsubscriptexpectation0 or subscriptsubscript0otherwisesubscriptsubscriptsubscriptotherwise cases0~~if~ ā f(y_pred,Ļ)ā% Ļ<0 or _yā y_pred ā f(y,Ļ)ā% Ļ>0\\ - ā f(y_pred,Ļ)āĻ _yā y_pred % ā f(y,Ļ)āĻ cases start_ROW start_CELL 0 if divide start_ARG ā f ( yitalic_p r e d , Ļ ) end_ARG start_ARG ā Ļ end_ARG < 0 or āy ā y start_POSTSUBSCRIPT p r e d end_POSTSUBSCRIPT divide start_ARG ā f ( y , Ļ ) end_ARG start_ARG ā Ļ end_ARG > 0 end_CELL start_CELL end_CELL end_ROW start_ROW start_CELL - divide start_ARG ā f ( yitalic_p r e d , Ļ ) end_ARG start_ARG ā Ļ end_ARG āy ā y start_POSTSUBSCRIPT p r e d end_POSTSUBSCRIPT divide start_ARG ā f ( y , Ļ ) end_ARG start_ARG ā Ļ end_ARG end_CELL start_CELL end_CELL end_ROW Saliency map take the gradient of all choices into account. But loose too much information by zeroing out many of the gradient. Although all these methods can be used to assess the vulnerability of LLMs, we highlight FI for its distinct advantages. Effectiveness: A quantitative comparison of these measures is provided in Section 4, while their computational complexities are discussed in Appendix A.2. Theoretical rigor: In particular, only FI possesses a reparameterization invariance property (see Appendix A.3), which further distinguishes it by enhancing interpretability. 4 External perturbations analysis In this section, we first demonstrate the effectiveness of FI in identifying vulnerable locations in both vision and language inputs though guided attack. Then, we conclude the section with a finding from cross-modal analysis. Identify Fragile Pixels We conduct the attack process on the MMbench dataset Liu et al. (2024c), a comprehensive benchmark designed to evaluate various multi-modal capabilities of VLMs. For a fair comparison, we identify the top 10 pixels using different stability measures and assess the modelās performance after masking out the corresponding pixels. Table 1: Accuracy on the MMBench dataset after masking out top ten pixels in images identified by different measures. The first block shows summary accuracy across all models and methods. Model Method Action Recognition Attribute Recognition Celebrity Recognition Function Reasoning Qwen VL FI (Ours) 0.320 0.402 0.673 0.411 Jacobian 0.668 0.587 0.906 0.604 Saliency 0.782 0.525 0.873 0.639 Random 0.812 0.550 0.881 0.683 Original 0.814 0.549 0.882 0.686 Qwen2.5 VL-3B FI (Ours) 0.720 0.735 0.780 0.723 Jacobian 0.731 0.752 0.797 0.755 Saliency 0.745 0.761 0.797 0.774 Random 0.882 0.931 0.957 0.928 Original 0.890 0.946 0.959 0.930 Qwen2.5 VL-7B FI (Ours) 0.768 0.750 0.796 0.723 Jacobian 0.778 0.768 0.815 0.755 Saliency 0.792 0.777 0.815 0.774 Random 0.891 0.944 0.951 0.925 Original 0.890 0.946 0.959 0.930 Identify Vulnerable Embedding Dimensions We conduct attack on pure-text LLMs to verify the effectiveness of our approach in identifying vulnerable embedding dimensions. Specifically, we follow the token embedding attack methods proposed in Liu et al. (2024a) and Fort (2023). More concretely, we compute the stability measure for each embedding dimension and select the top 0.1% most sensitive dimensions(Ļ) as identified by the metrics. We then apply a gradient-based attack strategy following Fort (2023), perturbing the selected dimensions in the direction of āāĻlogā”Pā¢(ypredā£x,Īø)subscriptāconditionalsubscriptpred- _Ļ P(y_pred x,Īø)- āĻ log P ( ypred ⣠x , Īø ). From both Table 1 and Table 2, we observe the following: (I) Stability measures are effective in identifying vulnerable input dimensions (i.e., pixels in images and dimensions in embeddings). Notably, LLMs are generally robust to random perturbations and such perturbations rarely lead to significant performance degradation. In contrast, perturbations guided by stability measures consistently result in substantial drops in performance. (I) Among all the stability measures evaluated, FI proves to be the most effective: masking pixels or perturbing dimensions identified by FI leads to the largest observed decline in performance. Table 2: Comparison of accuracy in the MMLU dataset after perturbing the same number of dimensions in the embedding space identified using different measures. Model Method Business Geo Culture Law Pythia 1B Saliency 0.278 0.272 0.210 0.243 Jacobian 0.273 0.264 0.201 0.241 Random 0.301 0.368 0.237 0.246 FI (ours) 0.270 0.261 0.195 0.236 SNIP 0.297 0.281 0.226 0.242 Original 0.303 0.370 0.240 0.247 Qwen2.5 3B Saliency 0.677 0.637 0.632 0.560 Jacobian 0.665 0.641 0.625 0.560 Random 0.805 0.781 0.781 0.672 FI (ours) 0.656 0.620 0.610 0.547 SNIP 0.783 0.663 0.665 0.563 Original 0.810 0.800 0.785 0.673 Qwen2.5 7B Saliency 0.756 0.789 0.709 0.725 Jacobian 0.764 0.782 0.717 0.720 Random 0.852 0.884 0.802 0.735 FI (ours) 0.748 0.780 0.705 0.713 SNIP 0.757 0.791 0.710 0.727 Original 0.856 0.890 0.810 0.737 Figure 3: A case study utilizing FI for cross-modal analysis. In the same example, the bottom-left image shows how Aggresive and Safe prompts affects the FI distribution on the image. Effect of Prompting on Pixel Vulnerability While the significant impact of prompt design on VLM performance is well-recognized (Zhou et al., 2022), and carefully crafted prompts are known to even jailbreak these models (Shayegani et al., 2023), a quantitative analysis of this cross-modal influence ā specifically, how prompting affects the processing and stability of visual input ā remains largely unexplored. Our study aims to bridge this gap by investigating how varying prompt instructions influence the sensitivity of VLMs to visual perturbations. Specifically, we examine two types of prompts: ⢠Aggressive Prompts: Designed to encourage the model to consider every detail in the image, potentially increasing sensitivity to noise. ⢠Safe Prompts: Intended to focus the model on salient entities and relationships, potentially enhancing robustness by ignoring irrelevant details. We computed the FI value for each pixel and visualized the resulting distributions under different prompt settings, as illustrated in Figure 3. Our main findings are as follows: (I) Prompt choice has a substantial impact on the stability of individual pixels within the image. As shown on the left of Figure 3, aggressive prompts shift the FI distribution toward higher values, resulting in a marked increase in both the mean and maximum FI values. This suggests that the model becomes more sensitive to pixel-level perturbations throughout the image. In contrast, safe prompts significantly shift the FI distribution toward lower values, indicating reduced sensitivity and improved stability against perturbations in less relevant regions. (I) Vulnerability remains even with careful prompt design. Although safe prompts generally reduce FI values, they do not fully guarantee model stability, as outliers with large FI values persist. As shown in the right column of Figure 3, even when applying the safe prompt, masking out the two pixels with the highest FI values still leads to incorrect model predictions. This result underscores the persistent challenge of achieving robustness in VLMs and demonstrates the effectiveness of the FI measure for identifying vulnerable regions. Our findings contribute to the growing body of literature on cross-modal interactions in VLMs, offering a stability-centric perspective that complements existing behavioral and attributional analyses. Importantly, this framework can inform the development of more robust multimodal systems and prompt design strategies for safety-critical applications. 5 Internal perturbations analysis In this section, we first conduct an parameter sparsification experiment to demonstrate the effectiveness of the FI. We then apply the FI measure to mitigate parameter interference during model merging, showcasing its potential for guiding LLM improvement. 5.1 Parameter sparsification (a) business (b) geography (c) culture (d) law Figure 4: Comparing the accuracy in the MMLU dataset of Qwen2-7B when parameters have been sparsified at different rates. We conduct experiments on multiple-choice problems from the MMLU dataset (Hendrycks et al., 2020) and sequence generation tasks from Alpaca-Eval (Dubois et al., 2024) to examine how these perturbations impact two key capabilities of large models: knowledge retention and instruction-following. Details of both experimental setups are provided in Figure A.1. As shown in Figure 2 and Figure 4 , sparsifying just 2ā3% of the high-FI parameters significantly degrades the modelās knowledge capacity, leading to catastrophic forgetting and hallucinations, with performance dropping by up to 75%. A similar trend is observed in Table 6 at around the 10% sparsity level. In contrast, models remain relatively robust against random sparsification, often exhibiting nearly identical behavior even after 5% sparsification. These findings demonstrate FIās effectiveness in identifying fragile parameters and further support the inherent structure within the parameter matrix, aligning with recent observations on model brittleness (Ma et al., 2023; Wei et al., 2024; Yu et al., 2024). 5.2 FI-Guided Parameter Protection in Model Merging Model merging is a technique for acquiring domain-specific knowledge by combining models from different domains, thereby reducing the computational cost of additional fine-tuning (see (Yang et al., 2024c) for a review). However, a persistent challenge is that merging parameters introduces perturbations that can hinder a modelās ability to retain previously learned information. To address this, we use FI to identify parameters susceptible to forgetting and exclude them from the merging process. We demonstrate that FI can be seamlessly integrated into mainstream model merging methods, including Average Merging (Wortsman et al., 2022), Task Arithmetic (Ilharco et al., 2022), and TIES (Yadav et al., 2024). Additionally, we include DARE (Yu et al., 2024) as a competing baseline for completeness. We consider merging two models, A and B, both fine-tuned from the same base model. Let ĪøAsubscript _AĪøitalic_A, ĪøBsubscript _BĪøitalic_B, and ĪøBasesubscriptBase _BaseĪøBase denote the parameters of models A, B, and the base model, respectively. We first introduce the merging methods and then demonstrate how FI can be integrated into these algorithms to mitigate perturbation effects. Average Merging Average merging obtains the merged model by averaging ĪøAsubscript _AĪøitalic_A and ĪøBsubscript _BĪøitalic_B, resulting in parameters ĪøAvg=ĪøA+ĪøB2subscriptAvgsubscriptsubscript2 _Avg= _A+ _B2ĪøAvg = divide start_ARG Īøitalic_A + Īøitalic_B end_ARG start_ARG 2 end_ARG. Task Arithmetic Task arithmetic constructs ātask vectorsā by subtracting a base model from each task-specific model and then merges these vectors linearly before adding back the base model ĪøTask=ĪøBase+γā¢(Ī“A+Ī“B)subscriptTasksubscriptBasesubscriptsubscript _Task= _Base+γ( _A+ _B)ĪøTask = ĪøBase + γ ( Ī“italic_A + Ī“italic_B ), where Ī“A=ĪøAāĪøBasesubscriptsubscriptsubscriptBase _A= _A- _BaseĪ“italic_A = Īøitalic_A - ĪøBase and similarly for Ī“Bsubscript _BĪ“italic_B. Both Average Merging and Task Arithmetic modify all parameters in models A and B, potentially degrading performance by disturbing their most sensitive parameters. To address this, we employ a protection strategy that preserves these vulnerable parameters while merging only the less critical ones. Specifically, we identify the top k%percentk\%k % of high-FI parameters in both models and record their locations in ĪAsubscriptĪ _AĪitalic_A and ĪBsubscriptĪ _BĪitalic_B. Then, for each layer in both ĪøTasksubscriptTask _TaskĪøTask and ĪøAvgsubscriptAvg _AvgĪøAvg, we revert parameters at locations in ĪAā©ĪBāsubscriptĪsuperscriptsubscriptĪcomplement _Aā© _B Īitalic_A ā© Īitalic_Bā to their original values from ĪøAsubscript _AĪøitalic_A, and parameters at locations in ĪBā©ĪAāsubscriptĪsuperscriptsubscriptĪcomplement _Bā© _A Īitalic_B ā© Īitalic_Aā to their original values from ĪøBsubscript _BĪøitalic_B. TIES (TrIm, Elect Sign) operates in two steps. First, it reduces redundancy by setting a fraction of the ātask vectorsā Ī“Asubscript _AĪ“italic_A and Ī“Bsubscript _BĪ“italic_B to zero. Then, for each remaining entry, it retains the weight from the vector with the larger absolute value. FI-guided protection can be incorporated into both steps. In the first step, we protect Ī“Asubscript _AĪ“italic_A at locations ĪAsubscriptĪ _AĪitalic_A and Ī“Bsubscript _BĪ“italic_B at ĪBsubscriptĪ _BĪitalic_B from being trimmed. In the second step, entries within ĪAsubscriptĪ _AĪitalic_A are preserved as Ī“Asubscript _AĪ“italic_A, while those in ĪBsubscriptĪ _BĪitalic_B remain as Ī“Bsubscript _BĪ“italic_B, regardless of their absolute values. We merged Qwen2.5-Math-7B (Yang et al., 2024a) and HuatuoGPT-o1-7B (Chen et al., 2024), as both models are further fine-tuned from the same base model, Qwen2.5-7B (Yang et al., 2024b). We evaluate the performance of the merged models on math and health subjects within the MMLU benchmark (Hendrycks et al., 2020). From Table 3, we observe the following: (1) Across all merging methods, FI-guided protection generally enhances the performance of the merged models in both domains. For example, the Average model merging method with FI-guided protection yields approximately a 1% improvement in both the Math and Health domains. (2) Furthermore, TIES with FI protection applied in its first stage performs the best among all merging methods, whereas DARE does not perform well in this setting. FI-protect Math Health Mean Qwen2.5 Math-7B / 0.616 / / Huatuo o1-7B / / 0.724 / Average Without 0.534 (-8.2%) 0.514 (-21.0%) 0.524 With 0.543 (-7.3%) 0.522 (-20.2%) 0.533 Task Without 0.577 (-3.9%) 0.597 (-12.7%) 0.587 With 0.573 (-4.3%) 0.598 (-12.6%) 0.586 TIES Without 0.565 (-5.1%) 0.596 (-12.8%) 0.581 With I 0.583 (-3.3%) 0.606 (-11.8%) 0.595 With I 0.566 (-5.0%) 0.601 (-12.3%) 0.584 DARE Task / 0.573 (-4.3%) 0.589 (-13.5%) 0.581 DARE TIES / 0.560 (-5.6%) 0.588 (-13.6%) 0.574 Table 3: Performance of merging Qwen2.5-Math-7B and HuatuoGPT-o1-7B. The āMeanā column reports the average accuracy across both Math and Health tasks. Blue and cyan percentages indicate the performance drop for the āWithoutā and āWithā variants comparing to the original model, respectively. Figure 5 uses average merging as an example. The results indicate that as the percentage of protected parameters increases, the performance of the merged models initially improves but later declines, highlighting a trade-off in FI-guided protection. Protecting a small proportion of parameters with the highest FI helps mitigate performance degradation caused by parameter conflicts. However, a high percentage of protection may lead to forgetting issues in both domains. To determine the optimal protection percentage, we conduct a hyperparameter search on the validation set. More details can be found in Appendix A.4. Figure 5: Accuracy of average-merged models with FI-guided protection across both domains for different protection percentages k. 6 Conclusion, Limitations & Discussion In summary, we introduced a stability measure, FI, to systematically identify the fragility of LLMs and VLMs (Breach in the shield). Through experiments under both internal and external perturbations, we demonstrate the effectiveness of our proposed method. Our work constitutes an initial attempt to leverage sensitivity measures for improving model performance, focusing primarily on their application to model merging at the inference stage. While our study provides insights into the potential of such measures, we believe that further research is warranted to explore their utility in enhancing model training. 7 Limitations Our method relies on gradient information and is not applicable to āblack-boxā models that do not expose internal parameters or gradients to users. In such cases, text-only approaches like Influence Function (Koh and Liang, 2017) are more suitable. References Ashkboos et al. (2024) Saleh Ashkboos, Maximilian L Croci, Marcelo Gennari do Nascimento, Torsten Hoefler, and James Hensman. 2024. Slicegpt: Compress large language models by deleting rows and columns. arXiv preprint arXiv:2401.15024. Atil et al. (2024) Berk Atil, Alexa Chittams, Liseng Fu, Ferhan Ture, Lixinyu Xu, and Breck Baldwin. 2024. Llm stability: A detailed analysis with some surprises. arXiv preprint arXiv:2408.04667. Bai et al. (2023) Jinze Bai, Shuai Bai, Shusheng Yang, Shijie Wang, Sinan Tan, Peng Wang, Junyang Lin, Chang Zhou, and Jingren Zhou. 2023. Qwen-vl: A frontier large vision-language model with versatile abilities. arXiv preprint arXiv:2308.12966. Brown et al. (2020) Tom Brown, Benjamin Mann, Nick Ryder, Melanie Subbiah, Jared D Kaplan, Prafulla Dhariwal, Arvind Neelakantan, Pranav Shyam, Girish Sastry, Amanda Askell, and 1 others. 2020. Language models are few-shot learners. Advances in neural information processing systems, 33:1877ā1901. Chen et al. (2024) Junying Chen, Zhenyang Cai, Ke Ji, Xidong Wang, Wanlong Liu, Rongsheng Wang, Jianye Hou, and Benyou Wang. 2024. Huatuogpt-o1, towards medical complex reasoning with llms. arXiv preprint arXiv:2412.18925. Cook (1986) R Dennis Cook. 1986. Assessment of local influence. Journal of the Royal Statistical Society Series B: Statistical Methodology, 48(2):133ā155. Dinh et al. (2017) Laurent Dinh, Razvan Pascanu, Samy Bengio, and Yoshua Bengio. 2017. Sharp minima can generalize for deep nets. In International Conference on Machine Learning, pages 1019ā1028. PMLR. Dubey et al. (2024) Abhimanyu Dubey, Abhinav Jauhri, Abhinav Pandey, Abhishek Kadian, Ahmad Al-Dahle, Aiesha Letman, Akhil Mathur, Alan Schelten, Amy Yang, Angela Fan, and 1 others. 2024. The llama 3 herd of models. arXiv preprint arXiv:2407.21783. Dubois et al. (2024) Yann Dubois, BalĆ”zs Galambosi, Percy Liang, and Tatsunori B Hashimoto. 2024. Length-controlled alpacaeval: A simple way to debias automatic evaluators. arXiv preprint arXiv:2404.04475. Fort (2023) Stanislav Fort. 2023. Scaling laws for adversarial attacks on language model activations. arXiv preprint arXiv:2312.02780. Frantar and Alistarh (2023) Elias Frantar and Dan Alistarh. 2023. Sparsegpt: Massive language models can be accurately pruned in one-shot. In International Conference on Machine Learning, pages 10323ā10337. PMLR. Gehman et al. (2020) Samuel Gehman, Suchin Gururangan, Maarten Sap, Yejin Choi, and Noah A Smith. 2020. Realtoxicityprompts: Evaluating neural toxic degeneration in language models. arXiv preprint arXiv:2009.11462. Gong et al. (2024) Xinyu Gong, Jianli Zhang, Qi Gan, Yuxi Teng, Jixin Hou, Yanjun Lyu, Zhengliang Liu, Zihao Wu, Runpeng Dai, Yusong Zou, and 1 others. 2024. Advancing microbial production through artificial intelligence-aided biology. Biotechnology Advances, page 108399. Hendrycks et al. (2020) Dan Hendrycks, Collin Burns, Steven Basart, Andy Zou, Mantas Mazeika, Dawn Song, and Jacob Steinhardt. 2020. Measuring massive multitask language understanding. arXiv preprint arXiv:2009.03300. Ilharco et al. (2022) Gabriel Ilharco, Marco Tulio Ribeiro, Mitchell Wortsman, Suchin Gururangan, Ludwig Schmidt, Hannaneh Hajishirzi, and Ali Farhadi. 2022. Editing models with task arithmetic. arXiv preprint arXiv:2212.04089. Kaushal et al. (2023) Ayush Kaushal, Tejas Vaidhya, and Irina Rish. 2023. Lord: Low rank decomposition of monolingual code llms for one-shot compression. arXiv preprint arXiv:2309.14021. Koh and Liang (2017) Pang Wei Koh and Percy Liang. 2017. Understanding black-box predictions via influence functions. In International conference on machine learning, pages 1885ā1894. PMLR. Lee et al. (2018) Namhoon Lee, Thalaiyasingam Ajanthan, and Philip HS Torr. 2018. Snip: Single-shot network pruning based on connection sensitivity. arXiv preprint arXiv:1810.02340. Li et al. (2024) Junyi Li, Jie Chen, Ruiyang Ren, Xiaoxue Cheng, Wayne Xin Zhao, Jian-Yun Nie, and Ji-Rong Wen. 2024. The dawn after the dark: An empirical study on factuality hallucination in large language models. arXiv preprint arXiv:2401.03205. Liu et al. (2024a) Chris Yuhao Liu, Yaxuan Wang, Jeffrey Flanigan, and Yang Liu. 2024a. Large language model unlearning via embedding-corrupted prompts. arXiv preprint arXiv:2406.07933. Liu et al. (2024b) Haotian Liu, Chunyuan Li, Qingyang Wu, and Yong Jae Lee. 2024b. Visual instruction tuning. Advances in neural information processing systems, 36. (22) Jiashuo Liu, Jiajin Li, Peng Cui, and Jose Blanchet. Stability evaluation of large language models via distributional perturbation analysis. Liu et al. (2018) Qiang Liu, Lihong Li, Ziyang Tang, and Dengyong Zhou. 2018. Breaking the curse of horizon: Infinite-horizon off-policy estimation. Advances in neural information processing systems, 31. Liu et al. (2024c) Yuan Liu, Haodong Duan, Yuanhan Zhang, Bo Li, Songyang Zhang, Wangbo Zhao, Yike Yuan, Jiaqi Wang, Conghui He, Ziwei Liu, and 1 others. 2024c. Mmbench: Is your multi-modal model an all-around player? In European conference on computer vision, pages 216ā233. Springer. Luo et al. (2025) Yingfeng Luo, Tong Zheng, Yongyu Mu, Bei Li, Qinghong Zhang, Yongqi Gao, Ziqiang Xu, Peinan Feng, Xiaoqian Liu, Tong Xiao, and 1 others. 2025. Beyond decoder-only: Large language models can be good encoders for machine translation. arXiv preprint arXiv:2503.06594. Ma et al. (2023) Xinyin Ma, Gongfan Fang, and Xinchao Wang. 2023. Llm-pruner: On the structural pruning of large language models. Advances in neural information processing systems, 36:21702ā21720. Men et al. (2024) Xin Men, Mingyu Xu, Qingyu Zhang, Bingning Wang, Hongyu Lin, Yaojie Lu, Xianpei Han, and Weipeng Chen. 2024. Shortgpt: Layers in large language models are more redundant than you expect. arXiv preprint arXiv:2403.03853. Novak et al. (2018) Roman Novak, Yasaman Bahri, Daniel A Abolafia, Jeffrey Pennington, and Jascha Sohl-Dickstein. 2018. Sensitivity and generalization in neural networks: an empirical study. arXiv preprint arXiv:1802.08760. Ouyang et al. (2025) Shuyin Ouyang, Jie M Zhang, Mark Harman, and Meng Wang. 2025. An empirical study of the non-determinism of chatgpt in code generation. ACM Transactions on Software Engineering and Methodology, 34(2):1ā28. Peng et al. (2024) Sheng Y Peng, Pin-Yu Chen, Matthew Hull, and Duen H Chau. 2024. Navigating the safety landscape: Measuring risks in finetuning large language models. Advances in Neural Information Processing Systems, 37:95692ā95715. Perez et al. (2022) Ethan Perez, Saffron Huang, Francis Song, Trevor Cai, Roman Ring, John Aslanides, Amelia Glaese, Nat McAleese, and Geoffrey Irving. 2022. Red teaming language models with language models. arXiv preprint arXiv:2202.03286. Qi et al. (2024) Xiangyu Qi, Kaixuan Huang, Ashwinee Panda, Peter Henderson, Mengdi Wang, and Prateek Mittal. 2024. Visual adversarial examples jailbreak aligned large language models. In Proceedings of the AAAI Conference on Artificial Intelligence, volume 38, pages 21527ā21536. Shayegani et al. (2023) Erfan Shayegani, Yue Dong, and Nael Abu-Ghazaleh. 2023. Jailbreak in pieces: Compositional adversarial attacks on multi-modal language models. arXiv preprint arXiv:2307.14539. Simonyan et al. (2013) Karen Simonyan, Andrea Vedaldi, and Andrew Zisserman. 2013. Deep inside convolutional networks: Visualising image classification models and saliency maps. arXiv preprint arXiv:1312.6034. Sun et al. (2023) Mingjie Sun, Zhuang Liu, Anna Bair, and J Zico Kolter. 2023. A simple and effective pruning approach for large language models. arXiv preprint arXiv:2306.11695. Touvron et al. (2023) Hugo Touvron, Thibaut Lavril, Gautier Izacard, Xavier Martinet, Marie-Anne Lachaux, TimothĆ©e Lacroix, Baptiste RoziĆØre, Naman Goyal, Eric Hambro, Faisal Azhar, and 1 others. 2023. Llama: Open and efficient foundation language models. arXiv preprint arXiv:2302.13971. Wei et al. (2024) Boyi Wei, Kaixuan Huang, Yangsibo Huang, Tinghao Xie, Xiangyu Qi, Mengzhou Xia, Prateek Mittal, Mengdi Wang, and Peter Henderson. 2024. Assessing the brittleness of safety alignment via pruning and low-rank modifications. In Proceedings of the 41st International Conference on Machine Learning, pages 52588ā52610. Wortsman et al. (2022) Mitchell Wortsman, Gabriel Ilharco, Samir Ya Gadre, Rebecca Roelofs, Raphael Gontijo-Lopes, Ari S Morcos, Hongseok Namkoong, Ali Farhadi, Yair Carmon, Simon Kornblith, and 1 others. 2022. Model soups: averaging weights of multiple fine-tuned models improves accuracy without increasing inference time. In International conference on machine learning, pages 23965ā23998. PMLR. Xie and Lukasiewicz (2023) Zhongbin Xie and Thomas Lukasiewicz. 2023. An empirical analysis of parameter-efficient methods for debiasing pre-trained language models. arXiv preprint arXiv:2306.04067. Yadav et al. (2024) Prateek Yadav, Derek Tam, Leshem Choshen, Colin A Raffel, and Mohit Bansal. 2024. Ties-merging: Resolving interference when merging models. Advances in Neural Information Processing Systems, 36. Yang et al. (2024a) An Yang, Baosong Yang, Beichen Zhang, Binyuan Hui, Bo Zheng, Bowen Yu, Chengyuan Li, Dayiheng Liu, Fei Huang, Haoran Wei, and 1 others. 2024a. Qwen2. 5 technical report. arXiv preprint arXiv:2412.15115. Yang et al. (2024b) An Yang, Beichen Zhang, Binyuan Hui, Bofei Gao, Bowen Yu, Chengpeng Li, Dayiheng Liu, Jianhong Tu, Jingren Zhou, Junyang Lin, and 1 others. 2024b. Qwen2. 5-math technical report: Toward mathematical expert model via self-improvement. arXiv preprint arXiv:2409.12122. Yang et al. (2024c) Enneng Yang, Li Shen, Guibing Guo, Xingwei Wang, Xiaochun Cao, Jie Zhang, and Dacheng Tao. 2024c. Model merging in llms, mllms, and beyond: Methods, theories, applications and opportunities. arXiv preprint arXiv:2408.07666. Yu et al. (2024) Le Yu, Bowen Yu, Haiyang Yu, Fei Huang, and Yongbin Li. 2024. Language models are super mario: Absorbing abilities from homologous models as a free lunch. In Forty-first International Conference on Machine Learning. Zhang et al. (2024) Xinyu Zhang, Huiyu Xu, Zhongjie Ba, Zhibo Wang, Yuan Hong, Jian Liu, Zhan Qin, and Kui Ren. 2024. Privacyasst: Safeguarding user privacy in tool-using large language model agents. IEEE Transactions on Dependable and Secure Computing. Zheng et al. (2025a) Tong Zheng, Lichang Chen, Simeng Han, R Thomas McCoy, and Heng Huang. 2025a. Learning to reason via mixture-of-thought for logical reasoning. arXiv preprint arXiv:2505.15817. Zheng et al. (2025b) Tong Zheng, Yan Wen, Huiwen Bao, Junfeng Guo, and Heng Huang. 2025b. Asymmetric conflict and synergy in post-training for llm-based multilingual machine translation. arXiv preprint arXiv:2502.11223. Zhou et al. (2022) Kaiyang Zhou, Jingkang Yang, Chen Change Loy, and Ziwei Liu. 2022. Learning to prompt for vision-language models. International Journal of Computer Vision, 130(9):2337ā2348. Zhu et al. (2007) Hongtu Zhu, Joseph G Ibrahim, Sikyum Lee, and Heping Zhang. 2007. Perturbation selection and influence measures in local influence analysis. 35:2565ā2588. Zhu et al. (2011) Hongtu Zhu, Joseph G Ibrahim, and Niansheng Tang. 2011. Bayesian influence analysis: a geometric approach. Biometrika, 98(2):307ā323. Appendix A Appendices A.1 Detail of Parameter sparsification experiment Experiment on MMLU We conduct experiments on the multiple-choice problems from the MMLU (Hendrycks et al., 2020) dataset, using Qwen2-7B. We take the cross-entropy loss, i.e., f=ālogā”Pā¢(y=ypred|x,Īø)conditionalsubscriptpredf=- P(y=y_pred|x,Īø)f = - log P ( y = ypred | x , Īø ), as the target function, and calculate the FI value according to Theorem 3.3. In this setup, we treat the task as a 4-class classification problem with the possible classes being "A," "B," "C," and "D". Experiment on Alpaca-Eval We use the Alpaca-eval validation set (Dubois et al., 2024), a widely adopted benchmark, and conduct experiments with various open-source models, including LLaMA2, LLaMA3 (Touvron et al., 2023), and Qwen2 (Bai et al., 2023), across different sizes. We report two metrics: ROUGE-1 (comparing to pre-sparsity responses) and length-control winning rate (LCWR), comparing to GPT-3.5 Turbo. Higher scores are better for both metrics. To estimate the average FI for sequence generation, we use the fixed-context approach with L=55L=5L = 5. For each sample z, we estimate lā¢(z)subscriptFI_l(z)FIitalic_l ( z ) by generating N=1010N=10N = 10 responses, truncating them at position lā11l-1l - 1. These truncated sequences are used to approximate the conditional expectation by computing the sample average. The per-token FI values are then aggregated using Equation 3 to obtain seqLā¢(z)subscriptsuperscriptseqFI^L_seq(z)FIitalic_Lseq ( z ), which is averaged across all samples to estimate the overall FI. A.2 Computation complexity analysis Let n denote the number of samples, p denote the dimension of perturbation (p=33p=3p = 3 for pixel-wise computations and p=11p=1p = 1 for parameter-wise computations), and d represent the total number of pixels or parameters. Computational Complexity Analysis: ⢠Jacobian-Norm: ā¢(nā¢pā¢d)O(npd)O ( n p d ), arising from gradient computation per pixel/parameter. ⢠Saliency-Map: Identical to Jacobian-Norm, ā¢(nā¢pā¢d)O(npd)O ( n p d ). ⢠FI-inverse: ā¢(nā¢p3ā¢d+nā¢pā¢d)superscript3O(np^3d+npd)O ( n p3 d + n p d ), with ā¢(p3)superscript3O(p^3)O ( p3 ) from inverse matrix computations and ā¢(nā¢pā¢d)O(npd)O ( n p d ) from gradient calculations. ⢠FI-cSVD (our method): ā¢(nā¢p2ā¢r0ā¢d+nā¢pā¢d)superscript2subscript0O(np^2r_0d+npd)O ( n p2 r0 d + n p d ), where ā¢(p2ā¢r0)superscript2subscript0O(p^2r_0)O ( p2 r0 ) stems from the compact SVD used to compute matrix inversion efficiently. In practical scenarios: 1. Parameter-wise stability: Since individual parameters have dimension p=11p=1p = 1, the FI calculation reduces to scalar inversion, thus the complexity simplifies to ā¢(nā¢pā¢d)O(npd)O ( n p d ), matching Jacobian-Norm and Saliency-Map. 2. Pixel-wise stability (image data): Given that each pixel has dimension p=33p=3p = 3 (RGB), the FI calculation involves compact SVD for a 3Ć3333Ć 33 Ć 3 matrix. Theoretically, this makes our method about 9999 times slower compared to baseline methods. However, in practical implementation, our approach is only approximately 2222 times slower. The table below presents the average time required to compute FI, Saliency Map, and Jacobian Norm for a single image using Qwen2VL-7B. All results are averaged over 100 images and measured on an A100-80G GPU. Method Time (s) FI 0.3828 Saliency-Map 0.1964 Jacobian-Norm 0.1939 Table 4: Empirical computation times for different methods. A.3 Reparametrization Invariance of FI The proposed FI measure has the property of transformation invariance. Theorem A.1 (Reparametrization invariance). Suppose that Ļitalic-ĻĻĻ is a diffeomorphism of Ļ. Then, Fā¢IĻā¢(Ļ0)subscriptsubscript0FI_Ļ( _0)F Iitalic_Ļ ( Ļ0 ) is invariant with respect to any reparameterization corresponding to Ļitalic-ĻĻĻ. Specifically, let Ļ~ā¢(t)=ĻāĻā¢(t),Ļ~0=Ļā¢(Ļ0),formulae-sequence~italic-Ļsubscript~0italic-Ļsubscript0 Ļ(t)=Ļ Ļ(t), Ļ_0=Ļ( _0),over~ start_ARG Ļ end_ARG ( t ) = Ļ ā Ļ ( t ) , over~ start_ARG Ļ end_ARG0 = Ļ ( Ļ0 ) , we have Fā¢IĻ~ā¢(Ļ~0)=Fā¢IĻā¢(Ļ0).subscript~subscript~0subscriptsubscript0FI_ Ļ( Ļ_0)=FI_Ļ( _0).F Iover~ start_ARG Ļ end_ARG ( over~ start_ARG Ļ end_ARG0 ) = F Iitalic_Ļ ( Ļ0 ) . The detailed proof can be found in shu2019sensitivity. Theorem A.1 establishes that Fā¢IĻā¢(Ļ0)subscriptsubscript0FI_Ļ( _0)F Iitalic_Ļ ( Ļ0 ) is invariant under any diffeomorphic (e.g., scaling and spinning) reparameterization of the original perturbation. This invariance property is not shared by other measures, such as Jacobian norm (Novak et al., 2018), Cookās local influence measure (Cook, 1986), and Sharpness (Novak et al., 2018). For instance, consider a perturbation of the form α+Īā¢Ī±Īα+ α + Πα, where α is a subvector of (xā¤,Īøā¤)ā¤superscriptsuperscripttopsuperscripttoptop(x ,Īø ) ( x⤠, θ⤠)ā¤. If we apply a scaling reparameterization αā²=Kāαsuperscriptā²direct-productα =K αⲠ= K ā α, where K is a scaling vector and ādirect-product ā denotes element-wise multiplication, then the Jacobian norms change: āJā¢(α)āF=[āi(āfāαi)2]1/2ā āJā¢(αā²)āF.subscriptnormsuperscriptdelimited-[]subscriptsuperscriptsubscript212subscriptnormsuperscriptā²\|J(α)\|_F= [ _i ( ā fā _i% )^2 ]^1/2ā \|J(α )\|_F.ā„ J ( α ) ā„F = [ āi ( divide start_ARG ā f end_ARG start_ARG ā αitalic_i end_ARG )2 ]1 / 2 ā ā„ J ( αⲠ) ā„F . In contrast, the FI measure remains unchanged. Such a reparameterization does not alter the function itself but may affect the measure values, potentially weakening the correlation between perturbation and performance degradation. A similar discussion can be found in (Dinh et al., 2017). A.4 Detail of FI-guided protection in model merging Table 5: Searched ranges of hyperparameters of model merging methods Hyper parameter Search Ranges of Hyperparameters Protecting ratio k [1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%] Weight parameter γ in Task Arithmetic & TIES [0.3, 0.4, 0.5, 0.6, 0.9, 1.0] A.5 Proof of Theorem 3.3 Proof. We apply Taylor expansion to fā¢(Ļā¢(t))f(Ļ(t))f ( Ļ ( t ) ) at the point Ļā¢(t)Ļ(t)Ļ ( t ): fā¢(Ļā¢(t)) f(Ļ(t))f ( Ļ ( t ) ) =fā¢(Ļā¢(0))+āfā¢(Ļ0)ThĻ0ā¢t+12ā¢(hĻ0Tā¢Hfā¢(Ļ0)ā¢hĻ0+āfā¢(Ļ0)Td2ā¢Ļā¢(0)/dā¢t2)ā¢t2+oā¢(t2),absent0superscriptsubscriptāsubscript0subscriptāsubscript012superscriptsubscriptāsubscript0subscriptsubscript0subscriptāsubscript0superscriptsubscriptāsubscript0superscript20superscript2superscript2superscript2 =f(Ļ(0))+ā_f (Ļ_0 )^Th_% Ļ_0t+ 12 (h_Ļ_0^TH_f (Ļ_0% )h_Ļ_0+ā_f (Ļ_0 )^Td^2% Ļ(0)/dt^2 )t^2+o (t^2 ),= f ( Ļ ( 0 ) ) + āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT hitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT t + divide start_ARG 1 end_ARG start_ARG 2 end_ARG ( hitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPTT Hitalic_f ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT hitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT d2 Ļ ( 0 ) / d t2 ) t2 + o ( t2 ) , where āfā¢(Ļ0)=āfā¢(Ļ)/āĻ|Ļ=Ļ0subscriptāsubscript0evaluated-atsubscript0ā_f (Ļ_0 )=ā f(Ļ)/ .ā% Ļ |_Ļ=Ļ_0āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT = ā f ( Ļ ) / ā Ļ |Ļ = Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and Hfā¢(Ļ0)=ā2fā¢(Ļ)/āĻā¢āĻT|Ļ=Ļ0subscriptsubscript0evaluated-atsuperscript2superscriptsubscript0H_f (Ļ_0 )=ā^2f(Ļ)/ .āĻ% āĻ^T |_Ļ=Ļ_0Hitalic_f ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT = ā2 f ( Ļ ) / ā Ļ ā Ļitalic_T |Ļ = Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT. From the definition of SCsubscriptS_CSitalic_C, SC2ā¢(Ļt,Ļ0)superscriptsubscript2subscriptsubscript0S_C^2( _t, _0)Sitalic_C2 ( Ļitalic_t , Ļ0 ) can be approximated as SC2ā¢(Ļt,Ļ0)=t2ā¢hĻ0Tā¢GĻ0ā¢hĻ0+oā¢(t2).superscriptsubscript2subscriptsubscript0superscript2superscriptsubscriptāsubscript0subscriptsubscript0subscriptāsubscript0superscript2S_C^2(Ļ_t,Ļ_0)=t^2h_Ļ_0^TG_ _0% h_Ļ_0+o (t^2 ).Sitalic_C2 ( Ļitalic_t , Ļ0 ) = t2 hitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPTT Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT hitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT + o ( t2 ) . Based on lāHĖopitalās rule, the stability measure FI from Equation1 can be rewritten as: Ļā¢(Ļ0)=maxhĻā”hĻTā¢āfā¢(Ļ0)āfā¢(Ļ0)Tā”hĻhĻTā¢GĻ0ā¢hĻ.subscriptsubscript0subscriptsubscriptāsuperscriptsubscriptāsubscriptāsubscript0superscriptsubscriptāsubscript0subscriptāsuperscriptsubscriptāsubscriptsubscript0subscriptāFI_Ļ (Ļ_0 )= _h_Ļ h% _Ļ^Tā_f (Ļ_0 )ā_f (Ļ_% 0 )^Th_Ļh_Ļ^TG_ _0h_Ļ.FIitalic_Ļ ( Ļ0 ) = maxitalic_h start_POSTSUBSCRIPT Ļ end_POSTSUBSCRIPT divide start_ARG hitalic_Ļitalic_T āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT hitalic_Ļ end_ARG start_ARG hitalic_Ļitalic_T Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT hitalic_Ļ end_ARG . We then reparameterize Ļ to Ļ~=GĻ0ā1/2ā¢Ļ~superscriptsubscriptsubscript012 Ļ=G_ _0^-1/2 ~ start_ARG Ļ end_ARG = Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 Ļ. According to Theorem A.1, the stability measure FIFI remains invariant under this reparameterization Fā¢IĻā¢(Ļ0)=Fā¢IĻ~ā¢(Ļ~0)=argā¢maxhĻ~ā”hĻ~ā¤ā¢GĻ0ā1/2ā¢āfā¢(Ļ0)āfā¢(Ļ0)ā¤ā”GĻ0ā1/2ā¢hĻ~hĻ~ā¤ā¢hĻ~.subscriptsubscript0subscript~subscript~0subscriptargmaxsubscriptā~superscriptsubscriptā~topsubscriptsuperscript12subscript0subscriptāsubscript0subscriptsuperscriptātopsubscript0subscriptsuperscript12subscript0subscriptā~superscriptsubscriptā~topsubscriptā~FI_Ļ( _0)=FI_ Ļ( Ļ_0)= *% arg\,max_h_ Ļ h_ Ļ G^-1/2_Ļ% _0 _f( _0)ā _f( _0)G^-1/2_ _0h% _ Ļh_ Ļ h_ Ļ.F Iitalic_Ļ ( Ļ0 ) = F Iover~ start_ARG Ļ end_ARG ( over~ start_ARG Ļ end_ARG0 ) = start_OPERATOR arg max end_OPERATORh start_POSTSUBSCRIPT over~ start_ARG Ļ end_ARG end_POSTSUBSCRIPT divide start_ARG hover~ start_ARG Ļ end_ARG⤠G- 1 / 2Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT āā¤f ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT G- 1 / 2Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT hover~ start_ARG Ļ end_ARG end_ARG start_ARG hover~ start_ARG Ļ end_ARG⤠hover~ start_ARG Ļ end_ARG end_ARG . The maximization problem is now in the form of a Rayleigh quotient, which attains its maximum when hĻ~subscriptā~h_ Ļhover~ start_ARG Ļ end_ARG is proportional to GĻ0ā1/2ā¢āfā¢(Ļ0)superscriptsubscriptsubscript012subscriptāsubscript0G_ _0^-1/2 _f( _0)Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT. Substituting back into the Rayleigh quotient, we find: Ļā¢(Ļ0)subscriptsubscript0 _Ļ( _0)FIitalic_Ļ ( Ļ0 ) =(GĻ0ā1/2ā¢āfā¢(Ļ0))Tā¢GĻ0ā1/2ā¢āfā¢(Ļ0)āfā¢(Ļ0)Tā”GĻ0ā1/2ā¢(GĻ0ā1/2ā¢āfā¢(Ļ0))(GĻ0ā1/2ā¢āfā¢(Ļ0))Tā¢(GĻ0ā1/2ā¢āfā¢(Ļ0))absentsuperscriptsuperscriptsubscriptsubscript012subscriptāsubscript0superscriptsubscriptsubscript012subscriptāsubscript0superscriptsubscriptāsubscript0superscriptsubscriptsubscript012superscriptsubscriptsubscript012subscriptāsubscript0superscriptsuperscriptsubscriptsubscript012subscriptāsubscript0superscriptsubscriptsubscript012subscriptāsubscript0 = (G_ _0^-1/2 _f( _0) )^% TG_ _0^-1/2 _f( _0) _f( _0)^TG_% _0^-1/2 (G_ _0^-1/2 _f( _0) )% (G_ _0^-1/2 _f( _0) )^T (G_ _0% ^-1/2 _f( _0) )= divide start_ARG ( Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT )T Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 ( Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT ) end_ARG start_ARG ( Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT )T ( Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 / 2 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT ) end_ARG =āfā¢(Ļ0)TGĻ0ā1ā¢āfā¢(Ļ0)āfā¢(Ļ0)Tā”GĻ0ā1ā¢āfā¢(Ļ0)āfā¢(Ļ0)TGĻ0ā1ā¢āfā¢(Ļ0)absentsuperscriptsubscriptāsubscript0superscriptsubscriptsubscript01subscriptāsubscript0superscriptsubscriptāsubscript0superscriptsubscriptsubscript01subscriptāsubscript0superscriptsubscriptāsubscript0superscriptsubscriptsubscript01subscriptāsubscript0 = _f( _0)^TG_ _0^-1 _f(% _0) _f( _0)^TG_ _0^-1 _f( _0)% _f( _0)^TG_ _0^-1 _f( _0)= divide start_ARG āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT end_ARG start_ARG āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT end_ARG =āfā¢(Ļ0)TGĻ0ā1ā¢āfā¢(Ļ0).absentsuperscriptsubscriptāsubscript0superscriptsubscriptsubscript01subscriptāsubscript0 = _f( _0)^TG_ _0^-1 _f( _0% ).= āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPTT Gitalic_Ļ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT- 1 āf ( Ļ start_POSTSUBSCRIPT 0 ) end_POSTSUBSCRIPT . This concludes the proof. ā A.6 Additional experiment results on parameter sparsification Table 6: Performance of Different Models Based on Criteria with Full Value and Sparsity Percentages Model Criteria Full 6% Sparsity 8% Sparsity 10% Sparsity 12% Sparsity FI-High Random FI-High Random FI-High Random FI-High Random Llama2-13B Rouge-1 1.0 0.52 0.59±0.02plus-or-minus0.590.020.59± 0.020.59 ± 0.02 0.4 0.43±0.06plus-or-minus0.430.060.43± 0.060.43 ± 0.06 0.18 0.68±0.01plus-or-minus0.680.010.68± 0.010.68 ± 0.01 0.05 0.19±0.03plus-or-minus0.190.030.19± 0.030.19 ± 0.03 LCWR 0.43 0.38 0.41±0.03plus-or-minus0.410.030.41± 0.030.41 ± 0.03 0.29 0.34±0.07plus-or-minus0.340.070.34± 0.070.34 ± 0.07 0.09 0.42±0.0plus-or-minus0.420.00.42± 0.00.42 ± 0.0 0.01 0.08±0.05plus-or-minus0.080.050.08± 0.050.08 ± 0.05 Llama3-8B Rouge-1 1.0 0.46 0.52±0.04plus-or-minus0.520.040.52± 0.040.52 ± 0.04 0.21 0.41±0.06plus-or-minus0.410.060.41± 0.060.41 ± 0.06 0.09 0.25±0.04plus-or-minus0.250.040.25± 0.040.25 ± 0.04 0.04 0.12±0.03plus-or-minus0.120.030.12± 0.030.12 ± 0.03 LCWR 0.42 0.4 0.38±0.01plus-or-minus0.380.010.38± 0.010.38 ± 0.01 0.12 0.30±0.03plus-or-minus0.300.030.30± 0.030.30 ± 0.03 0.0 0.12±0.01plus-or-minus0.120.010.12± 0.010.12 ± 0.01 0.0 0.01±0.01plus-or-minus0.010.010.01± 0.010.01 ± 0.01 Llama2-7B Rouge-1 1.0 0.44 0.56±0.01plus-or-minus0.560.010.56± 0.010.56 ± 0.01 0.25 0.45±0.02plus-or-minus0.450.020.45± 0.020.45 ± 0.02 0.06 0.33±0.02plus-or-minus0.330.020.33± 0.020.33 ± 0.02 0.0 0.21±0.02plus-or-minus0.210.020.21± 0.020.21 ± 0.02 LCWR 0.42 0.32 0.4±0.0plus-or-minus0.40.00.4± 0.00.4 ± 0.0 0.12 0.35±0.01plus-or-minus0.350.010.35± 0.010.35 ± 0.01 0.0 0.19±0.05plus-or-minus0.190.050.19± 0.050.19 ± 0.05 0.0 0.1±0.03plus-or-minus0.10.030.1± 0.030.1 ± 0.03 Qwen2-7B Rouge-1 1.0 0.09 0.41±0.05plus-or-minus0.410.050.41± 0.050.41 ± 0.05 0.01 0.30±0.09plus-or-minus0.300.090.30± 0.090.30 ± 0.09 0.01 0.31±0.06plus-or-minus0.310.060.31± 0.060.31 ± 0.06 0.01 0.15±0.02plus-or-minus0.150.020.15± 0.020.15 ± 0.02 LCWR 0.41 0.03 0.35±0.03plus-or-minus0.350.030.35± 0.030.35 ± 0.03 0.02 0.25±0.1plus-or-minus0.250.10.25± 0.10.25 ± 0.1 0.03 0.20±0.05plus-or-minus0.200.050.20± 0.050.20 ± 0.05 0.03 0.08±0.02plus-or-minus0.080.020.08± 0.020.08 ± 0.02 Qwen2-1.5B Rouge-1 1.0 0.18 0.4±0.13plus-or-minus0.40.130.4± 0.130.4 ± 0.13 0.16 0.32±0.02plus-or-minus0.320.020.32± 0.020.32 ± 0.02 0.05 0.28±0.08plus-or-minus0.280.080.28± 0.080.28 ± 0.08 0.05 0.23±0.02plus-or-minus0.230.020.23± 0.020.23 ± 0.02 LCWR 0.14 0.03 0.07±0.04plus-or-minus0.070.040.07± 0.040.07 ± 0.04 0.04 0.02±0.02plus-or-minus0.020.020.02± 0.020.02 ± 0.02 0.0 0.04±0.0plus-or-minus0.040.00.04± 0.00.04 ± 0.0 0.0 0.02±0.02plus-or-minus0.020.020.02± 0.020.02 ± 0.02