Paper deep dive
Scalability of Morality: A Particle-Based Numerical Study on the Decoupling of Law and Ethics in Large-Scale Populations
Amir Arslan Haghrah, Amir Aslan Haghrah
Intelligence
Status: succeeded | Model: Gemma-4-26B-A4B | Prompt: intel-v1 | Confidence: 93%
Last extracted: 7/9/2026, 8:21:33 AM
Summary
This study introduces a particle-based computational framework to investigate how morality scales in expanding populations. It demonstrates that when population size exceeds individual memory capacity, local re-encounter probabilities drop, neutralizing decentralized peer-to-peer accountability. This structural dilution triggers a non-linear phase transition where global behavioral norms decouple from ethical baselines and drift toward a minimalist legal floor, characterized by path-dependent hysteresis and irreversible moral decay.
Entities (8)
Relation Signals (5)
Monte Carlo ensemble simulations → demonstrates → Non-linear phase transition
confidence 96% · Monte Carlo ensemble simulations demonstrate a distinct, non-linear phase transition as the population scales (N→∞).
Particle-based computational framework → models → Decoupling of law and ethics
confidence 95% · This study introduces a particle-based computational framework to investigate the scalability of morality and the systemic decoupling of formal law from decentralized social ethics in expanding populations.
Population size (N) → outpaces → Finite memory capacity (L)
confidence 94% · When the population metric outpaces memory capacity (N ≫ L), the local re-encounter probability drops as O(L/N).
Cyclic scale experiments → exposes → Path-dependent hysteresis loop
confidence 93% · Furthermore, cyclic scale experiments expose a prominent, path-dependent hysteresis loop, mathematically formalizing the non-Markovian inertia and irreversible nature of moral decay in self-organizing social systems.
Drop in local re-encounter probability → neutralizes → Decentralized peer-to-peer accountability
confidence 92% · This structural dilution neutralizes decentralized peer-to-peer accountability, causing global behavioral norms to decouple from moral baselines and drift toward a minimalist legal floor.
Cypher Suggestions (0)
No Cypher suggestions yet.
Abstract
Abstract:This study introduces a particle-based computational framework to investigate the scalability of morality and the systemic decoupling of formal law from decentralized social ethics in expanding populations. While micro-societies reinforce ethical conduct through local reciprocity, macroscale systems introduce anonymity that strains cognitive memory limitations. We model individual agents as discrete particles with finite memory capacities ($L$) and dynamically evolving, stochastic choice profiles ($\mu$) regulated by non-linear social pressure switches. Monte Carlo ensemble simulations demonstrate a distinct, non-linear phase transition as the population scales ($N \to \infty$). When the population metric outpaces memory capacity ($N \gg L$), the local re-encounter probability drops as $\mathcal{O}(L/N)$. This structural dilution neutralizes decentralized peer-to-peer accountability, causing global behavioral norms to decouple from moral baselines and drift toward a minimalist legal floor. Furthermore, cyclic scale experiments expose a prominent, path-dependent hysteresis loop, mathematically formalizing the non-Markovian inertia and irreversible nature of moral decay in self-organizing social systems.
Tags
Links
- Source: https://arxiv.org/abs/2606.27039v1
- Canonical: https://arxiv.org/abs/2606.27039v1
Trouble viewing inline? Open PDF directly →
Full Text
50,089 characters extracted from source content.
Expand or collapse full text
Scalability of Morality: A Particle-Based Numerical Study on the Decoupling of Law and Ethics in Large-Scale Populations Amir Arslan Haghrah, Amir Aslan Haghrah Abstract This study introduces a particle-based computational framework to investigate the scalability of morality and the systemic decoupling of formal law from decentralized social ethics in expanding populations. While micro-societies reinforce ethical conduct through local reciprocity, macroscale systems introduce anonymity that strains cognitive memory limitations. We model individual agents as discrete particles with finite memory capacities (L) and dynamically evolving, stochastic choice profiles (μ) regulated by non-linear social pressure switches. Monte Carlo ensemble simulations demonstrate a distinct, non-linear phase transition as the population scales (N→∞N→∞). When the population metric outpaces memory capacity (N≫LN L), the local re-encounter probability drops as (L/N)O(L/N). This structural dilution neutralizes decentralized peer-to-peer accountability, causing global behavioral norms to decouple from moral baselines and drift toward a minimalist legal floor. Furthermore, cyclic scale experiments expose a prominent, path-dependent hysteresis loop, mathematically formalizing the non-Markovian inertia and irreversible nature of moral decay in self-organizing social systems. 1 Introduction The dual nature of social order is anchored in the distinction between legality and morality [16]. While ethics represents the internal normative framework that guides individual conduct toward collective well-being, legal systems serve as the external, codified apparatus designed to enforce social cooperation through a structured regime of rewards and punishments [13, 4]. Historically, the stability of the human civilisation has relied on the alignment of these two pillars. However, a critical divergence exists: legal frameworks are often ”minimalist” by design, ensuring only the prevention of gross violations while leaving a vast territory of ”unethical but legal” behaviors open for exploitation [27, 5]. In such spaces, individuals may derive significant utility from actions that, while technically permissible under the law, fundamentally erode the ethical fabric of the community [25]. The efficacy of moral regulation is inherently sensitive to the architecture of the society in which it operates. In small-scale, high density social configurations, ethical behavior is self-reinforcing through a decentralized system of individual sanctions and rewards [22]. In these ”micro-societies,” the high probability of recurrent interactions creates a robust feedback loop of accountability [1]; the social cost of an unethical act, manifested through reputation loss or direct retaliation, often outweights the immediate gain [21]. However, as societies scale toward massive, complex systems, a fundamental phase transition occurs [9]. This study addresses the ”Scalability of Morality” by examining how population growth decouples formal law from ethical behavior. As the number of interacting particles increases, the anonymity afforded by system scale significantly dilutes individual-level sanctions [28]. In large-scale systems, stochastic interactions ensure an unethical actor is unlikely to re-encounter the same victim, neutralizing individual retribution networks [23]. Consequently, the burden of maintaining social integrity shifts entirely to the legal framework, which is structurally ill-equipped to govern nuanced ethical conduct. Utilizing a numerical framework grounded in particle dynamics and Monte Carlo simulations, we quantify this decoupling as a distinct nonequilibrium phase transition [26]. Identifying the mathematical threshold at which moral accountability collapses is a vital engineering necessity for designing resilient governance structures in an era of unprecedented global connectivity and urban density. Following this conceptual framework, the remainder of this introduction provides a concise review of existing literature to contextualize our contribution, followed by a summary of the paper’s novelty and structural organization. The quantitative study of the social behavior has its roots in Sociophysics, where the interactions of individuals are modeled using the principles of statistical mechanics [15]. Early models, such as the Ising model or Schelling’s segregation model, demonstrated that macroscopic social patterns emerge from microscopic rules of interaction [6]. Recent advancements in Computational Social Science [18] have expanded these concepts to study the evolution of cooperation under varying structural and systemic complexities [12]. Building upon this paradigm, recent literature in Computational Social Science has shifted toward decoding how systemic architecture and communication limits modify cooperative equilibria. Large-scale global data sets and advanced data processing tools have allowed researchers to document a profound tension between growing institutional complexity and localized social cohesion [7]. While macro-level structures frequently rely on generalized institutional rules to sustain large populations, behavioral experiments integrated with agent-based models reveal that decentralized monitoring functions efficiently only when individual interactions remain visible and localized [4]. When these interactions are stretched across expansive, sparse digital or physical networks, the cognitive and structural limits of informal social control become a primary bottleneck for collective stability [10]. To precisely isolate these dynamics, modern studies frequently turn to specialized multi-agent modeling frameworks. In these environments, autonomous decision-making entities are governed by heterogeneous behavioral rule specifications, allowing for the precise analysis of emergent phenomena across varying spatio-temporal scales [14]. Recent simulations mapping the evolution of pre-state communities indicate that while decentralized information exchange and gossip-driven reputation assessment can suppress opportunism in small configurations, long-term systemic scalability invariably triggers a structural reliance on centralized, coercive frameworks [29]. Consequently, computational modelers have increasingly positioned agent-based frameworks as vital simulation sandboxes to test the bounds of institutional policy and public governance prior to structural deployment [3]. This collapse of localized moral accountability under rapid scaling can be formally mapped to the physics of phase transitions. Within nonequilibrium statistical mechanics, collective human behavior is modeled as an active matter system where microscopic adjustments in density or proximity induce spontaneous macroscopic symmetry breaking [26]. Simulating spatial evolutionary games via Monte Carlo updates reveals that shifting from localized lattices to massive, well-mixed networks yields sharp bifurcations between universal cooperation and systemic defection [20]. Furthermore, multi-population mean-field theories confirm that when agent density thresholds cross a critical value, the global order parameter experiences an abrupt jump [8]. Our paper contributes directly to this intersection by using particle dynamics to mathematically bound the operational threshold at which decentralized moral regulation structurally fails. Previous research in the field of Sociophysics has successfully used Monte Carlo methods to model phenomena like opinion formation and cultural dissemination. These studies demonstrate that the ”global state” of a system is highly sensitive to the local ”interaction density” [15]. Furthermore, the literature on Agent-Based Modeling (ABM) has explored how local feedback loops, such as individual rewards and punishments, maintain system cooperative stability [11]. However, a significant gap remains: most models assume a uniform enforcement of rules, failing to account for the divergent dynamics between centeralized legal constraints and stochastic individual interactions [17]. While the ”Anonymity Effect” in large-scale systems has been qualitatively discussed in social sciences, it lacks a rigorous numerical characterization regarding how group size actively alters individual behavioral strategic payoffs [24]. Specifically, how the scaling of a system (N⟶∞N ∞) physically ”dilutes” the impact of local individual sanctions, leading to a decoupling from global legal frameworks, is a problem of systemic stability that has yet to be fully quantified through the lens of particle-based simulations and nonequilibrium spin models [23, 2]. The evolution of complex Cyber-Physical-Social Systems (CPSS) underscores the difficulty of modeling human behavior under deep systemic uncertainty [30]. Because data-centric learning models often struggle to integrate domain expertise with behavioral ambiguities, contemporary literature emphasizes specialized mathematical frameworks—such as fuzzy logic and decentralized architectures—to bridge the gap between microscopic agent behaviors and macroscopic social trends [30]. Moving away from black-box predictive models, this study introduces a particle-based stochastic framework that maps non-linear social pressures directly onto state-space dynamics. Our approach builds upon the established paradigm of modeling competitive socio-behavioral dynamics as non-linear, parameter-driven systems rather than relying on empirical big data. As demonstrated by Li et al. [19], formalizing social phenomena like opinion propagation as evolving network parameters offers deep theoretical insights into stability, convergence, and attractor dominance. Extending this methodology to normative dynamics, we treat individual ethical propensities (μ) and motivational velocities (ρ) as coupled dynamic parameters to evaluate how macroscale structural changes drive a self-organizing population toward competing, path-dependent equilibria. This article departs from qualitative sociological narratives and static optimization models, offering a dynamic, numerical prespective on social ethics. The primary innovations are: • Formalization of Decoupling Dynamics: We define a mathematical framework to observe how the correlation between ethical behavior and legal compliance weakens as a function of system entropy and population scale. • Stochastic Interaction Modeling: Using a particle-based approach, we simulate the probability of encounter-based sanctions, providing a more realistic representation of social accountability than traditional mean-field approximations. • Identification of the Critical Scaling Point: We identify the numerical threshold where the decentralized ”ethical feedback” fails, leaving the system’s integrity solely dependent on centralized (legal) control. The remainder of this paper is organized as follows: In Section 2, the system architecture and modeling are presented, detailing the mathematical formulation of the particle system, the transition rules for rewards and punishments, and the Monte Carlo setup. In Section 3, we provide the simulation results and sensitivity analysis, which analyze the behavior of the model under varying population densities and examine the decoupling phenomenon. Section 4 offers a systemic discussion, interpreting the findings within the context of control theory and social stability. Finally, Section 5 summarizes the results and outlines future directions for large-scale social modeling. 2 System Architecture and Modeling The proposed model simulates a society as a complex system of interacting particles, where each particle represents an individual agent capable of making stochastic decisions regarding ethical conduct. The architecture is designed to capture the interplay between individual memory, personal motivation, and the overarching social environment. 2.1 Agent Architecture and the Ethical State Space Each agent in the system, defined as a Particle, is characterized by an internal state vector known as the Ethical View. Unlike static models of cooperation, the Ethical View is a dynamic construct that evolves based on the agent’s history of interactions. The core of this state is defined by the parameter μ[k]∈[0,1]μ[k]∈ [0,1 ], representing the Unethical Behavior Probability in kthk^th iteration. At any given interaction in kthk^th iteration, an agent chooses an unethical path with probability μ[k]μ[k] and an ethical path with probability 1−μ[k]1-μ[k]. The evolution of μ[k]μ[k] is not arbitrary but is governed by a secondary internal variable, ρ[k]ρ[k], which represents the Cumulative Motivational Rate of Change. Mathematically, the tendency of an agent to shift toward or away from ethical behavior is dictated by the magnitude and sign of ρ[k]ρ[k]. To ensure that the probability μ[k]μ[k] remains bounded within the logical [0,1] [0,1 ] interval, we employ a logistic-style growth equation for its update: μ[k+1]=μ[k]+Δρ[k]μ[k](1−μ[k])μ[k+1]=μ[k]+ ρ[k]μ[k] (1-μ[k] ) (1) where Δ is the simulation time step. This formulation ensures that the ethical state demonstrates ”saturation” effects near the boundaries, representing the increased resistance to change when an agent is either deeply entrenched in ethical conduct or fully committed to unethical behavior. 2.2 Memory Dynamics: The L-Length Finite Horizon A critical feature of the system architecture is the inclusion of a finite memory for each agent, implemented through a Whitelist and a Blacklist. Each particle maintains these lists with a maximum length L, representing the cognitive or social limits of an individual’s memory regarding past encounters. • Whitelist: Contains identifiers of agents who have previously acted ethically toward the particle. • Blacklist: Contains identifiers of agents who have committed unethical acts toward the particle. As a particle system, this memory structure introduces a ”Sallen-Key” style filtering effect on social interactions: as the population N grows significantly larger than the memory capacity L, the probability of an agent re-encountering a specific individual from their memory diminishes. This structural constraint is the primary mechanism through which population scaling decouples individual ethical sanctions from global social behavior. 2.3 Interaction Logic and Stochastic Feedback Loops The societal dynamics are driven by pairwise interactions between particles, selected through a random shuffling mechanism to simulate the stochastic nature of social encounters. During each interaction between a primary agent (p1p_1) and a secondary agent (p2p_2), the behavior of p1p_1 is determined stochastically by its current unethical probability μ[k]μ[k]. The system distinguishes between two primary categories of social encounters, each governed by specific feedback functions that map the event into a change in the motivational rate, accumulating directly into ρ: 1. Initial Interactions: If p1p_1 is not present in p2p_2’s memory records (residing in neither the whitelist nor the blacklist), the interaction is treated as an unestablished social encounter. The system handles this regime through the following flattened sequence: • Ethical Action: If p1p_1 acts ethically, p2p_2 appends p1p_1 to its localized Whitelist. • Unethical Action: If p1p_1 acts unethically, p2p_2 appends p1p_1 to its localized Blacklist. • Rate Accumulation: The motivational feedback rates for both participating particles are updated concurrently via the functions η (for the primary actor) and γ (for the secondary observer or victim). 2. Re-interactions (Reciprocity): If p1p_1 is already structurally recognized by p2p_2, the encounter bypasses initial learning and triggers a Reciprocal Feedback mechanism governed by the function λ. The direction, sign, and magnitude of this feedback are modeled using a single level of structural possibilities: • Altruistic Gifting (λpos<0 _pos<0): Occurs during a positive re-interaction when p1p_1 acts ethically while present on p2p_2’s whitelist within a cooperative society. Ethical behavior is actively rewarded or ”gifted” by peers, driving down the motivational rate ρ and decreasing future unethical choices. • Ethical Deprivation (λpos>0 _pos>0): Occurs during a positive re-interaction under predatory socio-economic conditions. Here, an agent experiences an ethical deprivation cost, missing out on the material gains of exploitation while witnessing unpunished unethical actors being rewarded, paradoxically accelerating the agent’s transition toward opportunism. • Neutral Indifference (λpos=0 _pos=0): Occurs during a positive re-interaction where routine ethical conduct yields no motivational modification, representing an environment where compliance is taken for granted without generating structural social capital or friction. • Negative Punishment Dynamics (λneg _neg): Occurs during a negative re-interaction when p1p_1 acts unethically and is already recognized on p2p_2’s blacklist. This triggers a strong, stabilizing negative reinforcement value designed to simulate localized social punishments, ostracization, or peer-enforced reputational sanctions. 2.4 Mapping Functions and Social Pressure The internal update functions (η, γ, λ) serve as the core control parameters of the dynamical system. To model the non-linear influence of the collective macroscopic environment on microscopic individual motivation, we utilize a Hyperbolic Tangent-based Switch Function for η and γ: f(x;c,s,a,b)=(b−a)2[1−tanh(s(xc−1))]+af (x;c,s,a,b )= (b-a)2 [1- (s ( xc-1 ) ) ]+a (2) In this formulation, x represents the Mean Unethical Probability (μmean _mean) of the society, c defines the critical threshold center, s dictates the scaling slope of the transition, while a and b represent the lower and upper asymptotic bounds, respectively. This formulation acts as a non-linear regulatory mechanism that tracks environmental feedback: • Sub-threshold Regime (micro-scale or clean societies): When μmean<c _mean<c, the function approaches its upper asymptotic bound b. For an agent committing an unethical act, this yields a positive motivational return accumulating on ρ, indicating that low systemic corruption provides an open, unpoliced window where individual temptation can freely accelerate. • Super-threshold Regime (macro-scale or corrupt societies): As the society passes the critical threshold c, the widespread deterioration of social norms activates environmental resistance or formal systemic pressures. The function drops toward its lower asymptotic bound a (a<0a<0), introducing a negative feedback force on ρ, that suppresses further unconstrained growth of individual unethical probabilities. 2.5 Damping and Decay of Motivation To prevent the system from experiencing infinite motivational growth and to simulate the ”cooling off” period after a social interaction, we introduce a Damping Factor ϕφ. Following each iteration of the society, every agent’s motivational rate ρ is updated as: ρ[k+1]=ϕρ[k]ρ[k+1]=φρ[k] (3) where 0<ϕ<10<φ<1. This ensures that while an interaction can cause a sharp spike in an agent’s tendency to change their behavior, this impulse decays over time unless reinforced by further social contact. 2.6 Global Iteration and Social Shuffling The temporal evolution of the society is modeled as a series of discrete time-steps. In each iteration, a ”Social Shuffling” process occurs to ensure the stochasticity of encounters. The population N is randomly permuted, and a subset of particles, defined by the number of interactions in each iteration, is paired for collision. This mechanism ensures that the interaction density is controlled and that no agent is artificially locked into a fixed neighborhood, mirroring the fluid nature of mass societies. Following the pairing phase, the system performs a global update of all internal states. The calculation of the Global Ethical Metric (μmean _mean) is performed at the end of each step: μmean[k]=1N∑i=1Nμi[k] _mean[k]= 1N _i=1^N _i[k] (4) This value serves as the primary state variable of the society, acting as the global feedback signal that influences individual motivational functions in subsequent steps. 2.7 Monte Carlo Framework and Statistical Convergence To account for the high sensitivity of the system to initial conditions and the stochastic nature of the hyperbolic tangent switch functions, a Monte Carlo (MC) methodology is employed. A single execution of the code (running for 20,00020,000 iterations) represents one realization of the social trajectory. To achieve statistical significance and minimize the impact of ”stochastic outliers,” the simulation is repeated across a large ensemble of independent runs. Each realization begins with: • Heterogeneous Initialization: Initial ethical probabilities (μ) and damping rates (ϕφ) are assigned using a random distribution, ensuring the population is not a monolith but a collection of diverse ethical profiles. • Variable Memory Allocation: The memory length L for each particle is slightly randomized to reflect varying cognitive capacities for social record-keeping. The results from these multiple runs are averaged to compute the Expected Social Behavior. This ensemble averaging allows us to identify the ”Steady State” of the society and observe the Phase Transition from ethical stability to legal decoupling. By sweeping the population parameter N across different Monte Carlo sets, we can precisely quantify the point at which the society’s memory-based ethical enforcement collapses, leaving the global ethical state to drift toward the ”unethical but legal” attractor. 2.8 Sensitivity and Stability Analysis Finally, the Monte Carlo setup is used to perform a sensitivity analysis on the mapping functions η, γ, and λ. By adjusting the slopes (s) and centers (c) of the switch functions, we observe how resilient the society is to ”ethical shocks.” This numerical approach allows us to treat ethics not as a philosophical absolute, but as a self-organizing system whose stability is a function of interaction frequency, individual memory, and population scale. 3 Simulations and Analysis To validate the theoretical framework of moral scalability and the decoupling of law from ethical behavior, we executed a series of numerical experiments utilizing the detailed multi-particle stochastic simulation platform. The simulation parameters, including the specific mappings for the switch functions (η,γ,λη,γ,λ), the damping rate ϕφ, and initialization thresholds, were tuned to explore the phase space of the system. Specifically, the baseline multi-agent ensemble is structured around a heterogeneous population of N particles, partitioned into two distinct sub-populations to capture structural socio-behavioral diversity: Pi∈Group 1 (N1=78N),governed by η1,γ1,λ1Group 2 (N2=18N),governed by η2,γ2,λ2P_i∈ casesGroup 1 (N_1= 78N),&governed by _1, _1, _1\\ Group 2 (N_2= 18N),&governed by _2, _2, _2 cases (5) In which, i is the individual index. The continuous temporal domain is discretized using a uniform simulation time step of Δ=0.1 =0.1, with each unique realization spanning a relaxation horizon of 20,00020,000 global iterations to guarantee asymptotic convergence into a statistical steady state. Individual cognitive storage limits are governed by a base memory length parameter of L, subjected to a randomized scaling allocation modeled as: Li=⌊(0.75+0.25)L⌋L_i= (0.75+0.25U )L (6) where U represents a random number in [0,1] [0,1 ] interval with uniform continuous distribution used to induce structural heterogeneity across individual memory horizons. Particle properties are initialized dynamically using stochastic profiles bounded by maximum operational thresholds for the initial state attributes: μi[0]∼μmax,μmax=0.3ϕi∼ϕmax,ϕmax=0.8 cases _i[0] _maxU,& _max=0.3\\ _i _maxU,& _max=0.8 cases (7) while setting the initial motivational acceleration rate to a uniform rest state of ρi[0]=0 _i[0]=0 for all particles. During each discrete temporal step, a collision density of ⌊N/2⌋ N/2 random pairings is executed via a uniform shuffling index sequence, establishing a high-fidelity stochastic network landscape that prevents localized spatial locking. The core microscopic behavioral updates are dictated by the environmental switch parameters evaluated at the end of each iteration loop. For Group 1, representing the standard sociological segment, the mapping profiles are defined as: η1(ℐ)=0.0,ℐ=ethicalf(μmean;0.5,2.0,−0.2,0.2),ℐ=unethical _1(I)= cases0.0,&I=ethical\\ f( _mean;0.5,2.0,-0.2,0.2),&I=unethical cases (8) γ1(ℐ)=0.0,ℐ=ethicalf(μmean;0.5,2.0,−0.2,0.2),ℐ=unethical _1(I)= cases0.0,&I=ethical\\ f( _mean;0.5,2.0,-0.2,0.2),&I=unethical cases (9) λ1(ℐ)=0.0,ℐ=ethical−0.6,ℐ=unethical _1(I)= cases0.0,&I=ethical\\ -0.6,&I=unethical cases (10) In which, ℐ∈ethical,unethicalI∈ \ethical,\ unethical \ is the interaction type and f(x;c,s,a,b)f(x;c,s,a,b) represents the hyperbolic tangent switch function defined in equation (2). Conversely, Group 2 models a highly specialized behavioral segment characterized by asymmetric sensitivity and tighter behavioral bounds, operating under the following functional mapping criteria: η2(ℐ)=0.0,ℐ=ethicalf(μmean;0.5,2.0,−0.1,0.0),ℐ=unethical _2(I)= cases0.0,&I=ethical\\ f( _mean;0.5,2.0,-0.1,0.0),&I=unethical cases (11) γ2(ℐ)=−0.025,ℐ=ethicalf(μmean;0.5,2.0,−0.05,0.05),ℐ=unethical _2(I)= cases-0.025,&I=ethical\\ f( _mean;0.5,2.0,-0.05,0.05),&I=unethical cases (12) λ2(ℐ)=−0.1,ℐ=ethical−0.6,ℐ=unethical _2(I)= cases-0.1,&I=ethical\\ -0.6,&I=unethical cases (13) 3.1 Parametric Sensitivity and Systemic Bifurcation Analysis To systematically evaluate the structural resilience of the proposed multi-particle framework against severe environmental perturbations, we subjected the mixed population to a coordinated, non-localized behavioral disruption, an ”ethical shock.” Right after initialization, a randomized fraction (≈50%≈ 50\%) of the aggregate particle ensemble had their individual unethical probability dynamically inflated via μi←min(0.95,μi+0.40) _i← (0.95, _i+0.40), artificially simulating a sudden collapse of legal deterrence or an injection of systemic corruption. The capacity of the decentralized moral ecosystem to absorb, damp, or catastrophically cascade this structural trauma was explored across five independent Monte Carlo parameter regimes sweeping the center (c) and slope (s) dimensions of the non-linear hyperbolic tangent regulator f(x;c,s,a,b)f(x;c,s,a,b) defined in equation (2). The ensemble trajectory averages over independent stochastic realizations reveal a stark, structurally defined bifurcation of global steady states (μmean _mean), illustrated comprehensively in Figure 1. Figure 1: Ensemble average trajectory profiles. The emergent asymptotic behaviors can be classified into three distinct structural regimes dictated by the systemic sensitivity thresholds: 1. Hyper-Sensitive Restorative Regimes (c=0.3,s=1.0c=0.3,s=1.0 and c=0.4,s=1.5c=0.4,s=1.5): Under these configurations, the society features an early-onset regulatory trigger. Despite the substantial initial shock magnitude, the global state vector rapidly relaxes downward within the first T≈1,200T≈ 1,200 steps, converging cleanly onto an ethically optimized steady state (μmean≈0.20 _mean≈ 0.20 and μmean≈0.28 _mean≈ 0.28, respectively). Because the center threshold c is low, minor levels of collective deviation are sufficient to trigger strong negative shifts in the motivational acceleration rate (ρ≪0ρ 0) through the microscopic functions η and γ. Individual opportunistic behavior is fiercely penalized during early encounters, driving the network toward complete macro-stabilization before temporary interactions can harden into unpunished habits. 2. Balanced Baseline Regime (c=0.5,s=2.0c=0.5,s=2.0): This parameterization acts as a symmetrical median. The socio-behavioral system displays standard tipping-point characteristics, absorbing part of the shock but showing bounded structural degradation. The ensemble curve settles stably around an intermediate equilibrium of μmean≈0.38 _mean≈ 0.38. The system remains homeostatic, but its capacity to fully restore the pre-shock ethical equilibrium is partially degraded due to the sharper transition boundary. 3. Rigid Insulated Cascading Regimes (c=0.6,s=3.0c=0.6,s=3.0 and c=0.7,s=4.0c=0.7,s=4.0): In direct contrast to the sensitive variants, these settings trigger an irreversible behavioral cascade. The system trajectories spike abruptly to high-entropy, low-trust attractors, stabilizing permanently at μmean≈0.49 _mean≈ 0.49 and μmean≈0.60 _mean≈ 0.60. Because the collective switch center c is set high and the slope s is steep, the peer-enforced moral regulation mechanism remains entirely dormant and ”insulated” during the early stages of widespread decay. Particles freely exploit the absence of local social sanctions. By the time the escalating collective metric crosses the high threshold boundary, opportunistic conduct has already normalized into a dominant, absorbing strategy across the network, locking the society into a permanent legal-moral decoupling phase. Crucially, the narrow standard error bounds and the rapid transition to a bounded, oscillatory statistical steady state for all trajectories after approximately 2,500 iterations demonstrate that long-term societal configuration stability is an explicit function of early-stage feedback sensitivity rather than late-stage rigid resilience. 3.2 Phase Transition and The Decoupling Threshold The first simulation scenario systematically swept the population size N across a geometrically spaced range from 10110^1 to 10310^3 to observe the macro-societal steady state under a fixed base memory capacity L=20L=20. For each population scale, a Monte Carlo ensemble of 2020 independent realizations was executed for 20,00020,000 iterations to ensure asymptotic convergence. Figure 2: Phase transition and the decoupling threshold diagram. As illustrated in the phase diagram, Figure 2, the system exhibits distinct behavioral regimes separated by a critical scaling threshold. In small-scale regimes (N≤50N≤ 50), the steady-state mean unethical probability (μmean _mean) remains suppressed at a baseline level (μmean≈0.33−0.36 _mean≈ 0.33-0.36). Within this domain, individual memory capacity is commensurate with the network size, yielding a high probability of local re-interactions. This local density enables the decentralized feedback loop driven by λ to actively penalize defectors and sustain a highly ethical norm. However, as the population scales past N=50N=50, a prominent non-linear phase transition occurs. The steady-state unethical footprint experiences a sharp, monotonic increase, climbing rapidly toward an upper saturation bound (μmean>0.60 _mean>0.60 at N=1000N=1000). Because the interaction volume is normalized proportionally to N to preserve constant social encounter frequencies per agent, this behavioral decay is isolated entirely as a function of network scale. As N→∞N→∞, the structural limit of individual memory (L) is overwhelmed by the vast space of anonymous agents. The individual retribution mechanism (λ) is effectively diluted to zero as the probability of a secondary encounter vanishes, forcing the system to decouple from its moral foundation and drift toward the minimalist, ”unethical but legal” attractor governed solely by the centralized legal floor. 3.3 Memory Horizon Dynamics To map the interdependencies between cognitive constraints and population scaling, we conducted a bivariate parameter sweep across a matrix of population sizes N∈[10,1000]N∈[10,1000] and base memory limits L∈[2,50]L∈[2,50]. The resulting three-dimensional state space surface is mapped in Figure 3. Figure 3: Memory Horizon Effects on Social Saturated States. The numerical surface demonstrates that the moral state of the particle system is fundamentally governed by the ratio LN LN. When L is small (e.g., L=2L=2), the system is highly sensitive to anonymity, causing the ethical fabric to collapse even at minimal population sizes (N≥20N≥ 20). Conversely, expanding the memory horizon to L=50L=50 significantly broadens the domestic stability domain, shifting the critical phase transition ”elbow” toward larger population thresholds. Crucially, the 3D topology reveals that while increasing individual memory capacity can delay the onset of moral decoupling, it cannot entirely prevent it. As N approaches the macroscale limit (10310^3), all memory contours asymptotically converge toward the maximum saturation zone (μmean≈0.60 _mean≈ 0.60). This bounds the limits of localized, memory-driven social policing, proving that decentralized moral enforcement scales sublinearly against geometric population growth. 3.4 Hysteresis and Path Dependency Analysis To investigate the systemic resilience and irreversibility of social norms, we subjected the society to a cyclic macro-state trajectory, executing an expansion path (Forward Path) immediately followed by a contraction path (Backward Path). The population was scaled up from N=50N=50 to N=1000N=1000, relaxed to steady state at each stage for 20,00020,000 iterations, and subsequently pruned back down to N=50N=50. (a) (b) Figure 4: Irreversibility and Memory of Social Norms. The resulting path-dependent response is captured in Figure 4. The simulation reveals a tight, yet persistent hysteresis loop between the expansion and contraction trajectories. Intriguingly, the directional behavior of this loop exhibits bistability depending on the micro-structural configuration of the active population: • Lagged Recovery (Upper Backward Path): Under standard parameter bounds where the damping rate ϕ=0.8φ=0.8 allows for extended motivational memory, the backward path remains slightly above the forward path. This implies a moral ”inertia” or scarring effect: once an anonymous mass society has normalized unethical conduct, reducing the physical group size does not instantly restore the moral consensus, as the historical inflation of internal unethical probabilities (μi _i) resists rapid realignment. • Accelerated Recovery (Lower Backward Path): In contrast, when the interaction dynamics accentuate aggressive, highly localized penalties or when pruning preferentially eliminates high-μ outliers (simulating legal banishment or migration out of corrupt centers), the backward path loops below the forward path. In this regime, the contraction of scale dramatically magnifies the density of the remaining whitelist networks, allowing the micro-societal feedback mechanism to regenerate structural ethics at an accelerated rate compared to its original assembly. This bistable non-reversible trajectory mathematically formalizes the assertion that macroscale social morality is not a static state function, but a historically contingent dynamical system. 4 Systemic Discussion The numerical results presented in Section 3 offer a quantitative foundation for what we define as the Scalability of Morality. By analyzing the macroscopic states of our particle system, we can evaluate the structural breakdown of decentralized social control through the lens of complex systems and control theory. This discussion contextualizes the critical implications of the ”decoupling effect” and explores the fundamental limits of relying purely on formal legal frameworks to maintain social integrity in expanding populations. 4.1 The Anonymity Phase Transition as an Information-Deterioration Phenomenon From a systems engineering perspective, the transition of μmean _mean from a suppressed, stable state to a high-entropy saturation zone (as observed in Figures 2 and 3) can be modeled as an informational decay within the societal network. In small-scale populations (N≤30N≤ 30), the system operates with high feedback fidelity. Because the population size N is well within the boundaries of the individual memory capacity L, the topological connections formed via the agent blacklists and whitelists remain dense and highly correlated. Every pairwise interaction carries a high probability of structural reciprocity (λ), meaning that an agent’s current behavioral choice directly influences their future state via immediate local feedback. As N→∞N→∞, the interaction topology undergoes a fundamental phase transition. Even though the frequency of encounters per particle remains invariant (since the number of interaction is proportional to N), the probability that any two specific particles collide more than once scales as (L/N)O(L/N), rapidly approaching zero. Consequently, the local feedback loops governed by the reciprocity function λ are structurally neutralized. The ”Anonymity Effect” acts as a high-frequency noise filter that dilutes the historical memory of interactions across the network. In this macro-scale regime, individual particles no longer face localized, memory-driven retribution from their peers. The decentralized, bottom-up regulatory mechanism, which naturally enforces ethical conduct in micro-societies, collapses due to a structural lack of information retention. Ethics, which relies heavily on localized accountability and the preservation of reputation, is thus stripped of its evolutionary feedback mechanism when a system scales past its cognitive storage limits. 4.2 The Structural Inadequacy of the Centralized Legal Floor When decentralized moral feedback is neutralized by population scaling, the burden of preserving systemic stability shifts entirely onto the centralized legal apparatus. In our model, this legal domain is represented by the global baseline dynamics within the mapping functions η and γ. However, as demonstrated by the high saturation values (μmean>0.60 _mean>0.60) at N=1000N=1000, a centralized framework is structurally ill-equipped to compensate for the erosion of decentralized ethics. This systemic failure stems from a fundamental design constraint: legal systems are inherently minimalist and non-adaptive. To prevent tyrannical over-regulation, formal laws are codified to police a restricted subset of critical societal disruptions (gross violations), leaving a vast spectrum of interpersonal conduct completely unmonitored. When the local social pressure governed by the switch functions (η,γη,γ) flattens out due to anonymity, agents quickly discern that they can optimize their individual utility by operating in the ”unethical but legal” domain. Because a centralized legal system lacks the multi-point, fine-grained observation capacity of a decentralized peer-to-peer memory network, it cannot enforce the nuances of moral obligation. The legal framework provides a rigid static boundary rather than a dynamic restoration force. Thus, as a society expands, the progressive decoupling of law and ethics forces the global state toward a minimalist equilibrium, a state where agents are legally compliant but ethically detached, severely undermining the long-term resilience of the collective. 4.3 Hysteretic Memory and the Irreversibility of Social Decay The emergence of a path-dependent hysteresis loop (Figure 4) highlights a vital temporal characteristic of societal dynamics: the state of social morality is non-Markovian. The system’s macro-state is not merely a reflection of its current scale N, but is deeply contingent upon its historical trajectory. When a society undergoes scale expansion (Forward Path), the gradual transition into anonymity inflates the individual unethical probabilities (μi _i) across the population. Once these state variables saturate near their upper bounds, they establish a new, highly resilient behavioral norm. If the society subsequently experiences a contraction in scale (Backward Path), such as urban decentralization or community fragmentation, the system does not simply retrace its steps along the forward curve. The ”moral scarring” represents a systemic inertia where agents, having adapted to an anonymous, low-trust environment, continue to employ defensive or opportunistic strategies even when placed back into a high-density, micro-scale network. This asymmetry proves that rebuilding social trust and moral consensus requires significantly higher energy expenditures (stronger localized incentives, structural re-alignment, or extended relaxation times) than the energy that caused the initial decay. The non-reversible nature of this loop implies that policies aimed at mitigating ethical decline must focus heavily on proactive structural preservation, as treating a highly saturated, decoupled macro-system after a collapse is computationally and socio-structurally inefficient. 5 Conclusion In this study, we have provided a rigorous, particle-based numerical characterization of the ”Scalability of Morality”. By modeling individual agents as discrete particles with finite memory capacity (L) and stochastic choice profiles (μ), we successfully simulated the complex feedback loops that govern the evolution of cooperation and ethical conduct in varying population scales. Our Monte Carlo simulations confirmed a distinct, non-linear phase transition: as the population size N expands significantly beyond the individual memory capacity (N≫LN L), the probability of recurrent local encounters scales as (L/N)O(L/N), effectively approaching zero. This structural dilution neutralizes decentralized retribution mechanisms (λ). Consequently, the enforcement of social norms completely decouples from localized moral networks, forcing the macroscopic system to drift toward the minimalist, low-trust attractor governed solely by the centralized legal floor. Furthermore, our bivariate parameter sweeps mapped out the exact 3D topology of memory horizon dynamics, demonstrating that while expanding individual memory capacities can delay the onset of moral decay, it cannot fundamentally prevent it under geometric population growth. Crucially, the cyclic scale experiments exposed a profound, path-dependent hysteresis loop. This non-Markovian property proves that macroscale social morality is highly contingent upon its historical trajectory; once an anonymous mass society has normalized opportunistic conduct, reducing the physical group size does not instantly restore the structural ethical consensus due to the intrinsic inertia of agent state probabilities. Future directions for large-scale social modeling will focus on the integration of heterogeneous legal frameworks where enforcement probability is also modeled as a dynamic, spatially constrained variable rather than a centralized constant. Additionally, we aim to explore the impact of structured network topologies, such as scale-free or small-world architectures, to determine if specific clustering methods can isolate sub-communities from the global anonymity effect, thereby preserving decentralized moral feedback loops even within macroscale populations. This computational framework serves as an engineering step toward designing more resilient, self-organizing governance structures for modern high-density, interconnected populations. References [1] R. Axelrod and W. D. Hamilton (1981) The evolution of cooperation. science 211 (4489), p. 1390–1396. Cited by: §1. [2] A. Baronchelli (2018) The emergence of consensus: a primer. Royal Society open science 5 (2). Cited by: §1. [3] M. Belfrage, F. Lorig, and P. Davidsson (2024) Simulating change: a systematic literature review of agent-based models for policy-making. In 2024 Annual Modeling and Simulation Conference (ANNSIM), p. 1–13. Cited by: §1. [4] J. Bentham (1789) An introduction to the principles of morals and legislation. T. Payne and Son, London, UK. Note: Reprinted by Oxford University Press, 1996 Cited by: §1, §1. [5] A. B. Carroll (1991) The pyramid of corporate social responsibility: toward the moral management of organizational stakeholders. Business horizons 34 (4), p. 39–48. Cited by: §1. [6] C. Castellano, S. Fortunato, and V. Loreto (2009) Statistical physics of social dynamics. Reviews of modern physics 81 (2), p. 591–646. Cited by: §1. [7] R. Conte, N. Gilbert, G. Bonelli, C. Cioffi-Revilla, G. Deffuant, J. Kertesz, V. Loreto, S. Moat, J. Nadal, A. Sanchez, et al. (2012) Manifesto of computational social science. The European Physical Journal Special Topics 214 (1), p. 325–346. Cited by: §1. [8] P. Contucci, I. Gallo, and G. Menconi (2008) Phase transitions in social sciences: two-population mean field theory. International Journal of Modern Physics B 22 (14), p. 2199–2212. Cited by: §1. [9] R. I. Dunbar (1992) Neocortex size as a constraint on group size in primates. Journal of human evolution 22 (6), p. 469–493. Cited by: §1. [10] R. Eynon and N. Gillani (2025) Computational social science and critical studies of education and technology: an improbable combination?. Learning, Media and Technology, p. 1–15. Cited by: §1. [11] P. Gupta, Q. Zhong, H. Yakura, T. Eisenmann, and I. Rahwan (2025) The role of social learning and collective norm formation in fostering cooperation in llm multi-agent systems. arXiv preprint arXiv:2510.14401. Cited by: §1. [12] T. A. Han, Z. Song, T. Cimpeanu, M. H. Duong, M. Krellner, V. Capraro, and M. Perc (2025) Cooperation versus social welfare. Physics of Life Reviews. Cited by: §1. [13] H. L. A. Hart (1961) The concept of law. Oxford University Press. Cited by: §1. [14] R. Jamali and S. Lazarova-Molnar (2026) Agent-based modeling and simulation for economic markets: a comprehensive review of applications, challenges, and opportunities. Journal of Simulation, p. 1–50. Cited by: §1. [15] M. Jusup, P. Holme, K. Kanazawa, M. Takayasu, I. Romić, Z. Wang, S. Geček, T. Lipić, B. Podobnik, L. Wang, et al. (2022) Social physics. Physics Reports 948, p. 1–148. Cited by: §1, §1. [16] I. Kant (2017) The metaphysics of morals. Cambridge University Press. External Links: Document, ISBN 9781107451353 Cited by: §1. [17] E. Koid, H. Van Der Voort, and M. Warnier (2025) Agent-based modeling for data-driven enforcement: combining empirical data with behavioral theory for scenario-based analysis of inspections. Data & Policy 7, p. e6. Cited by: §1. [18] D. M. Lazer, A. Pentland, D. J. Watts, S. Aral, S. Athey, N. Contractor, D. Freelon, S. Gonzalez-Bailon, G. King, H. Margetts, et al. (2020) Computational social science: obstacles and opportunities. Science 369 (6507), p. 1060–1062. Cited by: §1. [19] S. Li, L. Jin, and S. Li (2022) Modeling and analysis of competitive behavior in social systems. IEEE Transactions on Computational Social Systems 10 (3), p. 1347–1355. Cited by: §1. [20] S. N. Menon, V. Sasidevan, and S. Sinha (2018) Emergence of cooperation as a non-equilibrium transition in noisy spatial games. Frontiers in Physics 6, p. 34. Cited by: §1. [21] M. A. Nowak (2006) Five rules for the evolution of cooperation. science 314 (5805), p. 1560–1563. Cited by: §1. [22] E. Ostrom (1990) Governing the commons: the evolution of institutions for collective action. Cambridge University Pres. Cited by: §1. [23] M. Perc, J. J. Jordan, D. G. Rand, Z. Wang, S. Boccaletti, and A. Szolnoki (2017) Statistical physics of human cooperation. Physics Reports 687, p. 1–51. Cited by: §1, §1. [24] M. Pereda, V. Capraro, and A. Sánchez (2019) Group size effects and critical mass in public goods games. Scientific reports 9 (1), p. 5503. Cited by: §1. [25] K. Pistor (2019) The code of capital: how the law creates wealth and inequality. Princeton University Press. Cited by: §1. [26] D. Sarker, Y. Zhang, L. K. Perry, D. S. Messinger, and C. Song (2026) Alignment phase transition in socially driven motion. Science advances 12 (12), p. eadz6890. Cited by: §1, §1. [27] C. D. Stone (1975) Where the law ends: the social control of corporate behavior. Harper & Row, New York, NY. Cited by: §1. [28] A. Traulsen, C. Hauert, H. De Silva, M. A. Nowak, and K. Sigmund (2009) Exploration dynamics in evolutionary games. Proceedings of the National Academy of Sciences 106 (3), p. 709–712. Cited by: §1. [29] R. Vasellini et al. (2025) Applications of agent based models to complex social systems: coordinated human behaviour from the bottom up. Cited by: §1. [30] T. Wang, Y. Zhu, P. Ye, W. Gong, H. Lu, H. Mo, and F. Wang (2022) A new perspective for computational social systems: fuzzy modeling and reasoning for social computing in cpss. IEEE Transactions on Computational Social Systems 11 (1), p. 101–116. Cited by: §1.